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SYNOPSIS 


Unevenness in the track induces vibrations in moving vehicles. 
Aircraft ground maneuvers over uneven runway induce dynamics in 
the air frame. Excessive level of vibration causes discomfort 
among passengers and crews, < produces undesirable movement of the 
cargo and may undermine the vehicle instrumentations and control 
mechanism. The dynamics may lead to fatigue and failure of 
structural components . The pavement also experiences dynamic 
loading, undergoes deflection and stress cycling and the 
associated wear and tear. 

Rapid growth of traffic and high speed transportation has 
brought a significant change in the design of vehicles and their 
component parts. To allow for more passengers and pay loads, 
vehicle size has increased. More advanced materials and 
construction techniques have evolved causing greater flexibility 
in the vehicles. Modern aircraft with its less rigid structures 
and higher takeoff and landing speeds presents the best example. 
Under these conditions, analysis of vehicle dynamic behaviour has 
assumed greater importance for assessing ride quality, design of 
vehicle structures, control of vibration and optimization of 
suspension systems . 

Moving vehicles exert fluctuating forces on the track pavement 
caused by its vibration. The track may undergo deformation of 
varying degree depending on the type of its construction. This is 



particularly important for ropeways, other types of guideways and 
aircraft operations from carrier decks. 

In the present study, analytical methods have been developed to 
study the vehicle-track system dynamics induced by track 
unevenness . The track unevenness is assumed to be a nonhomogeneous 
random process specified by a generalised power spectral density 
function. The mean track profile may be variable to model defects 
like slope, bumps, dips, uneven and differential settlements, 
craters, repair patches, etc. The track is assumed to deflect 
vertically under vehicle imposed loads. It is modelled as a beam 
of uniform mass and flexural rigidity supported by elastic and 
dissipative subgrade. 

Aircraft takeoff is an accelerating while landing is a 
decelerating run. Variable ground motion pattern of the vehicle 
has been modelled by a power series with finite terms. 
Coefficients of the polynomials can be selected to depict any 
ground run pattern of the vehicle. Consequent to the variable 
vehicle velocity and nonhomogeneous characteristics of the track, 
input sensed by the vehicle wheels become nonstationary in time. 

The vehicle has been modelled as a combination of discrete and 
continuous elements with different wheel arrangements. Vehicle 
wheels, tyres and part of the suspension systems are lumped as 
unsprung masses. Aircraft fuselage and wings are modelled as beam 
elements with variable mass, stiffness and damping distributions. 
The wings may carry arbitrary numbers of concentrated masses at 
arbitrary locations. 
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The variation of the sectional properties of the nonuniform 
beam has been modelled by power series . This representation is 
general and any type of variation can be accommodated with desired 
accuracy. A method has been outlined to find the exact solution 
for nonuniform beam vibration. Natural frequencies and mode shape 
functions have been obtained in closed form for bending, torsion 
and bending- torsion coupled dynamics of variable section beams. 

The proposed theory has been validated by comparison with 
experimental values and available published results. Experiments 
have been conducted to obtain resonant frequencies in bending for 
mass loaded tapered beams with different support conditions. 

Coupled bending- torsion free vibration of a wing of a small 
military aircraft has been studied. Available data for wing 
sectional properties at discrete span locations have been fitted 
to the polynomial series and natural frequencies are obtained by 
the present theoretical approach. Available experimental results 
on the natural frequencies of the same wing are compared with the 
theoretical results. 

Vehicle-track dynamic equations have been developed. The 
continuous elements of the model involve partial differential 
equations with space variable coefficients. Equations of motion 
corresponding to lumped elements are ordinary differential 
equations. In the present approach, the equations of motion of the 
continuous systems are first reduced to a set of linear ordinary 
differential equations in time dependent generalised coordinates 
using normal mode expansion techniques. The set of ordinary 





differential equations together with the equations of lumped 
masses are then recast into first order state space form. These 
equations are completely uncoupled using modal analysis. 

Analytical solution for the response has been obtained and 
expressions for second order statistics have been developed. The 
following models, in order of increasing complexity, have been 
considered: 

(1) Single Point Input Heave Model 

(2) Two Points Input Heave-Pitch Model 

(a) Rigid Fuselage. 

(b) Flexible Fuselage in Bending. 

(3) Three Point Input Heave-Pitch-Roll Model 

(a) Flexible Fuselage in Bending and Rigid Roll, 

(b) Flexible Fuselage in Bending and Torsion. 

In the next part of the study, nonlinear behaviour of the 
suspension parameters is included in the analysis. Combination of 
air spring, orifice and coulomb damping have been treated with 
single point input heave model including tapered flexible wings. 
An instantaneous linear approach has been proposed. The method 
derives the benefit of linear theory while at the same time 
retains the original nonlinear behaviour. 

System equations for the nonlinear model have been developed 
wherein the suspensions forces are nonlinear functions of 
instantaneous shock strut stroke and velocity. The system 
equations are integrated numerically with simulated track profile 
samples. Suspension forces are determined at discrete time steps 
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from known response of original nonlinear systems . Derivatives of 
the forces evaluated at discrete time steps are utilized to 
represent ''instantaneous spring and damping constant' ' . With 
these instantaneous suspension parameters, step by step 
application of the analogous linear model has been carried out to 
obtain second order response statistics. Reliability of the 
solution has been assessed by comparison with Monte Carlo 
simulated results. 

A passenger aircraft has been selected in the study for 
generating numerical results. Response statistics have been 
obtained for constant velocity taxi, accelerating takeoff and 
decelerating landing runs. Some major conclusions drawn from the 
studies are 

(1) Mean response are primarily dependent on the track mean 
profile while response covariance are primarily influenced by 
track second order statistics. 

(2) Vehicle and track response are sensitive to the vehicle 
forward speed. For the vehicle, track and subgrade parameters 
selected in the study, the evaluated track deformation is 
smaller in magnitude. 

(3) Unsymmetrical distribution of roughness across the track 
causes a significant system response not only in the roll mode 
but also in the heave and pitch modes due to strong coupling 
in the system. 

(4) Response magnitude of aircraft, in general, is found to be 
more at landing touch down compared to taxi and takeoff runs . 



xxx 

(5) The modelling predicts higher response magnitudes for rigid 
vehicles as compared to flexible vehicles. 

(6) Bending flexibility of aircraft type of vehicle is important 
for accurate description of the response behaviour. The low 
tread makes the torsional flexibility less significant for 
the prediction of vehicle response. 

(7) Consideration of rigid vehicle model may be adequate for 
preliminary design of undercarriages and vehicle structural 
components . 

(8) General response behaviours in instantaneously linearised 
model and Monte Carlo simulation show close similarity. 

(9) The proposed instantaneously linearised scheme uses much less 
computer processing time and storage space as compared to 
Monte Carlo simulation. 

The analytical methods presented can handle the coupled 
vehicle- track dynamics problem with any general description of 
vehicle forward motion. Spatial nonhomogenity of the track 
roughness can be included in the analysis . The method uses linear 
suspension and track behaviour. Flexible vehicle body along with 
slender attachments can be incorporated in the system model . A 
general variation of cross sectional properties of flexible 
systems modelled as a beam can be handled by the analytical 
approach developed for the variable section beam. The results of 
the nonuniform beam vibrations are in good agreement with the 
experimental and other published results. 

A general combination of spring and damping nonlinearity in 


Xxxi 


vehicle suspension can be tackled by the proposed instantaneously 
linear approach. The present approach differs from other 
linearised schemes in the way that instead of finding equivalent 
set of linear parameters at one time, it introduces 'instantaneous 
spring and damping coefficients' evaluated from the response of 
the original nonlinear system at discrete time instants close to 
each other. The proposed instant to instant linearisation scheme 
may be useful to other nonlinear systems where conventional 
methods are difficult to apply. 



CHAPTER 1 


INTRODUCTION 


1.1 TRACK INDUCED VEHICLE DYNAMICS 

All tracked vehicles in motion receive dynamic excitation 
from the track and are induced to oscillate. The primary cause of 
vehicle excitation is the surface unevenness of the track. Track 
surfaces have some irregularities however carefully they are 
prepared. There are several causes of track unevenness*, traffic 
loading, environmental effects, construction materials and built 
in construction irregularities . Newly constructed tracks may be 
poorly finished or may have many design features such as 
construction joints, thermal expansion joints etc. Tracks that 
have been in service for sometime often form localised distresses 
due to the application of heavy wheel load and environmental 
changes. These normally include cracks, bumps, potholes, and other 
defects. Tracks are supported by prepared sub grade /foundation 
system. Poor drainage, swelling soils, freeze -thaw cycles and 
nonuniform consolidation of subgrade may all contribute to surface 
roughness . 

Vehicles are complex structural systems with combinations of 
wheels and body connected by suspension elements . The shock and 
vibration caused by track unevenness are transmitted from wheels 
to the superstructure of the vehicle through the suspension 
system. The vertical irregularities due to variations in the track 
level in longitudinal and transverse directions induce in the 
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vehicle heave, pitch, roll, yaw, side sway and axial motions. The 
result is a time varying stress in the structural components. 
These may, at times, reach levels to cause structural failures. 
However, in general, stress levels are much lower than the 
critical value. The cumulative effects of the repeated stress 
cycling may be significant to bring about metal fatigue damage. 
The excessive level of vibration causes discomfort among 
passengers and crews, produce undesirable movement of the cargo 
and may undermine the vehicle instrumentations and the control 
mechanism. 

Suspension system has a major role in absorbing the shock and 
vibration effects on the vehicle. Different types of shock 
absorbers are used in the suspension systems to meet the 
requirements. Most ground vehicles and aircraft landing gear 
suspensions use oleopneumatic type absorbers. The oleopneumatic 
absorber, commonly termed as oleostrut, is a member composed of 
two telescopic cylinders. When the strut compresses or expands, 
oil inside the air tight cylinder is forced through an orifice 
from one cylinder to the other. This provides a damping effect 
while compression and expansion of air inside the air chamber 
produces spring action. 

A moving vehicle exerts a fluctuating component of force on 
the track due to vibratory motion induced by the surface 
unevenness. For a given track, the magnitude of the wheel-track 
force variation is a function of the track roughness, vehicle 
speed and suspension characteristics. The fluctuating nature of 
the dynamic wheel load, with increase in number of load 
repetitions, may cause track degradation -and fatigue damage. 
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For guided ground vehicles, the track is formed with slender 
members while for non guided vehicles the track is generally a 
paved surface. The track is supported by proper foundation systems 
whose construction techniques are different for different types of 
tracks. Nonguided vehicles such as automobiles and airplanes use 
roads and airport pavements in their ground maneuvers . These 
pavements are constructed with one or more layers of processed or 
unprocessed materials placed on the prepared subgrade. There are 
two categories of road and airport pavements- flexible and rigid. 
Flexible pavements typically consist of bituminous surface course, 
a base course of crushed or uncrushed aggregate and subbase 
course of lower quality material than in base course. Rigid 
pavements consist a of slab of cement concrete that rests on a 
prepared subgrade or subbase. This is reinforced with steel to 
control and minimise harmful effects of cracking. Concrete 
pavements are provided with longitudinal and transverse joints 
between the slabs. 

Rapid growth of traffic and high speed transportation has 
enhanced the importance of track induced vehicle dynamics. 
Analysis of vehicle dynamic behaviour is essential for assessing 
ride quality, design of vehicle structures, control of vibration 
and optimisation of suspension parameters. This has brought about 
improved vehicle-track design, lower maintenance costs and greater 
ride comfort . 

1.2 LITERATURE REVIEW 

Study of track induced dynamics of vehicle has assumed 
greater significance from the view point of passengers comfort, 
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cargo safety, fatigue damage analysis and controllability of the 
system. A review of earlier works on different aspects of the 
problem is being presented in this section. 

In early attempts to investigate vehicle dynamics, excitation 
from the ground has been modelled as sinusoidal waves, triangular 
waves or step functions. Response have been obtained for a single 
degree freedom vehicle model with linear suspension at constant 
velocity run with step input by Houbolt [1] and with sinusoidal 
input by Wong [2] . 

While these studies could provide a basis for various 
designs, they could not reflect the actual ride behaviour since 
surface profiles are rarely of simple forms. Soon it was realised 
that ground surface profiles should be more realistically 
described as a random process. The concept of the track profile 
being a random process, has induced many studies for the 
characterisation of input process in probabilistic sense. Walls et 
al. [3] carried out spectral analysis from the measured data of 
actual runway surface. Power Spectral Density (PSD) function of 
runway elevations have been determined in order to permit a 
description of the frequency characteristics of the runway 
roughness . 

Dodds and Robson [4] have shown that typical road surfaces 
may be considered as realisation of homogeneous and isotropic two 
dimensional Gaussian random process . The complete description of 
the roughness is provided by single autocorrelation function 
evaluated from any longitudinal profile. 

Different track roughness models have been proposed by. 
Sussman [5], Kamesh and Robson [6] and. Honda et al . [7]. Their 
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relationships cater for a wider frequency range. It has been shown 
that PSD is related to parameters called roughness and waviness 
index. These parameters can shift the roughness level up and down 
and thus a guideline for classification of roads in terms of 
roughness has been framed. 

Heath [8] developed cross-spectrum between heights of 
parallel profiles on a homogeneous and isotropic random surface. 
Formulae involving single integration have been derived which 
express the cross -spectrum in terms of PSD of single track. Limit 
behaviour of the cross -spectrum at small and large wave numbers 
have been found and discussed. 

Marcondes et al. [9] have shown that some modification of PSD 
of track profiles are necessary when road section contains 
concentrated roughness like pothole, discrete bump, rail road 
crossing and other inhomogeneities. From profile survey report of 
several road and airport pavements, they suggested different fits 
for different section of the pavement incorporating a 
discontinuity in frequency. 

Later Marcondes et al . [10] proposed a model of relationship 

between pavement PSD and parameter called International Roughness 
Index (IRI) . In order to establish the relationship, spectral 
analysis of different measured profiles was carried out whereas 
IRI was computed using mathematical model of quarter car 
simulation. 

A recent study by Iyengar and Jaiswal [11] has shown that 
train track irregularity can be- modelled as stationary Gaussian 
process which can be defined in terms of the PSD function. PSD. 
function obtained can be used to extract further information on 
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level crossing and peak statistics of unevenness data. 

Stochastic description of track input in terms of mean and 
correlation or PSD function expedited the research on vehicle 
dynamics using random vibration theory. Silsby [12] studied the 
effect of various landing gear parameters on the airplane response 
on runway roughness with a rigid mass vehicle model. Tung et al . 
[13] used a heave-pitch model of aircraft to study the response of 
aircraft to runway roughness while including landing gear 
nonlinear behaviour. Numerical integration scheme has been 
employed for deterministic input while statistical linearisation 
and perturbation method was used for stationary input. 

Kirk and Perry [14] studied the taxiing induced vibration in 
aircraft by the PSD method with two degrees of freedom heave 
model. This work has been extended to heave-pitch degrees of 
freedom model [15] . The aircraft fuselage has been idealised as a 
rigid beam. Suspension parameters have been linearised by 
equivalent linearisation technique. 

Prediction of vehicle response at constant velocity run is 
not always appropriate as it can be appreciated that vehicle 
velocity is commonly variable in time. Specifically to aircraft, 
landing and takeoff run take place in decelerating and 
accelerating conditions. This renders the track input experienced 
by the vehicle nonstationary. 

Virchis and Robson [16] conducted an analytical study on a 
simple vehicle model at accelerating runs . They concluded that 
nonstationary effects are not large and can be neglected as a 
first approximation. Sobczyk and Macvean [17] developed a method, 
for analysing the response of linear ■ vehicle travelling with 
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variable velocity. Both deterministic and random variation of the 
vehicle forward velocity have been considered. Numerical 
integration of system equations in time domain with track 
correlation as input yields the second order response statistics . 

Yadav and Nigam [18] investigated the response of lumped mass 
model of vehicle to nonstationary ground input using a formulation 
of the system equation in space domain. This results in the 
governing differential equations to have space dependent 
coefficients while the input process remains homogeneous in space 
domain. The evolutionary spectral approach is then used to 
determine the response statistics of vehicles in constant and 
variable velocity runs . 

More complicated ground velocity variations are incorporated 
in the study of vehicle dynamics with the concept of filter theory 
in order to simplify stochastic ground inputs. Kotb et al . [19] 

used a frequency domain formulation with modelling of the spectral 
density function as the output of a second order shaping filter 
under white noise input. Hammond and Harrison [20] have employed 
time domain solution for single degree of freedom linear vehicle 
model. State space form has been employed to represent the vehicle 
dynamic equations . Track excitation utilised a spatial shaping 
filter model . 

In another study by Harrison and Hammond [21] , a nonlinear 
model excited by nonstationary random process has been considered. 
The method is based on the technique of modelling the input 
process as a spatial shaping filter which may be linked to vehicle 
dynamic equations through velocity functions. The nonlinear, 
problem was overcome by statistical linearisation. 
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Ground profiles are usually nonhomogeneous in practice . This 
non homogeneity may be due to variable mean in presence of slope 
or defects like dip, bumps, fault or uneven settlement of subgrade 
and repair patches etc., even when the random part of the track 
unevenness is homogeneous in space. Yadav [22] obtained the 
response of two degrees of freedom heave vehicle model for 
nonhomogeneous track input by direct solution of the system 
equations. Yadav and Kapadia [23] and Yadav and Ramamoorthy [24] 
have used Monte Carlo simulation technique for the evaluation of 
response statistics of heave and heave-pitch models of aircraft. 
The effect of linkage dynamics and nonlinearity of the suspension 
characteristics have been considered. Track inputs were simulated 
in time domain with known statistics of input and ground velocity 
history. Response has been found by numerically integrating the 
system equations. Averaging across a large number of response 
samples yields the vehicle response statistics. 

Vehicle suspension system may be passive with spring and 
dashpot elements which do not require any external power source to 
operate. On the other way, active system generates suspension 
forces using servomechanism with external power source. Dahelberg 
[25] studied the response of a single degree of freedom vehicle 
model and developed an optimisation procedure to compare the 
behaviour under passive and active systems. The suspension 
parameters have been optimised under the constraint of ride 
comfort, road holding and suspension working space. The scheme has 
been extended to two degrees of freedom system [26] . The random 
road input was considered as zero mean random process. 

Narayanan and Raju developed stochastic optimal control 
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theory of nonstationary response of a single degree freedom [2 7] 
and two degree freedom vehicle model [2 8] . Vehicle has been 
idealised as rigid lumped mass/masses, suspension behaviour being 
linear. Effect of nonlinearity on active system performance has 
been examined in extension of the work [2 9] . Velocity squared 
quadratic damping and hysteresis type of stiffness nonlinearities 
have been considered in the analysis. Constant acceleration run 
and random road profiles have been assumed. 

Aircraft landing response with decelerating run on the uneven 
runway depends to a great extent on the geometric layout of the 
shock strut. Venketeshan [30] carried out a comparative study of 
the landing response of aircraft employing telescopic gear with 
linear and nonlinear characteristics. Relative merits of landing 
gears have been shown comparing their merits under optimum 
conditions. Study has been carried out with lumped aircraft model 
with zero aircraft forward velocity. It was found that performance 
of nonlinear gear is better than the linear one. 

Reddy et al . have analysed levered [31] and semilevered [32] 
suspension type landing gears by including kinematics of the 
linkage in the formulation. They have examined two degrees of 
freedom heave model while incorporating shock strut and tyre 
nonlinearities . 

Yadav and Ramamoorthy [33] have studied nonlinear landing 
gear behaviour at touchdown. A heave pitch model of aircraft with 
oleopneumatic shock strut in the articulated nose and the 
telescopic main gear has been considered. The system equations 
have been integrated numerically to obtain response in the . 
duration of impact. Sensitivity of the system response to 
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variations in some shock strut parameters has been investigated to 
bring out improvements in the performance. 

Yadav and Singh [34] have studied the braked landing response 
of aircraft . Optimal landing run without any skidding has been 
obtained with a one-step prediction of variable braking force with 
a linear two degrees of freedom aircraft model. Performance has 
been studied with different mean shapes for the random runway 
profile . 

Increasing traffic demand and operational speed of ground 
vehicles have brought a significant change in the design of 
vehicle. The vehicles have increased in size and load carrying 
capacity while at the same time becoming lighter and flexible due 
to use of advanced materials in their construction. Effects of 
structural flexibility necessitated the change in modelling 
considerations . 

Flexible models of guided ground vehicles have been analysed 
for response by Wilson and Biggers [35] . In their analysis, the 
vehicle model was considered to be a slender beam of uniform mass 
and stiffness property. Interaction of the first few flexural 
modes of the vehicle along with the rigid body modes have been 
considered. Track irregularity has been modelled as a stationary 
random process. 

Cherchas [36] studied the combined effects of vehicle and 
guideway flexibility in a mathematical model of magnetically 
levitated guided ground vehicle. The study also included the 
effects of vehicle suspension and propulsion pad translation and 
multiple guideway spans. He concluded that vehicle flexibility is. 
significant in certain frequency ranges.. 
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Hac [37] considered a discrete - continuous vibratory system 
at constant forward speed. The model used was to represent a 
model of vehicle with active suspension. The road irregularity has 
been assumed to be a stationary Gaussian zero mean process. He 
evaluated the effect of body elasticity on system dynamics and 
proposed a stochastic optimal control theory for vibration 
isolation in the system. 

Hac et al . [38] modified the control strategies for active 

vehicle suspensions with the goal of reducing high frequency 
vibration. Vehicle elasticity has been taken into account in the 
modelling and vehicle forward velocity has been assumed constant. 
The authors also considered the vibration control problem of 
flexible model with semiactive suspension [3 9] . They used 
adjustable dampers as force actuator in which the rate of energy 
dissipation can be varied either between two distinct levels or 
continuously within a given range of damping. Flexible vehicle 
model has been idealised as uniform beam with few flexural modes. 

In the study of flexible vehicles, many authors have modelled 
the vehicle body as uniform section beam for which modal 
parameters can be easily derived. However, in many practical 
situations, beam geometry varies along the span. This becomes 
sometimes necessary to cut the cost of the material or for 
structural and architectural reasons. The dynamics of nonuniform 
beam element is complex as the governing differential equation 
contains space dependent coefficients. Analytical solutions of 
nonuniform bending vibrations have been attempted by Conway and 
Dubil [40] , Sanger [41] and Gorman [42] with the help of Bessel, 
functions . Natural frequency and mode shapes have been obtained 
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for wedge and cone beams and linearly tapered beams with simple 
boundary conditions. Bessel function solutions for the normal mode 
function are restricted to a class of tapered beam and thus can 
not handle general situation. 

Wang [43] introduced hypergeometric series solution for 
nonuniform beam vibration. The method imposes less restriction on 
the type of taper but is too complicated for general use. 
Naguleswaran has applied Frobenious method for a beam with 
constant depth and linearly varying breadth [44] and wedge and 
cone section beams [45] . Natural frequencies for different support 
conditions and wide range of taper ratios have been obtained. 

Recently Abrate [46] has shown that for some nonuniform 
beams, differential equation of motion can be transformed into the 
equation of motion of a uniform beam. The eigen values then become 
same as for the uniform beam or rod. However, the method has 
serious restriction on the beam's taper and boundary conditions. 

A uniform section cantilever beam with tip mass has been 
analysed for free and forced vibration response by To [47] . Closed 
form solution for the response has been obtained for deterministic 
base excitation. 

The many analytical methods available are unable to handle 
the general variation of cross sectional property in a complex 
structural model of vehicle-*- Chehil and Jategaonkar [48] applied 
Galerkin technique for the solution of natural frequency of beam 
with variable section. 

Bapat et al. [49] have applied transfer matrix method to 
determine natural frequencies of uniform beams carrying several 
concentrated masses. Classical and. nonclassical boundary 
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conditions have been considered in the study. Kim and Dickinsc 
[50] employed Rayleigh Ritz method for the solution of nonunifoi 
beam with linear and parabolic taper. An orthogonal polynomia 
shape function has been assumed. Effect of concentrated masses ar 
different boundary conditions on the natural frequency have bee 
studied. 

Karabolis and Beskos [51] introduced finite element metho 
(FEM) in the static and dynamic analysis of tapered beams. Besko 
also applied the boundary element method [52] to the problem. Us 
of transfer matrix has been made by Chu and Pilkey [53] for th 
solution of nonuniform beam vibration. Wu and Lin [54] suggested 
combined numerical and analytical method for free vibratio 
solution of a beam carrying several concentrated masses. However 
only uniform section beam has been treated. 

Variation of cross sectional properties of a vehicle body ha 
been considered in its modelling by Gibson [55] . He used FEM fo 
the flat car body idealising it as an assemblage of uniform beai 
elements connected at the nodes . Response of the vehicle has bee 
obtained in sinusoidal environment to establish resonano 
frequency and decay rate. 

Khulief and Sun [56] presented an analysis of long wheel basi 
vehicles such as buses, rail cars and truck trailers. Vehicle'! 
elastic chassis has been represented by finite beam elements. Th< 
track irregularity has been modelled by discrete bumps am 
constant vehicle speed has been assumed. Semiactive contro 
technique has been used for the suppression of vehicle vibration 
FEM technique used in flexible vehicle vibration problem i: 
essentially a lumped parameter approach,- wherein the distribute! 
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physical properties of a structure are represented by a finite 
number of idealised substructures interconnected at finite number 
of grid points . Although the method can handle most of the 
practical cases, it requires refined discretization of vehicle 
domain for accurate results leading to a costly computation. 

Apetaur and Opicka [57] developed an analytical method for 
the study of the vehicle response. Vehicle body has been treated 
as a combination of rigid and elastic bodies . Nonlinearity in the 
suspension behaviour has been taken into account . Road input was 
modelled by stationary random process. Results indicated that 
structural flexibility may play an important role in the 
deformation characteristics of the system above certain frequency 
range . 

Valk and Pacejka [58] proposed a mathematical model of 
landing gear taking into account the lateral bending and torsional 
motion of the shock strut. They neglected the randomness in runway 
profiles. The model was applied to test the stability of the 
landing gears in a civil aircraft. 

Vibration of most of the vehicle and their components are 
characterised by a combination of bending deflections and 
torsional rotations. The motions usually remain coupled whenever 
shear center and mass center do not coincide. Aircraft wings, 
control surfaces and many sections in ships and other ground 
vehicles fall under this category. The practical importance of the 
coupled bending- tors ion vibrations of the beam element has induced 
many approaches for evaluation of the normal modes. Gonald [59], 
Timoshenko [60] , and Dokumaci [61] have obtained exact modal 
solutions for uniform section beams with simple boundary 
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conditions. The methods have been applied to thin walled beams, 
aircraft wings and turbine blades. 

Bishop and Price [62] developed an exact solution of 
coupled-bending torsion of a beam including the effect of rotary- 
inertia and shear deformation. Orthogonality condition for the 
coupled normal modes have been derived and used to express the 
equation of motion in terms of normal coordinates. Bishop [63] 
extended the study on coupled bending - torsion of beam including 
warping . 

Different numerical and approximate methods have been tried 
for the solution of coupled bending-torsion modes of beams. 
Anderson and Houbolt [64] used Rayleigh's energy principle for the 
vibration of aircraft wings. Both unswept and swept wing's natural 
frequency and mode shape have been determined. Rao and Carnegie 
[65] used Galerkin technique for bladeing vibration in turbine. 
Mei [66] and Klausbruckner et al . [67] introduced FEM technique 

for coupled bending torsion vibration of beams. Klaushruckner in 
his formulation has taken the effect of warping. The theoretical 
results have been compared with the experimental results for 
channel section beams. Hallaur and Liu [68] and Banerjee [69] have 
developed dynamic stiffness matrix for coupled bending- torsion of 
beam elements. Dynamic stiffness matrix has been obtained by 

directly solving the governing differential equations. The theory 

/ 

developed for an individual member has been applied to multiple 
span beams and more general planar assemblage of beams . 

Eslimy-Isfahany et al. [70] developed an analytical theory 
for the response of a cantilever aircraft wing coupled in bending 
and torsion. Normal mode theory has been applied for the response 
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solution. Both deterministic and stationary random excitations 
have been compared. 

Moving vehicle exerts fluctuating force on the track which 
vary both in space and time. The structure undergoes deformation 
of varying degree depending on its construction. Track design is 
identified as a critical area of potential capital savings in any 
transport system. As track becomes lighter and more flexible, 
material cost is reduced though it causes complex vehicle- track 
interaction and affect ride quality. 

Blejwas [71] developed a procedure for simulating the dynamic 
interaction between traversing vehicle and support structure. 
Lagrange multipliers have been used to satisfy constraint 
relations between the interacting systems . Response has been 
obtained with a point mass and spring-mass-damper system over 
smooth and sinusoidal track at constant vehicle speed. 

Duff [72] considered the vibration analysis of infinite rail 
road track due to moving vibratory load passing over it. The rail 
road track has been modelled as infinite beam resting on Winkler 
foundation. Laplace and Fourier transform have been used to obtain 
the response. 

Yadav and Upadhayay [73] carried out analytical study on 
nonstationary dynamics of rigid heave model of train and flexible 
track over inertial foundation during variable velocity runs. 
Track has been assumed to be a finite length beam resting on 
elastic foundation. Track mean line was considered variable in 
vertical plane superimposed by random roughness. The above model 
was then extended to heave-pitch, heave-roll [74] and heave 
-pitch-roll models [75] . Comparison of model behaviour has 
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revealed strong influence of pitch and roll modes on the bounce 
motion of the vehicles while track response was sensitive to the 
types of sleepers used. 

Hardy and Cebon [76] have analysed the response of continuous 
road pavement to moving vehicle load. Response have been obtained 
in time domain for linear two degrees of freedom heave model using 
convolution theory . Road has been idealised as slender beam 
resting on Winkler foundation having distributed spring and 
damping. Road unevenness has been assumed to have a deterministic 
form represented by discrete steps. 

Hiltumen and Bush [77] developed a testing and analysis 
technique to determine airfield pavement deflection under moving 
aircraft loads. Deflections have been obtained both on the surface 
and at depths below the surface. Geophones were employed in actual 
field measurements to gather information of the propagated waves 
caused by dynamic effects of moving aircraft loads. 

Anant et al . [78] presented a dynamic analysis of rigid 
airport pavement with discontinuous joints. Aircraft model was 
assumed to be single degree of freedom mass, spring and dashpot 
system. Runway pavement has been idealised as thin elastic plates 
resting on a Winkler foundation. FEM technique was used to obtain 
vehicle and pavement dynamic behaviour. Peak dynamic deflection of 
the pavement has been found to depend significantly on the mass 
inertia of the aircraft. 

Cai et al. [79] studied the dynamic interaction between the 
vehicle and guideways of high speed transportation system. 
Flexible guideway has been modelled as a Bernouli Euler beam and 
the dynamic model of the vehicle consisted of a single mass 
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yielding a two degrees of freedom heave model with primary and 
secondary linear suspensions . Modal analysis technique has been 
used for the response calculation. It was shown that for a certain 
range of vehicle-guideway parameters, the interaction between 
guideway and vehicle can be ignored. In this, situation 
representation of vehicle model by single moving mass has been 
found to be adequate. 

In another study, Cai et al . [80] considered the model of 
multi-car and multi-load vehicle travelling at constant speed over 
randomly profiled double span guideway. Coupled effects of 
vehicle-guideway interactions over a wide range of vehicle speed 
and various guideway parameters have been investigated. Multicar 
vehicles have been found to provide better ride comfort than 
single car because of inter car constraints. The analytical 
techniques developed are suitable for analysis of vehicle/guideway 
interactions in maglev systems . 

Dynamic tyre force applied to road surface by heavy vehicles 
can cause premature road failure. Potter and Cole et al . [81] 
devised a measuring system of dynamic tyre force of commercial 
vehicles . Measured dynamic wheel load of instrumented vehicles has 
been used to investigate the effects of vehicle configurations, 
vehicle vibration and road strength on the prediction of road 
damage potential . 

Laib [82] investigated the dynamics of off-road vehicles on 
deformable soil. Here the original ground profile is altered 
significantly by the vehicle passing over the rough terrain. The 
vehicle transfer function has been used to compute effective 
ground profile which generates the actual vibratory motion in the 
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vehicle. The method has been illustrated for several soil type 
and for number of vehicles . 

Although technical literatures contain substantial amount c 
studies relating to various aspects of vehicle dynamics to trac 
irregularities, the following have not been adequately covered o 
not completely addressed - 

(1) Nonstationary excitation due to variable velocity run hav 
been studied mostly on lumped mass rigid vehicle models 
Analytical formulation for the flexible vehicle-track syste 
response with a general case of vehicle ground motio 
including spatially nonhomogeneous track roughness do no 
appear in the literatures . 

(2) In most flexible system, vehicle body has been modelled a, 
an elastic beam with constant cross sectional properties. I 
many vehicles flexible attachments are also present, om 
example being aircraft wings. These may undergo appreciabli 
dynamics induced by their support motion. The dynamics o 
flexible attachment of the vehicles has been ignored ii 
earlier studies. 

(3) The flexible structures are distributed parameter systems 
Their cross sectional properties and damping characteristic: 
may be nonuniform. In the solution of linear flexible vehicle 
model, mode superposition technique has been tried from time 
to time . Accurate knowledge of natural frequency and mode 
shape functions are essential requirements of the method. Ii 
many studies, FEM and other numerical techniques have beei 
used to obtain the vehicle modal parameters when cross 
sectional geometry is variable in space. Analytical techniques 
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for general variable section beams are not fully addressed in 
literatures 

(4) Flexible vehicle models with multiaxle - multiwheeled 
configurations including combined effect of rigid body 
heave-pitch-roll motion and coupled bending- torsional elastic 
modes have not been tackled for nonstationary track 
excitations . 

(5) Nonlinear suspension characteristics have been handled mostly 
by statistical linearisation and perturbation techniques for 
lower degrees of freedom model under stationary excitation and 
in more general case by Monte Carlo simulation. Statistical 
linearisation method may induce appreciable error when 
response probability function is different than the form 
assumed, which may happen in most practical cases of 
nonstationary input. Perturbation technique, on the other 
hand, has been applied with expansion of suspension forces in 
polynomial form. Nonlinear terms of the series are scaled by 
small parameters. Though theoretically valid, this brings out 
a restriction during application on its use to only weakly 
nonlinear systems. Monte Carlo simulation can take care of all 
the complexities involved in the problem. However, it becomes 
less attractive in terms of computer processing time and 
storage . 

1.3 PRESENT WORK 

In the present study an attempt has been made to fill some of 
the gaps mentioned above. Analytical methods have been developed 
to obtain response statistics of vehicle-track system induced by 
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track unevenness. The track unevenness is assumed to be a 
nonhomogene ous random process specified by a generailsed PSD 
function. A variable mean of the track profile superimposed over 
zero mean random process has been considered. Track is assumed to 
deflect vertically under the vehicle imposed loads. It is modelled 
as a beam of uniform mass and flexural rigidity supported by 
elastic and dissipative subgrade. 

Vehicle forward velocity has been assumed to be variable with 
time. Vehicle ground motion is modelled by a polynomial series 
with finite terms. Coefficients of the polynomial can be properly 
adjusted to depict any condition of the vehicle run. With the help 
of vehicle forward motion history, track input to vehicle can be 
obtained as a function of time . 

Vehicle has been modelled as a combination of discrete and 
continuous elements. Vehicle wheels, tyres and part of the 
suspension systems are lumped as unsprung masses . Vehicle body 
contains slender attachments. Flexible vehicle body and slender 
attachments are modelled as beam elements with variable mass, 
stiffness and damping distributions. The attachments are allowed 
to carry arbitrary number of concentrated masses at arbitrary 
locations . 

Sectional properties of the nonuniform beam has been modelled 
by power series . This representation is most general and any type 
of variation can be accommodated with desired accuracy. A method 
is outlined to find exact solution of nonuniform beam vibration. 
Natural frequencies and mode shape functions have been obtained in 
closed form for bending, torsion and bending- torsion coupled . 
dynamics of variable section beam. 
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As the natural frequencies and mode shapes are very important 
to understand the system behaviour, validation of the proposed 
theory has been sought by comparison with available published 
results and by conducting experiments . 

Governing differential equations for the vehicle- track 
systems have been developed. Continuous elements of the model 
involve partial differential equations with space variable 
coefficients. Equations of motion corresponding to lumped elements 
are ordinary differential equations. 

In the present analysis, the equations of motion of the 
continuous systems are first reduced to a set of linear ordinary 
differential equations in time dependent generalised coordinates 
using normal mode expansion techniques . The set of ordinary 
differential equations together with the equations of lumped 
masses are then recast into first order state space form. These 
equations are completely uncoupled using modal analysis. 

Analytical solution for response is obtained and expressions 
for second order statistics have been developed. The approach 
presented does not put any restriction on the nature of damping. 
The present method is applied to obtain the response statistics of 
transport aircraft with linear suspension parameters. Following 
models in order of increasing complexity have been tried: 

(1) Single Point Input Heave Model 

(2) Two Points Input Heave -Pitch Model 

(a) Rigid fuselage 

(b) Flexible fuselage in bending 

(3) Three Point Input Heave-Pitch-Roll Model 

(a) Flexible fuselage in bending with rigid roll 
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(b) Flexible fuselage in bending and torsion 

In the next part of the study, nonlinear behaviour of the 
suspension parameters is included in the analysis. Combination of 
air spring, orifice and coulomb damping has been treated with 
single point input heave model including tapered flexible wings. 
An instantaneous linear approach has been proposed. The method 
derives the benefit of linear theory while at the same time 
retains original nonlinear behaviour by incorporating 
instantaneous change in slope of the suspension forces at 
discrete operating points in time. 

System equations for the nonlinear model have been developed 
wherein the suspensions forces are nonlinear functions of 
instantaneous stroke length and velocity. The system equations 
are integrated numerically with digitally simulated track profile 
samples using Shinozuka's algorithm [83]. Suspension forces are 
determined at discrete time steps from known response of original 
nonlinear systems. Derivatives of the forces at discrete time 
steps are evaluated which are used to represent ' instantaneous 
spring and damping constant' . With these instantaneous suspension 
parameters, step by step application of the analogous linear 
model has been carried out to obtain second order response 
statistics. Reliability of the solution has been assessed by 
comparison with Monte Carlo simulated results . 



CHAPTER 2 


SYSTEM DESCRIPTION 

The system consists of a vehicle travelling over an uneven 
track supported by elastic and dissipative subgrade. Models have 
been developed for the system elements to study the coupled 
vehicle- track dynamics. These are being described below. 

2.1 VEHICLE FORWARD MOTION. 

Any general forward velocity pattern of vehicle may be 
described by the position of its center of gravity (c.g.) along 
the track with reference to a fixed datum. Let the position of the 
vehicle c.g. x c from the origin at time t be given by a 
polynomial of degree m as 

m k 

x (t) = E a, t (2.1) 

c k=0 K 

This form can easily be used, as mentioned above, to represent 

different conditions of vehicle forward motions with suitable 

selection of the coefficients of the polynomial. For example, a 

constant velocity travel requires a 2 = a 3 = ... = a m = 0 with a 0 

and a^ having non zero values. The vehicle with uniform 

acceleration or deceleration, can be represented with a 3 = a 4 = 

.... = a„ = 0 with a A , a n and a„ being non zero. In a situation 
m 0 1 2 J 

where acceleration or deceleration is variable, the coefficients 
a^ , a 4 , etc. are retained. 

A vehicle has, in general, front and rear axles. Location of 
the front and rear axle from the vehicle c.g can be obtained when 
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wheel base length is known. Thus one has 


X 1 (t) 

= x c ( t ) + 

111 k 

1 1 “ ^ a lk fc 

(2.2) 

x 2 (t) 

= x c <t) - 

m k 

^2 ~ ^ a 2k^ 
k=0 K 

(2.3) 

where a 1Q 

= a 0 + 1 1 

' a 20 =a 0 _1 2 

(2.4) 

and a lk = 

a 2k = a k' 

k=l, 2 , . . . , m 

(2.5) 


where and x^ are the instantaneous location of the front and 
rear axles measured from a reference point, 1 ± and 1 are the 
distances of the front and rear axles from the vehicle c.g. and 
wheel base length is the summation of 1 and 1 

JL Z 


2.2 TRACK ROUGHNESS 

The vertical height of the track along the longitudinal axis 
measured with respect to a flat datum at a distance x from the 
reference station can be represented by 


h(x) = h m (x) + h R (x) (2.6) 

where h (x) is a deterministic function describing the track mean 
and h_ (x) is a zero mean random process. The mathematical 
representation of h_ (x) using Stieltje's form is 


h R (x) 


= J dSj^Q) exp(jfix) 


(2.7) 


where j is the imaginary unit equal to v^T, ft is the spatial 

frequency and dS, (ft) is an orthogonal function with the 

n R 
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properties 


E [ dS, (Q) ] = p, (Q) dfi = 0 

R R 


( 2 . 8 ) 


and E 


ds h< n i> 


as^ 2 ) 


]■ 


T_ ( 

h R h R 


^1 ' ^2 ^ 


dfA dn. 


(2.9) 


where ^ (Q^,fl 2 ) is the generalised PSD function of the track 
R R 

random unevenness. Asterisk denotes complex conjugate. 

The Fourier transform of the track random roughness can be 
expressed as 


h R (fi) = - J h R (x) exp(-jQx) dx 


( 2 . 10 ) 


Using eq.(2.7) in eq.(2.10) 


V Q) = -£r 


WJ uu 

J { J dSh (Qf) exp ( j^iX) j- exp(-jQx) dx 


“00 -00 


Interchanging the order of integration and using the definition of 
Dirac delta function, it can be shown that 


h R (Q) 


I 



(n ± ) s (n-Q-j^ 


( 2 . 11 ) 


where 8 ( . ) is the Dirac delta function 

During variable velocity run the vehicle senses a nonstationary 
track input. The transformation of frequency from space to time 
domain would be governed by the relation 

0 ) = n V(t) ( 2 . 12 ) 
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where w is the temporal frequency and V(t) is the vehicle forward 
velocity at time t. This gives a nonlinear mapping of the temporal 
frequency cj into the wave number £2. 

Let the relation between the frequency and wave number for the 
track PSD be 


W 1 = ^ 1 V 1 an< ^ w 2 = ^ 2 V 2 


(2.13) 


where V 1 and V 2 are the vehicle forward velocities at time 
instants and t 2 respectively. Now from the definition of 

generalised PSD, one has 


$ h h^ n i' n 2 } 

ix K. 


= e | h R (Q 1 ) h* (n 2 ) | J 

00 

{I as h E (r l )5(r l- [J l>}{ 


= E 


dS h E <r 2> 5 ( r 2' n 2> 


00 00 


= J | E dS^ r 1 ) dS^{r n ) 


^2- 


] 5(r f n i 


) 5 (r 2 -£2 2 ) 


00 00 


-J J W r i' r 2 ,6(r l-V V l> 5(r 2- “2 1 v 2 )dr l dr ; 


- 00 - 00 


= “l ' V l' “2 ' V 2> 


(2.14) 


The function h m (x) describing the variable mean track is 
assumed as a power series 
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(2.15) 


£ 

h (x) = V h. 
m . z -u 1 

i=l 


x 


This representation is general and many shapes can be 
represented with this form accurately. Making use of eq.(2.1) in 
eq.(2.15), one has the mean track height experienced by the 
vehicle wheel at instant t 


h (t) 
m 



(2.16) 


Noting that, a polynomial of degree m raised to a power i, is a 
polynomial of degree m x i, eq.(2.16) can be expressed as [84], 


h m (t) =X h i £ c 
i=0 r=0 


i,r 


(2.17) 


where 


k = mi 


'i, 0 


C i/ r 



r 

E ( ki 

k=l 


r + k ) a, c. , for ri 1 (2.18) 

Kl 1 1 H KL 


The time derivative of h (t) 


can be written as 


£ k - 1 

h (t) = E h. E c. 

m i u 


i=l 


r=0 


i , r+1 


( r+1 ) t J 


(2.19) 


In case of vehicles with two axle and three wheel 
configurations, mean track levels at the various wheel locations 
would be 
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w .yvn 


( 2 . 20 ) 


k 


m2p‘ 


E (h.) E (c. ) t x 

i=0 p r=n 1 ' r 2 P 


r=0 


P=2 , 3 


( 2 . 21 ) 


where = (a 1Q ) 1 and 


(c. ) 


E (ki-r+k) a., (c. , K, ; r*l 


'i, r H a 1Q r ^ “Ik-i^-k^ll 


( 2 . 22 ) 


The expression for coefficients (c. ) _ are as ecr. 

x, r 2p 


(2.22) with 


lk 


(k=0 , 1 , . . . ) replaced by a pk where 


30 " 2 0 


a lk -a 2k -a 3k~ a k' k ~ 1 /2,...m 


The two numerals in the suffixes are used to denote the axle and 
the wheel locations. The first numeral is 1 for the front and 2 
for the rear axle. The second numeral becomes 1, 2 and 3 for 
centre, left and right wheels on a particular axel. 

The above track roughness model with the mean level and random 
parts, may be used with homogeneous as well as nonhomogene ous 
tracks. The analysis has been developed for the general 
nonhomogeneous unevenness, the homogeneous form being a particular 
case . 


2.3 VEHICLE-TRACK DYNAMIC MODEL 

Vehicle structural elements are elastic members . These are 
usually modelled as rigid masses to reduce the mathematical 
complexity of the analysis. However, idealising some slender 
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components as rigid masses may induce unacceptable errors . The 
present study attempts to use a combination of rigid and 
continuous idealisations of the vehicle. 

As mentioned in chapter 1, in the present study, the 
following linear suspension vehicle models with increasing order 
of details are considered : - 

(a) Heave Model 

(b) Heave-Pitch Model 
(i) Rigid fuselage 

(ii) Flexible fuselage in bending 

(c) Heave-Pitch-Roll Model 

(i) Flexible fuselage in bending with rigid roll 
(ii) Flexible fuselage in bending and torsion 
These linear models are being described below. 

2.3. 1 Single Point Input Heave model 

A model of the vehicle and the track is shown in Fig. 2.1. The 
wheel track contact is assumed to be at one point only. The 
aircraft fuselage is modelled as a rigid lumped sprung mass onto 
which are attached flexible left and right wings, designated by 
superscript / subscript L and R respectively. Each wing is 
considered to act as a cantilever beam having distributed mass, 
flexural rigidity and damping, that are variable along the span. 
Various stores and attachments carried by the wings are modelled 
as concentrated mass loadings. The landing gear assembly is 
idealised as a lumped unsprung mass with linear spring and 
viscous damper. Tyre stiffness and damping are also assumed 
linear. The runway track elastic behaviour is accounted for by 
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h R 




Fig. 2.1. Heave Model 



modelling the pavement as a Euler Bernouli beam with distributed 
flexural rigidity and mass taken constant along the span. The 
track is supported over elastic subgrade foundation having 
distributed stiffness and damping. These are assumed to be 
constant along the length of the track. The mass of the foundation 
is considered to incorporate the inertia effects under dynamic 
conditions . 

2.3.2 Two Point Input.* Heave-Pitch Model 

(a) Rigid Fuselage : - 

The aircraft fuselage is idealised as a rigid beam and wings as 
the flexible attachments. The vehicle and track models for the 
present analysis are illustrated in Fig. 2. 2 The track input 
through the wheel contacts are assumed to be at two points only. 
The vehicle body (fuselage) has pitching along with vertical 
bounce or heave motions. The two undercarriage shock absorbers and 
wheels are modelled as lumped masses with linear springs and 
dampers having only heave degrees of freedom. The wings are 
considered as variable taper cantilevered beams oscillating in 
coupled bending- torsion modes relative to fuselage. The wings 
may carry a number of concentrated masses placed at various 
locations along the span. The track is idealised as a continuous 
memember supported by elastic and dissipative subgrade. 

(b) Bending Flexible Fuselage 

The model of the vehicle and the track adopted is shown in 
Fig. 2. 3. The c.g of the aircraft is taken as the body fixed 


32 








origin. The fuselage is treated as a flexible continuous beam in 
bending with variable sectional properties. The absolute 
transverse displacement of the fuselage is composed of rigid 
body heave translation, pitch rotation of the longitudinal axis 
and elastic bending deformation relative to the rigid body 
displacement. The wings are modelled as cantilever beams of non 
uniform section having coupled bending - torsional motions. The 
front and rear undercarriage systems are modelled with linear 
shock strut springs and dampers. Tyre stiffness and dampings are 
treated as linear. The undercarriage masses are taken as lumped 
masses having only vertical heave motions. 

2.3.3 Three Point Input Heave-Pitch-Roll Model 
(a) Bending flexibile fuselage with rigid roll:- 

Aircraft model with tricycle landing gear arrangement is used 
(Fig. 2. 4). The front wheel follws the centre line profile while 
two main wheels follow the left and right profile of the track. 
The wheels of the vehicle in motion experience different track 
levels in the longitudinal and transverse direction. As a result, 
the vehicle is set to combined heave, pitch and roll motion. 

Vehicle body is considered, to be flexible in bending but 
rigid in torsional deformations. Flexible wings are considered as 
cantilever beams with roots attached to the vehicle body. The 
root motion of the beam induces dynamic loading on the span. Two 
main wheels are located under the wings. Wheels are treated as 
lumped masses, having only vertical degrees of freedom. 
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Inertial axis 



Fig. 2 . 3 . Heave-Pitch Model (Flexible fuselage) 









(b) Bending and torsionally flexibile fuselage :- 

The above three points input model is further improved by 
considering elastic torsional deformation of the vehicle body. 
Interactions of flexural and torsional modes of the vehicle body 
along with rigid body heave-pitch- roll motions have been 
considered. 


2.4 SYSTEM EQUATIONS 
2.4.1. Heave Model 

The system displacements are considered only. in the vertical 
plane and are designated positive upward (Fig. 2.1). The 
displacements z ^ and z 2 of the sprung and unsprung masses are 
measured from their respective static equilibrium positions. The 
transverse displacements of left and right wing w T and w„ are 
measured off their respective bending axes drawn from the roots . 
Deflection of the track mean line y at any instant t is taken 
transverse to its longitudinal axis . 

The equations of motion of two lumped masses are 


M s z l + c s ^ z i " z 2 ^ + K s ( z ]_ “ z 2 ) 


a 2 w, 


a w. 


3s l L E wL I wL (S L ) a ^2 


L S L = ° 


as R _ e wr i wr (s r ) 


SS R J S R=° 


(2.23) 
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M us z 2 + C us { z 2 - y < V b) - h ( x c > } 

+ K us |z 2 - y (x c ,t) - h (x c )} - C s ( z 1 - z 2 ) - K g (z 1 -z 2 )=0 (2.24) 

where M and M „ are the sprung and unsprung masses, K and C are 
s US s s 

the linearised spring and damping constant of the shock absorber. 

K and C are tyre stiffness and tyre damping coefficients. E I 
us us J w w 

are the wing's flexural rigidity, w is the transverse bending 
displacement of the cantilever wing. Subscript L or R is used for 
left and right wing. x c is the distance of the vehicle along the 
track and gives the wheel track contact location. 

The wings act like cantilever beams with base excitation as 
these are clamped at the fuselage and excited by the fuselage 
vertical displacement. The distributed aerodynamic lift as well as 
concentrated mass loadings at discrete locations are considered to 
act on the wing. The equations of motion of the left and right 
wing are. 


as! 


E wL I wL (S L } 


a 2 w, 


as: 


+ m wL (s L } 


a 2 w t 


at 


C wL (S L } 


aw T 


at 


= f w (s L' t} 


(2.25) 
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(2.26) 


o a 2 

s 2 r ^ w p i 

^2 L E ^R <S r’ -^2 ] + m wR (s R> 

R R 


3 "r , > 3W R 

7? + < W* ) — 


f S (S R' fc) 


where m w and c w are mass and viscous damping per unit span. The 
distributed applied force on the left wing beam is given by 


f (s, t)= 
w ' 


m wL <S L^ 2 1 + c wL (s L> Z 1 


- E m£{ w L (s,t) ”+ 2 ,)} 5 ( s T - sj + V^,(s,t) (2.27) 


L "k J 


where is kth concentrated mass loading on the wing at span 

location s, and V (s,t) is the aerodynamic lift force per unit 
length of the wing. These quantities can be particularised for the 
right wing by using proper subscript/superscript. 

There exists a similar expression for the impressed force on 
the right wing which can be obtained from eq. (2.27) by using 
subscript/superscript R. 

Assuming an elliptical distribution, lift on either wing can 


be expressed as [85] 


Vs.t) - L o v(t > 2 [ 1 - < s / V 2 ] 


1/2 


(2.28) 


where L Q = 


p C T S 
l _a L w 

2 7 t L . 


(2.29) 


in which p is the air density, C T is the lift coefficient, S is 

3. i-j f ** 

the wing surface area, and L w is the semi wing span. 
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The quantity V(t) is the aircraft forward velocity and wehave 


from eq. (2.1) , 


V 2 (t) 


- { x c !t> } " 


2 (m-1) 

D C 3 
r=0 


(2.30) 


2 

where ) 


and c' 


r a 


0 k=l 


S ( 3k - r) c^_ k 


a k < k +1 > a k + l 


r t 1 


(2.31) 


For vehicles where aerodynamic lift is not significant, the lift 
terms may be dropped from eq. (2.27) . 

The track is assumed to behave like a beam with continuous 
elastic support along its length. The wheel reaction acts as a 
moving dynamic excitation for the track beam. Its differential 
equation of motion for a single point input is given by 


E I 
P P 



m 


P 



d Y 

+ k f y = f (x, t) 

at y 


(2.32) 


where E I , m are the track flexural rigidity and distribued 
p p p 

mass; c f and k f are distributed viscous damping and spring 
constant for the foundation. These are assumed uniform along the 
track length. The applied force per unit length for one point 
input model is given by 
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f (X, t) 

hr 


Ks{ 


Y (x,t) - h (x) 


} 


+ K 


US 


{ 


- y (x,t) - h(x) 


S (x - x ) 
c 


(2.33) 


2.4.2 Heave -Pitch Model 
(a) Rigid Model 

The fuselage beam undergoes rigid body heave motion z and 
pitch rotation of the longitudinal axis © . Front and rear 
unsprung masses have only vertical degrees of freedom z- and z_ 

-L Zt 

measured from their static equilibrium positions . The wings 
attached to fuselage are assumed to have variable cross sections 
where the shear centre does not coincide with the centroid. The 
coupled bending -torsion oscillation of the wings is induced by 
fuselage heave-pitch motion. The transverse displacement and twist 
of the wings relative to the fuselage motion are denoted by w and 
0. All rotational displacements are assumed leading end up 
positive and linear displacements are assumed positive upward as 
indicated in Fig. 2.2 

The governing differential equations for the heave-pitch 
motions of the rigid fuselage are 


M s z + C sl (z * 1 1 9 Z - Z 1 > + K sl (z + HV z l> + C S2< Z - 1 2 e Z - Z 2> 


+ k s2* z " 1 2®z ' z 2* “ as. 


{ 


5 2 w_ 


E _ I T ( S T ) 
wL wL L 


as L J S L =0 


3Z E i 


E wR I wR (S R ) 


S 2 w t 




(2.34) 
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1 L S R=° 

wh.6ir6* C.«, K-,* C K 0 010 nose und Tnuin g 0 a.r‘ suspension dumping 

Sl Sl S A S A 

and spring constants respectively, I g is the mass moment of 
inertia of the fuselage, E w I w (s), G w J w (s) are flexural and 
torsional rigidities of the wing beam. Here 1 Q is the distance of 
the wing' s shear center at root measured parallel to longitudinal 
axis from the aircraft c.g. 

Equations of motion for the two unsprung masses are 


M ui z i + C si { z i ■ z ■ Vz } + K si{ Z i ' 2 
+ C ui{ ui {z.-htx^-y^. 


- 1.8 

1 z 


} 


t)j = 0 

i=l. 


2 


(2.36) 


i is equal to 1 for the nose wheel and 2 for the main wheel. 

The wing beam gets disturbance from the fuselage vertical 

motion and rotation at its support location. The equations of 

motion describing the coupled bending - torsion forced vibration 

of the wing beams are 
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J- 1 G J (s) -|® 
5S 1 w w ' 5s 


2 2 

t /„\ 3 8 / \ / 5 w 

W 8 ’ ~T5 + V s > e w (s) — 72 

ul ot 


& (s) + 
w 


e w 2 (s) c(s)l || + c w (s)e w (s) §| = e w (s) f w (s, t) -T w (s , t) 


(2.38) 


where I^ a (s) is the polar mass moment of inertia of the wing cross 

section about the shear center, m w (s) is the distributed mass 

along the span, e (s) is the off set of shear center from the 

elastic axis, c (s) is the distributed viscous damping to 

w 

translational motion and (s) is the distributed damping to 

w 

rotational motion. The impressed force f and torque T w per unit 
span are 


f w (s ' t) - - m ^ s) < 1 oK> - c w (s) < z "Vs* - m w (S) e w (s) e z 

+ e w (s) c w (s) ®z -Ji M k { z '- Vz + w <s ' t) } 5(s ’ s k> 

+ | M k e w ( S ) { 8 z +0 (S,t) } 5(s-s k ) + V w (S/t) (2.39) 

lr=1 v ' 
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r w (s - t) = - W 60 8 z + m w (s) e w (s) < 2 - 1 0 e z ) - 


# , ( s ) + 
w 


e w (s) C w (s) r e 2 + C w (s) e w (s) ( Z - 1 O 0 z ) - k ^ 1 I k 1 e z + e(s ' t} 


} 


5(s - s k ) + E M k e w (s) { 2 -1 0 e z + w(s ' fc) \ 5 ( s - s k } + M ac (s,t) 
k— X ^ 


(2.40) 


where p is the number of concentrated masses, M k and I k are the 

magnitudes of the kth concentrated mass and its moment of inertia 

about the shear center, V and M are the distributed lift force 

w ac 

and aerodynamic center moment per unit span of the wing. The left 
and right wing equations of motion are designated by placing 
appropriate subscripts L or R for the variables in eqs.(2.37) - 

(2.40) . 

The vertical lift distribution on the ( s , t ) is given by eq. 
(2.28) while aerodynamic center moment M (s,t) may be expressed 
as 


M ac (s,t) = M 0 c h (s)v2 (t) (2.41) 

P C 

, .. *a mac 

where M Q = ^ 

'C ' being the arodynamic centre moment coefficient and c, is 
mac xi 

the mean chord. 

The equation of motion of the track beam with foundation 
reaction taken into account is given by eq.(2.32). The track beam 
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is subject to two point loads and the applied loading per unit 
length can be written as 


fp(x,t)= - E { C u± ( z i - h(x) - y (x, t ) ) + K ui ( z i - h(x) 

y (x, t) } S (x - x i ) 


(2.42) 


(b) Flexible Vehicle body:- 

The vehicle body is treated as beam with variable cross 
section. The transverse displacement w^(r,t) at any station r is 
given by 
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f b (r,t) + — » b (r,t) 


(2.43) 


where is the flexural stiffness, m^ is the mass and c b is the 
viscous damping per unit length of the vehicle body. The impressed 
force f b and moment per unit length are 
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(2.44) 
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and 


M h (r,t) 



+ 


{ G wL J wL (S L ) 
{ G wR J wR (s R ) 



(2.45) 


The equations of motion of the two unsprung masses can be 
written as 


M ui*i + c ui { z i - - (' (x i' c) } + K ui { z i - h(x i> - 


C si { z r V 1 !- 10 } + K si { z i- Vh't’} " 


(2.46) 


(i=l, 2) 

The equations of motion of the non uniform flexible attachment in 
coupled bending - torsion are same as eqs . (2.37) and (2.38). The 
induced forces are however modified in view of the interaction of 
flexural modes of the fuselage. These are given as 


f w (s,t) « 


- m w (s) w b (l 0 ,t) - c w (s) w b (l Q ,t) 
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and 
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(2.48) 


The equation of motion for transverse oscillation of the 
track beam has the form same as in eq.(2.32) with the vertical 
force imposed on the track at the location of front and rear wheel 
is given by the expression (2.42) 


2.4.3 Heave-Pitch-Roll Model 

(a) Bending flexibile fuselage with rigid roll:- 

Transverse oscillation of the vehicle body can be described 
by eq.(2.43). However, the vertical force and bending moment per 
unit length induced in the vehicle body taking account of the 
rigid body rolling and interactions of the landing gear 
suspensions are given as 

f b <r , t ) = -[ C all { w b (r,t) -"ijjJ + K sll |w b (r,t) - Zll j jslr-y 
-J 2 [ C s2i{ »0 (r - t> Ai} + K s2i{ _Z 2i} ] S(r - 1 i> 
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(2.50) 


where C and K are the suspension damping and stiffness. As 

5 S 

mentioned earlier, two numerals in the suffixes are used to denote 
axle location and wheel location on the axle. The first numeral is 
1 for front and 2 for rear axle. Second numeral becomes 1 , 2 and 3 
for centre, left and right wheel respectively on the particular 
axle. The quantity w Q is the total deflection of the vehicle body 
at any transverse location from its longitudinal axis and is given 
by 


w Q (r, t) = w b (r,t) + g b (r) © r (t) 


(2.51) 


where g^r) is the distance of the point on the vehicle body at 
station r, measured transverse to the longitudinal axis. © r is the 
rigid body roll motion. The equation of motion governing the rigid 
body roll of the vehicle can be written as 
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where e 2 and e 3 are the distance of the wheels on the rear axle 
from the longitudinal axis of the vehicle body and are given by 


9b (1 2 )+s gl/ e 3-<9b (1 2 ,+s gR> 


The equations of motion for the front and the rear wheels are 


M ull z ll + C ull1 Z 11 - h ll (x l } ' y(x l' t} 
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(2.53) 


M u2i z 2i + C u2i{ z 2i - k 2i (x 2> ~ } 


+K u2i{ Z 2i- h 2i (x 2>- + C s2i{ Z 2i“ « 0 (1 i ' t > ““i (S gi' C > 


+ K s2ii Z 2i - W 0< 


li,t)-Wi(s gi ,t)} - 0 


(2.54) 


(i=2,3 ; w 2 =w L ; w 3 = w R ; S g2 - S gL ; S g3 - S gR ) 
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The coupled bending -torsional oscillation of the flexible 
attachment have the same form as eq. (2.37) and eq.(2.38) except 
that the applied vertical force and torque induced by fuselage 
motions and landing gear reaction are modified as 


f „ (s ' C) - - ra w (s) W t) - C w (B) W t) + V s,e w <s) -2 w 0 (1 0' t> 

ot 


+ V s,C w (s) ft v 'o (1 0- t) 'jj M k { w 0 (1 0’ t) + W(s ' t) } ‘ M k e w 


(s) 


» • -i ( * 

{ — 2 w^(l 0 ,t) + e(s,t) } J 5(s-s k ) - C s2i | w 0 (l 0 ,t) + 


at 


w (s, t) _z 2 i} +K s2i{ W 0 (1 0' t) + w(s,t)-z 2i }j 5(s-s gi ) + V w (s,t) 


(2.55) 


and 


rjs.t) = M ac (s,t)- i wa (s) 


at 


2 W 0 (1 0' t)+ m w (s)e w (s) W 0 (1 0' t) 


< V s) + e w (s)C w (s) >ft W 0 (1 0' t} +c w (s)e w (s)w 0 (1 0' t) 


P - a 2 


E I ] C { — 2 w o (1 0' t) + 0(s ' t) } - 

k =i K at u u 


p f ♦ • • • 

E M k 1 e w (s) w o (1 o ,t:) + 

k — 1 . 


5 (s-s k ) 


(2.56) 


The forces on the left and right wing can simply be written 
by placing subscript L and i=2 for left attachment; R and i=3 for 
right attachment in eqs.(2.55) and (2.56). 

It is assumed that the track beam ' deforms in the vertical 
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direction and amount of torsion on the track because of 
eccentricity of the wheels is small. The equation of motion for 
the transverse oscillation of beam is same as eq. (2.32) with 
impressed force given by 


“ C ull^ Z 11 ~ h n y(x,t) } + K u ii{ z n h n ^ 


y (x, t) 


>] 


5 (x-x 1 ) 


3 

- l [ c u2±< Ki ~ h 2 i‘ x) + 


Vi* z 2i - h 2i (x) -y (x ' t) > ] 5 < x - x i> 


(2.57) 


(b) Bending and torsionally flexibile fuselage 

In this model, the vehicle body is assumed to have elastic 
bending and torsional deformation induced by track roughness 
input. Total deflection of the vehicle body w Q (r,t) at any 
transverse location from its longitudinal axis is given by 


w Q (r , t) = w b (r,t) + 9 b (r) © b (r,t) (2.58) 

where © b (r,t) is the torsional deformation of the vehicle body 
measured positive counterclockwise from untwisted configuration. 

The transverse displacement in bending can be represented by 
eq.(2.43) where the expression for the impressed force f b has to 
be modified for the interactions of elastic torsional modes and 
thus w q from eq.(2.58) needs to be substituted in the expression 
f °r f b . 

The torsional deformation of the vehicle body is governed by 

the differential equation CENTRAL LlBHARI 

\ 1. T. f KANPUR 

125696 


51 



dr 


{ 


G b J b ^ 


d& t 

~dX~ 


}- 


W r) 


at 


7 - V r >' 


se l 

at 


-r b (r, t) 


(2.59) 


where G, I,, I, and are torsional stiffness, mass moment of 
b b br d 

inertia and rotational viscous damping per unit length of the 
beam. The applied torque r b per unit length is given by, 


r b (r,t) 




(r, t) 



w Q (r, t) 



e i (r) 


5 (r- 


V 


[{ 


E _ I x (s T ) 
wL wL L 


a 2 w x 


as: 




a 2 w 


E ^,1 _ (s„) 
wR wR R 


R 


as 


S R=° J 


S(r-1 2 ) (2.60) 

The equations of motion for the front and rear wheels can be 
obtained from eq.(2.53) and (2.54) after substituting w Q from 
eq. (2.58) . 

The coupled bending -torsional oscillation of the flexible 
attachment are induced by transverse displacement and rolling of 
the fuselage including the elastic modes. The governing equations 
of motion can be represented by eq.(2.37) and (2.38). Induced 
dynamic loading will be modified on account of elastic bending and 
twist and thus w Q from eq.(2.58) has to be used in forcing terms. 

Transverse oscillation of the track beam is given by 
eq.(2.32.) with forcing term as in eq.(2.57) 
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CHAPTER 3 


SOLUTION TECHNIQUE FOR LINEAR MODELS 

The system equations consist of a set of ordinary and partial 
differential equations . These can be solved with known initial and 
boundary conditions and input statistics . The linear models 
employed in describing the system yield coupled equations of 
motion. The present approach first reduces the partial 
differential equations of motion of continuous elements into a set 
of ordinary differential equations employing modal expansion 
technique. Subsequently, the complete set of system equations 
(including lumped mass equations of motion) are decoupled using 
complex modal analysis which are further processed to obtain 
response statistics . 

3.1 MODAL ANALYSIS OF FLEXIBLE STRUCTURES 

Flexible structures are characterised by parameters, such as 
mass, stiffness and damping that are functions of spatial 
variables. Continuous structures represent distributed parameter 
systems. The motion of a given structure is governed by partial 
differential equations to be satisfied over the domain of the 
structure and boundary conditions to be satisfied at the 
boundaries of the domain. The response of a -flexible structure may 
be obtained by modal analysis which requires the solution of the 
eigenvalue problem. For distributed parameter systems, there 
exists an infinite set of eigenvalues and corresponding eigen 
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functions. The displacement may be assumed in the form of a series 
of space -dependent eigenfunctions multiplied by time dependent 
generalised coordinates . This helps in transforming the partial 
differential equations into a set of independent ordinary 
differential equations. 

For most of the cases, the physical parameters of the flexible 
structures are not uniform in the domain of the structures . For an 
arbitrary variation of cross sectional properties, the closed form 
solution of the partial differential equations are not always 
feasible. The alternative is to use approximate solutions. This, 
however, may produce appreciable error in predicting natural 
frequency and displacement shape function specially in higher mode 
excitation. In the present study, a method has been outlined to 
obtain exact solution for the natural frequencies and mode shapes 
for a beam with arbitrary variation of cross section. This has 
been utilized for the modal analysis of the continuous elements in 
different vehicle models. The different steps involved in modal 
analysis are, solution of free vibration problem, development of 
orthogonality condition and discretization of the equations of 
motion of flexible structures in time dependent generalised 
coordinates . 

3.2 NATURAL FREQUENCIES AND MODE SHAPES 

3.2. 1 Bending Vibration 

In absence of external disturbance, the free vibration 
equation of non uniform beam can be written as 
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,2 r 


as 


ei (s) 


a 2 w 

as 2 J 


+ m(s) 


a 2 w 


+ C (s) 


aw 


at' 


at 


(3.1) 


Let the beam deflection be given by 

w ( s , t ) = W(s) exp (C+ jw)t (3.2) 

where W(s) is the mode shape function, o is the damped natural 
frequency and ^ is a factor related to the damping. 

Substitution of eq.(3.2) in the eq.(3.1) yields 



EI (s) 



+ m(s) (£+jo) 2 W (s) + c(s) (C+jw) W(s) - 0 

(3.3) 


Let. the variation of moment of inertia, cross sectional area 
and distributed damping be expressed as general polynomials whose 
coefficients can be matched to represent any desired variations. 
It is assumed that distributed damping is proportional to the 
distribution of crossectional area and hence to the mass. Thus we 

have 


I ( S ) = KO) f 1 (s) : 

A ( S) = A ( 0 ) f 2 (S) : 

C ( S ) = C (0) f 2 (s) 

where 1(0), A(0) , c(0) are the moment of inertia, cross sectional 

area and damping at the reference section. 

Substitution of eq.(3.4) in eg. (3.3), yields 


f, (s) = I Pi s' 
1 i=0 


f 2 (s) = £ q i 5 


i=0 


(3.4) 
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- , . d 4 W - 

f i (s) + 2 


df. 


d 3 W 


d 2 f x d 2 W 


ds ds' 


ds 2 ds 2 


A f 2 ( S ) W(S) 


(3.5) 


where X 4 = — [ (“ 2 ~ C 2 - DC) - jw( 2C + D) J (3.6) 

with D = c(0) / pA(O) 

Introducing a new independent variable z=As in eg. (3.5) 
transforms it to 


d 4 W d 3 W ^ ^1 d 2 W 

f i (z) -7JT + 2 — ~7~3~ + T2 T2 - f 2 (z) w = 0 


dz 


dz dz' 


dz dz" 


00 _ , 00 
where f 1 (z) = £ P k z ; f 2 (z) = I 2 
k=0 <-0 - 


„ pj- V- 

with p k = P k / X and q fc = / X 


(3.7) 

(3.8) 

(3.9) 


Let a series solution of the equation (3.7) be of the form 


W(z) = £ c. z 1 
i=0 


(3.10) 


Using W and its derivatives from eq.(3.10) along with 
eq.(3.8) in eq.(3.7) and following the product rule of two series 
[84] , one has 

00 

I ( P n + Q n + R n ~ S n } ^ = 0 (3.11) 

n=0 

where 
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n 


P n = I (n-k+l) (n-k+2) (n-k+3) (n-k+4) p k c n _ k+4 (3.12) 

k=0 
n 

Q n = £ 2 (k+!) (n-k+l) (n-k+2) (n-k+3) P k+1 c n _ k+3 (3.13) 

k=0 

n 

R n = Y ( k+1 ) ( k+ 2) (n-k+l) (n-k+2) P k+2 c n _ k+2 (3.14) 

k=0 


n 

S n = I ^k c n-k (3.15) 

k=0 

Since eq.(3.11) must be satisfied for every value of z, it implies 
that coefficient of each power of z must vanish. Thus we have 


c n+4 


1 

(n+1) (n+2) (n+3) (n+4) p 


- Y (n-k+l) (n-k+2) (n-k+3) 

k=l 


I ^k E n-k 


k=0 


(n-k+4) 


p k C n-k+4 


n 

- £ 2 (k+1) (n-k+l) (n-k+2) (n-k+3) p k c n _ k+3 

k=0 


n 

- £ (k+l) (k+2) (n-k+l) (n-k+2) p k+2 c n _ k+2 J (3.16) 

k=0 ( n = 0,1,2, . . .) 

Equation (3.16) above gives a recurrence relationship which 
can be used to determine the unknown coefficients of the series W 
except the first four coefficients c Q , c^, c 2 and . These have 
to be determined with the help of the beam boundary conditions . 
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After expanding eq.(3.16), it can be shown that all other 
coefficients cs ( i = 4,5,...) are linearly dependent on the first 
four coefficients c Q , c^, c 2 , . Thus one can write 

3 

. £ i = E L ik 5 k (3 - 17 > 

k=0 

where L j_] c = f ^ ( P r , q r . A) , i=4,5,...; r=l,2,... 

Expressing all coefficients of the series in terms of first four 
coefficients, it is now possible to have the solution in the form 


00 3 

W(S) = c Q + SjXs + c 2 x 2 s 2 + 5 3 x 3 s 3 ♦ E E L ik 5 k 


A^S 1 


(3.18) 


The above equation may be re expressed, 

3 

W(S) = E g. ( p v ,q v ,A,s) (3.19) 

i=0 1 K k 

The application of appropriate boundary conditions of the 
beam in eq.(3.19) results in the homogeneous equation 

B c = 0 (3.20) 

where the elements b^ of the matrix B contain W and its 
derivatives evaluated at beam boundaries. The element 5/ of vector 
c contains unknown coefficients c Q , c 1 , c 2 , c 3 . For non trivial 

solution of eq. (3.2 0) , the determinant of the matrix B needs to 
be equated to zero. The roots of the characteristic determinant 
gives multiple values of A which upon substitution in eq.(3.20), 
yields the vector c, which can be utilized for evaluation of the 
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jmode shape function W. The damped natural frequencies are obtained 
from eq . (3.6) 

3.2.2 Coupled Bending-Torsion Vibration 

The method developed for pure bending vibration is being now 
extended for a variable section beam in coupled bending-torsion 
piode . 

Assuming free undamped harmonic vibration in the normal modes 
<§ne has 


w((s,t) = W(s) exp (jut) ^ 

$.(;s,t) = <p( s) exp (jut) 

tyhere W(s) and <p ( s) are mode shape functions for the bending and" 
tiqrsi'onal' displacements and u is the undamped natural frequency of 
t|he beam in its natural mode . 

Substitution of eqs.('3.21) in coupled" differential equations 
qfr motion for the beam in free vibration, yields 


d 2 W 


4s 




EI.(S) 


ds' 


} - “ 2 { 


m(s)W(s) - e(s)m(s)^(s) 


= 0 


A AA 


ds V 


< GJ(s) 


"3s 


V 2 

- V + Ct) 


l a (s)0(s) - e (s)m(s) W(s) jj = 


C'3 . 22 ) 

('3\ 23) 


As in eq. (3.4) , variation of the sectional properties may b"e 
expressed ' as polynomial 1 functions whose coefficients can be 
selected to match a given variation. Let 
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00 


I (s) = I (0) f. (s) = 1(0) E p.s' 
1 i=0 


m (s) = pA(0) f, (s) = pA(0) E qj s' 
A i=0 


J(s) = J(0) f,(s) = j(0) E r. s' 
J i=0 


I a ( S ) . I a (0) £ 4 (s) = I„(0).I 0 U. s- 


e(s) = e(0) f(.(s)= e(0) E 9^ s ' 
5 i=0 


Substitution of equations (3.24) in eqs.(3.22) and (3.23) 


d ‘ ff, (s) -^5 


ds 




2 11 


ds' 


| - A 4 |f 2 (s) W ( s) - e(0)f 2 (s)f 5 (s)^(s) H ( 


3 § - | f 3 ( s) -df~} + { k a0 f 4 (s) 4>{S) " e(0)f 2 (s) w < s > f f 


,4 co 2 pA ( 0 ) „ El ( 0 ) u 2 I a (0) 

where A El (0) ; P “ GJ(0) ; k a0 pA(0) 


Again transforming from s to z coordinate, where 
eqs.(3.25) and (3.26) transform to 


f 1 (z) 


,4 t - t df., ,3 T - T d 2 f- ,2 T - 7 

d W , ,, ldW , 1 dW i nr r c / r,\ 

+ 2 -35- j - j - + 7-2 - 771 ' 1 w(z) f 2 (z) 


{ 


dz" dz" dz‘ 

e ( 0) f 2 ( z ) f 5 ( z ) 4>(z) | = 


(3.24) 

yields 

(3.25) 

(s)j = 0 

(3.26) 

(3.27) 

z = As, 

(3.28) 
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f 3 (z) 

d 2 d 

dz 2 

where 


W(z) = 

: W (AS) 

and 

00 

f 1 (z) 

= .E Pi 
i=0 


00 

f 3 ( Z > 

- 1 f i 
i=0 


00 

f 5 (z) 

= .E 9i ; 

3 r\ 


df 


3 dxp 


{ 


cE a£- + A M k : 0 f 4 (z) * {z) - 


e (0) f (z) f (z) W(z) 


} ■ 0 


(3.29) 


i=0 


l l 

i 2 ' ^1= q^A 


i=0 


(3.30) 


Assume the shape functions in the polynomial form 

00 00 , 

W(z) = £ b.z 1 ; l(z) = £ c. z 1 

i=0 i=0 


(3.31) 


Substituting W, <p and their derivatives along with egs . (3.30) in 

eqs. (3.28) and (3.29) and following the series product rules 
[84] , one has 


oo , n 


n 


£ 1 £ P k (n-k+l) (n-k+2) (n-k+3) (n-k+4) fc> n _ k+4 + 2 £ P k+1 (k+1) 


n=0 V k=0 


k=0 


n 


(n-k+1) (n-k+2) (n-k+3) b n _ k+ 3 + £ P k+2 (k+2) (n-k+1) 

k = 0 


n 


n k 


(n- 


-k+2) Vk+2 \l Q ^ b n-k + e(0) k 5 0 9 k-l C n-k 


n n 
z =0 


(3.32) 
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00 r n 

1 { Ei k 

Lo 1 k =° 


E r. (n-k+1) (n-k+2) c k2 + £ (k+1) r k+l (n k+1) c n-k+l 

In K k=0 


+ 6 - k " x2 


X 2 IB e(0)E E q-, g k _! 6 n _ k } z n = 0 (3.33) 

k=0 1=0 X K X J 


Since eqs . (3.32) and (3.33) must be satisfied for every value 
of z , it implies that coefficients of each power of z must vanish. 
Thus one has 


(n+4) ! p n L k=0 


n k 


E q k 6 k - e ( o) E E q k g k _i c 

- K n K k=0 1=0 


E p, (n-k+l) (n-k+2) (n-k+3) (n-k+4) ^ n _ k+4 

k=l 


2 £ Pk+l (k+1) (n ~ k+1) ( n ~ k+2) (n-k+3) b n _ k+3 

k=0 


- E P k+2 (k+1) (k+2) (n-k+1) (n-k+2) 6 n _ k+2 

k=0 


(3.34) 


— — I \ 2 t3 e (o) e E q k i k _i £ n _ k 

(n+2 ) ! r Q [ k=0 1=0 * K 1 n K 


- k «0 J 0 S k 5 n-k -J/k (n - k+1) (n - k+2) Sl-k+2 


-l < k+1) J k + 1 (n - k+1) 5 n-k + l 
k=0 


(3.35) 


( n = 0,1,2, 
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Equations (3.34) and (3.35) give recurrence relationships for the 
coefficients b^ and c^ (i=0, 1, 2, . . . ) . After expanding these two, 
it can be shown that coefficients b^ (i=4,5, ...) can be expressed 
in terms of the four basic coefficients b Q , b^ b 2 and b 3 of the 
series W and two basic coefficients c Qr c 1 of the series b- 
Similarly, c^ (i=2,3,...) can be expressed in terms of b Q , b^, b 2 , 
b 3 , 5 0 and c 1 . Thus one can write 


bi = E L ik 5 k + I M. 5 

1 k=Q IK K k=Q IK K 


3 , 1 , 

c i = E L ik b k + l M ik c. 
1 k=0 1K K k=0 1K 


; i=4 , 5 , 


; i=2 , 3 , 


(3 .36) 


The terms L^ k , Mj_ k and L^ k , M^ k are to be obtained from 
eqs.(3.34) and (3.35) in terms of known parameters p k , q k , r k , u k , 
g k , ^aO' e w (°) anc ^ the unknown frequency parameter A. From 

expressions (3.34) and (3.35), it can be seen when e(s) = 0, the 
coefficients b^ and c^ are independent of each other. This 
signifies uncoupling of the equations of motions. 

Solution for the free vibration in coupled bending torsion 
can now be expressed as 


W(s) = b Q + 


b 1 (As) 


b 0 (As) z + b.. (As) 


00 

+ E 

i=4 


3 

E : 

k=0 


J ik 


1 

E 

k=0 


M 


ik 



(3.37) 
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00 ( 3 / 1 / _ ^ 

<Ms) = c + c Us) + E E L ik b k + E m., c k (as) 

U 1 i=2 k=0 1JC K k=0 1JC 


(3.38) 


The six independent constants of integration can be found using 
six boundary conditions - four on W and two on <p . Use of the 
boundary conditions yields homogeneous equations of the form 

B «» = 0 (3.39) 
where the elements of the characteristic matrix B contain W, <p and 
their derivatives evaluated at beam boundaries . The vector « 
contains the unknown constants of integration b^ (i=0,l,2,3) and 
c ± (i=0 , 1) . 

As mentioned in case of bending vibration, the nontrivial 
solutions of the eq.(3.39) yields eigenvalues and eigenvectors 
that can be utilised to obtain natural frequencies and mode shapes 
from eqs.(3.27), (3.37) and (3.38) 

i 

3.3 ORTHOGONALITY OF THE MODES 

Given a set of eigenfunctions, the orthogonality of the normal 
modes can be proved. Let us consider the more general case as in 
the coupled bending - torsion vibration. 

Assuming the solutions of the eigen value problem in two 
distinct modes i and k as VL(s); <p^( s) (corresponding to natural 
frequency u>^) and W k (s); $ k ( s) (corresponding to natural frequency 
o k ) one has from eqs.(3.22) and (3.23) for the ith mode 

H 2 r d ? w 1 0 r 1 

— 2 EI(s) 2 = u i tn(s)W. (s) -e(s) m(s)<p. (s) (3.40) 

ds L ds J L 1 J 
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(3.41) 


d 

ds 


GJ ( S ) 


dtfK 

ds 


2 

1 


e (s)m(s) W. ( s) - 1 (s) <p ■ (s) 

1 OL X 




Multiplying eq.(3.40) by W k and eq.(3.41) by b k and integrating 
over the beam length, for any combination of clamped, hinged and 
free ends yields 


El (s) 


d 2 W i d 2 W K 


ds 2 ds 2 


XJ 

ds=co 2 J | m(s)W i (s)W k (s) 


- e (s) m(s) (s) W k (s) j- ds (3.42) 


and 


I 


G J ( s ) 


d<p ± d <*> k 
ds ds 


J_i 

ds = w 2 J | I a (s) 


- m(s) e(s)W i (s)^ k (s) J- ds (3.43) 


Similarly for the kth mode of vibration, replacing i by k in 
eqs . (3.40) and (3.41), multiplying eq.(3.40) by W ± and eq.(3.41) 
by cp ^ and performing integrations over the length of the beam, one 

has 


L 

' EI(s) 

0 


d 2 W ± d 2 W K 

T "~2 , 2 

ds ds 


ds 


L 

o 2 | { m(s)W k (s)W i (s) 
0 


- e (s)m(s) <£ k (s)W^ (s 



(3.44) 
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and 


d (/> . d$. ( 

J -as ds = “k j { V E) V s) V s> 

0 0 

- m(s) e (s)W k (s) <p^ (s) j- ds (3.45) 

Subtracting eq.(3.44) from eq.(3.42) and eq.(3.45) from 
eq. (3.43) , yield 

L L 

(w?- co^) J m(s) W k (s) W i (s) ds- to? J e (s)m(s) ( s) (s) ds 

0 0 

L 

+ to^ J e(s)m(s)<^ k (s)W i (s) ds = 0 (3.46) 

0 

and 

L 

(co?- co k ) J I a (s) <?!> k (s) ^(s) ds- to? J e(s)m(s)W i (s)^ k (s) ds 

0 0 

L 

+ co k J e(s)m(s)W k (s)«^ i (s) ds = 0 (3.47) 

0 


Combining eqs.(3.46) and (3.47), yields 


U 

(£J i" w k ] I { «i(s)W k <s)W.(s) + I a (s)<P i (s)<P k (s) 


e(s)m(s) ^ i (s)W ]c (s) +<£ k (s) W i (s) j |ds=0 


(3.48) 
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For it follows 

L 

| m (s) W k (s) W i (s) + I a (s)<^ i (s)^ k (s) - e(s)m(s) U i (s)W k (s) 

+ ^ k (s)W i (s)J | ds =0 (3.49) 

Equation (3.49) provides a combined orthogonality condition 
of normal modes in coupled bending - torsion vibration. 
Orthogonality conditions for independent bending and torsional 
modes can be obtained from eq.(3.49) as a special case of coupled 
modes for which e(s)=0. 

Thus for only bending of beam with <f> { s)=0, the orthogonality 
conditions may be stated as 

L 

S m(s) W. (s) W, (s) ds =0 if i*k (3.5 0) 

0 IK 

and for only torsion with W(s)=0, the orthogonality condition 
becomes 

L 

1 I (s) <f>.( s) <p, (s) ds =0 for i*k (3.51) 

q Ol X K 

The relations given by eqs.(3.50) and (3.51) agree with the 
results derived for independent bending and torsional vibrations 
in ref . [86] . 
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3.4 DISCRETIZATION PROCEDURE OF EQUATION OF MOTION OF CONTINUOUS 
ELEMENT 

3.4.1 Bending Vibration 

Assuming transverse displacement of the beam as 


w (s, t) =£ W i (s) 7 ? i (t) (3.52) 

i=l 


where W^s) and r? i (t) are the mode shape function and generalised 
coordinate associated with the ith normal mode. Substituting this 
in eq. (2.26) for general forced vibration equation of beam in 
bending and multiplying both side of the equation by W k one has 


d 2 W 


“ r A* ( Q ) 

£ -==- i EI(s) 2 ~ f W k (s) + m(s) w i^ s ) W k (s) 

i=l>- ds^ ( ds > 


+ C ( S ) W i (s)W k (s) 


(t) 


f (s, t)W k (s) 


(3 .53) 


Assuming damping proportional to the mass distribution, one 
has from eq. (3.3) 

EI(s) } ={ (w?-Ci-DCi) -jw i (2C i +D) |m(s)W i (s) (3.54) 

Introducing eq. (3.54) in eq. (3.53), integrating in the beam domain 
and invoking the orthogonality property of the normal modes 
from eq.(3.50), one has 
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(3.55) 


7] i (t)+ DTi.(t) +| (w? -Ci-DC i )-jw i (2C i +D) | 7J i (t) = Q i (t) 

1 = 1 , 2 , . . . 

where Q^(t) is the generalised force given by 

L 

Q i (t) = (l/Mg^) | f(s,t) w |s) ds 

0 


!3 .56) 


is the generalised mass associated with the ith normal mode 
and is given by 



L 

m ( s ) W? ( s ) ds 


(3.57) 


Here the damping term is decoupled as damping is assumed to be 
proportional to the mass. Considering the case of general damping 
and using undamped modes of free vibration, the discretized 
equations have coupled damping terms. The discretized equations 
for this case becomes 


7). (t) + £ D.. 7? (t) + OK 7). (t) = Q. (t) 

1 k=l K ii i 


!3 .58) 


i=l,2, 


L 

where = (1/M Jc(s) W^(s) W^Cs) ds (3.59) 

0 

3.4.2 Coupled-Bending Torsion Vibrations 

Let the transverse bending deformation and twist of the beam be 
expressed as 
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W ( S , t ) = I w, (s) T? i (t) 

i=l 


0 (s, t) = Z <Pj_ (s) 7]^ (t) 
i=l 


(3.60) 


Substituting eqs.(3.60) in eq.(2.37) for the general forced 
vibration equation of coupled bending -torsion of the beam , 
multiplying both sides by and integrating with respect to s 

from 0 to L 


l J[^{ 

i=l J L ds^ l 


El (s) 


d 2 W. 

x 


ds^ 


W k (s)T7i + m(s) w i (s) W k (s) 77. (t) 


m(s) e (s) (s) 


W k ( s) 


7) • (t) + C (S) 


W i (s) 


w k (s) 




- e (s) c (s) (s) w fc (s) T) i (t) j ds = J £ (s, t) W k (s) ds (3.61) 

0 

Again substituting eqs.(3.60) in eq.(2.38), multiplying both 
side by $ k and integrating in beam domain 


l [ J 

i=l L { 


d 

1 


GJ (s) 


dcp 


"i vw _ x a VD ' y i'°^k- 'i 


<t>,( S) 7). (t) - I (S) <(>. (s)<£, (S)7). (t) + 


? :(t) -{< 


m(s) e (s) W. (s) <£,(s) 77. (t) + e^(s) c(s)}- <p. (s) 0 v (s) 77 ^ (t) 


+ e(s)c(s)W i (s)^ k (s)7? i (t)l ds 
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| f(s,t) e (s) 0 k (s) - r(s,t) <P k (s) | ds 


(3.62) 


Using egs (3.22) and (3.23) in eqs.(3.61) and (3.62), one has 


l J m(s)W i (s) -e(s)m(s)0 i (s) |w k (s) 7? i (t)+ m(s) W ± (s) W k (s) 7] ± 

1=1 0 J 

-m(s) e(s)^ i (s)W k (s) 7?^+ c (s) W i (s) W k (s) i^-e (s) c (s) <p ± (s) (s) tj J ds 


= | f (s, t) W k (s) ds 


(3.63) 


i i B( I a (s) <^ i (s) -e(s)m(s) w i (s) j-^ k (s) T? i (t) 


+ I a (s) <f>^ (s) <p k (s) Tj i (t) - m ( s ) e ( s ) ( s ) 0 k ( s ) i? i (t) +4> ± {s) (p^is) 


jtf- (s) +c (s) e 2 (s) j- 7) i (t) -e (s) c (s) <^> k (s) W i (s) 7? i (t) j ds 

L 

= J | r(s,t)0 k (s)-f (s,t)e(s)0 k (s) | ds (3.64) 

0 

Combining eqs.(3.63) and (3.64) and invoking the orthogonality 
property of the coupled modes eq.(3.49), yields after 
simplification 


7ji (t) + l D ik \(t) + 7 ?i (t) = Q ± (t) 

JC — X 


„ (3.65) 
i=l / 2 r ... 

where the modal damping coefficients -D, generalised force Q 


71 



associated in ith mode are given by 


D.,_ = 


ik M. 


gi 


Jj 

J c(s) VL(s) W k (s) +| #(s) + e 2 (s)c(s)j <^(s) ^ k (s) 


e(s)c(s)j <f>^(s )W k (s) + W.(s)* k (s)} J ds 


(3.66) 


L 

Q i (t) ■ 4— J [ f <= 

where generalised mass M . is given as 
L 

M gi =J jm(s) W? (s) + I a (s)0?(s) - 2m(s) e(s) <p^ (s) (s)j-ds (3.68) 
0 


fi 


,t) W. (s) -e (s )<p. (s) 


} 


+r(s,t)^(s) ds 


(3.67) 


3.5 DISCRETIZED EQUATIONS FOR CONTINUOUS ELEMENTS OF VEHICLE 
MODELS 

Following the procedure discussed in the previous section, the 
discretized equations for the continuous elements in each model 
can be obtained. 

3.5. 1 Heave Model 

In this model, the continuous elements tapered wings and track 
are represented by distributed parameter beam element. 

Let represent the ith normal coordinate for the wing. The 

discretized equations for the left and right wing are (from 
eq.3 .55) 
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\,i (t) +Dl, wi (t) + { <“Li- <Li- D «Li> - 3“ Li (2C Li + D)} 


7 >wi (t) * 


and 


^ (t) +D VL (t > + 


=wi (t > 


(3.69) 


( “Ri- «Ri- D «Ri> 


' ju Ri (2 ?Ri + D) } m 




(3.70) 


where w. is the damped natural frequency of the wing and Q . is 
the generalised force associated with the ith mode. The 
generalised force and generalised mass can be obtained for the 
wing beam from eqs.(3.56) and (3.57) with proper substitution of 
left and right wing quantities. 

The track equation of motion (2.32) can be discretised in 
terms of the normal coordinates i) . as 


7] pi (t) + B T? pi (t) + { ( Wpi - BC ± ) - 3 w p i< 2 £i +B ) }v (t) 

= Q .(t) (3.71) 

where w p is the track damped natural frequency; £ is a factor 
related to the foundation damping and the constant B is c f/ m p - The 
displacement function is given by [8 6] 

(x) =Y^sin y^x + Y 2 cos y^x + Y^sinh y^x + Y^cosh y^x (3.72 

in which Y ±r y 2 , Y 3 and Y 4 are constants of integration and the 
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parameter r is expressed by 


4 

r = 


m 


El 

P P 


{ 


( 


2 J2 

w. - e ■ 


B?) 


j Wp (2e+B) 


} 


(3.73) 


where w^=k^/m p 

The eigen parameter r and constants of integration , Y 2 , y 3 
and Y 4 can be found with the appropriate end conditions of the 
track beam. The generalised parameters for the track beam can be 
obtained from eqs.(3.56) and (3.57) by using track- subgrade 
quantities . 


3.5.2 Heave -Pitch Model 
(a) Rigid Vehicle Body:- 

In this model, the partial differential equations of motion 
representing the coupled bending torsional motion are discretized 
to their time dependent generalised coordinates. The equations for 
left and right wing, thus become (from eq.3.65) 


7 >Li (t) 


k=l 


D L V 

lk 


77 Lk (t) 


+ G) 


Li 




(3.74) 


and 


7J Ri (t) 


+ E D® 

k=l 


ik ^Rk^ 


w Ri ^ {t) 


2* 

wi 


(t) 


(3.75) 


where the modal damping coefficients D, generalised force Q and 
generalised M for each wing can be evaluated using eqs.(3.66) - 
(3.68) . 

The track discretized equations would be same as in heave model 
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except generalised forces are to be evaluated with the expression 
for impressed force given for a two point input model in eq. (2.42) 


(b) Flexible Fuselage in bending: - 

The equation of motion for the ith normal coordinate of the 
flexible fuselage in bending can be written as 


”bi (t) 


7 , bk (t) 


H bl <t> = Q bi (t) 


(3.76) 


where T 7 bi represents normal coordinates in the ith mode. 

Here is the undamped natural frequency corresponding to 
the ith normal mode. D^. is the modal damping coefficient and Q bi 
is the generalised force. The damping term is allowed to remain 
coupled as no assumption has been made regarding the nature of 
damping. 


It may be mentioned that vehicle body performs two rigid body 
motions in heave and pitch which correspond to zero natural 
frequencies. Other finite non zero frequencies are associated with 
elastic bending modes. The natural frequencies and mode shape 
functions for bending displacement have been found analytically 
for a variable section beam. The displacement shape can be 


obtained 

from eq. (3.18) 

with (r+d c ) substituted 

for s 

for 

the 

vehicle 

domain -d£ r s 
c 

(L b ~d c ) where is the 

length 

of 

the 


fuselage and d c denotes the aircraft c.g distance from the tail 


end. 


Modal parameters for generalised force, mass and damping 
can be obtained from eqs.(3.56), (3.57) and (3.59) using fuselage 
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quantities . 

The modal equations for the left and right wing remain same 
as rigid model. However, the generalised forces should be 
calculated with inclusion of fuselage flexural modes in the 
induced force and torque given in eqs.(2.47) and (2.48). 

3.5. 3 Heave-Pitch-Roll Model 

(a) Flexible fuselage in bending with rigid roll:- 

The fuselage of the aircraft undergoes rigid body heave, 
pitch and roll motions along with elastic bending. The partial 
differential equation of motion representing transverse bending 
will have discretized form same as eq.(3.76). The generalised 
forcing term has to be evaluated on account of the rigid body 
rolling which effect the vertical motions of the two rear wheels. 

Similarly, the discretized equations of the left and right 
wing have the same form as in eqs.(3.74) and (3.75) except the 
generalised force calculation will be different for rolling of the 
airplane inducing different reactions in two main landing gears 
carried under the wings. The generalised forces are to be 
calculated using f and r in eqs . (2.55) and (2.56) . 

In track discretized equations, the generalised forcing term 
has to be evaluated for the dynamic load transfered on the track 
by three wheels in eq. (2.57) 

(b) Flexibile fuselage in bending and torsion: - 

The discretized equations for the bending of fuselage body is 
given by eq.(3.76) with generalised force evaluated after . 
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substituting the expression for vertical applied force to take 
torsional deformation into account. The equation for torsion is 


T >ti (t) ’W 11 + <3 - 77) 

i=l, 2 , . . . 

where the TJ t ^ are the generalised coordinates, ^ti. tors ^- ona ^ 

/ 

natural frequencies, and Q^are modal damping coefficients in 

torsional modes of the vehicle body. The vehicle body performs 
rigid body heave - pitch - roll motion and thus first two modes in 
bending represent heave and pitch motion as mentioned earlier 
where as first mode in torsion represent rigid roll. Rigid body 
modes correspond to zero natural frequencies of the flexible 
modes. For the vehicle body modelled as variable section beam the 
torsional mode T^ can be found as 


V r) - T 0 + T lV r+d c> + J 2 k | 0 M ik T k^ (r+d c ):L 


(3.78) 


for 


-d c s r s (L b -d c ) 


where torsional frequency parameters X fc defined as 


*?= “t I br (0> '' G b J b (0) 


(3.79) 


The integration constants Tq and T^ appearing in torsional 
mode shape function can be obtained with application of two 
boundary conditions for torsional motion in eq. (3.78) . The 
generalised modal quantities for force, mass and damping may be 
obtained from eqs . (3.56), (3.57) and (3.59) by substituting 
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appropriate rotational quantities . 

The discretized equations for the wing will remain same, 
except the generalised forces are to re evaluated using the 
expression for f and F from eqs.(2.55) and (2.56) with w Q 
substituted from eq.(2.58). The track modal equations in normal 
coordinates are same as that used in rigid heave-pitch-roll model. 

3.6 DECOUPLING OF SYSTEM EQUATIONS 

The system equations include the equations of motion for 
lumped mass elements and modal equations for flexible elements in 
the vehicle models. Theoretically, continuous elements have 
infinite number of normal modes. However, the higher modes do not 
contribute as much to the response as the lower modes . The number 
of significant modes can, therefore, be truncated to a finite 
size. Let there be n^ number of significant modes in transverse 
bending of vehicle body (including rigid body bounce and pitch 
motion) , n fc number of twisting modes of the vehicle body 
(including rigid body rolling) , n r vertical degrees of freedom for 
unsprung masses, n w and n^ elastic modes of the flexible wing and 
track. Then total degrees of freedom of the system would be 
n=Uk+n t +n r +2n w +np . The system equations then can be represented 
using matrices as 

• • • 

M q ( t ) + C q ( t ) + K q(t) = F(t) (3-80) 

where q(t) is the response vector, F(t) is the generalised force 
vector (nxl) and M, C, and K are system mass, damping and 
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stiffness matrices (nxn) respectively. 

Premultiplying both sides of the eg. (3.80) by M -1 and using 
an identity q(t) - q(t) =0, the eq. (3.80) can be cast into 2n 
form as 

p(t) + A p(t) = p(t) (3.81) 

f q(t) ] r m" 1 f ( t ) ) 

where p(t) = 4 \ • p(t) = \ - L and 

l q(t) J l o J 



(3.82) 


I being an identity matrix and 0 a null vector/matrix. 

Let os be the eigen values and the eigen vectors for the 

matrix A. Consider a linear transformation of the response vectors 

v(t) = U _1 p(t) (3.83) 

in which U is the modal matrix with eigen vectors arranged as 
its columns. It is now possible to uncouple [87] the system 
equations (3.81) as 


v.(t) + a. V. (t) = Ri<t) 


i=l, 2 , . . . , 2n 


where 


n 


n 


n 


R ± (t) 


- I S ir p r (t) - Ihz J, m rk F k 
r=l r=l k=l 


(3.84) 


(3.85) 


here u. are elements of the inverse of matrix U and m . are the 
lr 

elements of the inverse of matrix M. 
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The general solution of eg. (3.84) may be written as 


00 

v i (t) = X Qi exp( -as t) + J dS (R^ (w) ) (3.86) 

~00 

where X o ^ are the constants of integration to be determined from 
the initial condition, and H^(w,t) are the transient frequency 
response functions given by [87] 


H i (w,t) = 


30 ) + 


a i 


exp ( j cot ) - exp { 


os (t- 


V 


(3.87) 


It can be seen that as t Q -» - oo, for positive real part of a 
. . . 1 

approaches a limiting value of — . exp (j cot) 

30) + as 

The response in generalised coordinates now can be expressed 
by substituting eqs . (3.85) in eq.(3.86) and using eq.(3.82) 


2n 


Sm (t) =l 1 x oi’Vi exp ( -“i t) 
1=1 1 


2n n n _ 

iJl\i x Z*ir k ^ m rk J V w ' fc) dS(F k (w)) 


(3.88) 


m=l, 2 , . . . , n ; m.^ = m + n 


3.7 RESPONSE STATISTICS 


3.7.1 Mean Response 

Considering the expectation of eq.(3.88) the mean response 
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becomes 


V fc) = E IV" ] 


2n 

(t) X oi U m 1 i ex P < '“it* 

mi 1=1 i 


2n n n 


11 XX XX * p 

iSlN 1 r?i iir k=l” rk J Hi< “ -t> E L dS(Ft 


(u) ) 


1 X 1 = 1 , 2 , . . . , n 


where 


E dS (F k (o)) 


* 

4 « (w) 


(3 .89) 

(3.90) 


* 

in which i±„ (w) is the Fourier transform of the mean generalised 

r 

force process 4 p and defined by the integral 

00 

4* (w) = 2 tt J M F ex P ( " 3 ) dr (3.91) 

k k 

- 00 

The mean generalised force depends on the track profile 
characteristics as well as aerodynamic force and moments . The kth 
element of the mean generalised forcing vector can be expressed as 


(t) 


2 

= I 

p=l 


L< 


A-_ c h (t) + B 
kp s mp s 


kps 


h (t) 
mps 


} 


C k V(t) 


(3.92) 


in which h (t) , li (t) and V(t) can be obtained from 
mp s mp s 

egs . (2.2 0) , (2.21) and (2.30). The coefficients A k , B k and C k 

would be different for each model. Substituting for H^(aj,t) and 
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using eqs. (3.90)- (3.92) in eq.(3.89) 

integration over «, yields 


and performing the 


2n 2n n n 

1 V t) °i=l Xoi "“i 1 eXP ( t r ,Vi J^ir ,.L m rk x ik (t) 


i_l r= l k-l 

m=l, 2 , . . .n 


(3.93) 


where 


X ik (t) = 


2 3 

E E 

p=l S=1 


[ kps lps v T kps 2ps ' ' J j ' k 3 


Au. y-,«„(t) + B v _y o _(t)j ] + C v 3- 0 (t) 


(3.94) 

The coefficients A^., and and the components and 

3" 3 of the integral 1^ for different models are given below 


Heave Model 

In this model, input from the track is at single point for 
which p=s=l. Hence these two subscripts are ommited in the 
expression. 

Coeeficients are given as 


0 

for 

k = 1,3, 4, . . 

K 

for 

k =2 

US 


'"k-n (x c> 




1 > K us f or k = n i +1 ' n i +2 ' * • * ' n 

M , j. 

k-n 1 

where n. = 2n +2 
1 w 


(3.95) 
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The expression for B^. is obtained from the expression of A^. by 
replacing K us with C ug . 


for k = 1,2, n 1 +l,n 1 +2, . . . ,n ; 


°- r 

w9 k-2 0 


k-2 L 


(s L ) 1 - (s L / L wL ) 


for k = 3,4,...,n +2 ; 


(3.96) 


0 

M 9 ? 

k-r 0 


J <-(n,.2) (S R> f 1 - (S R ' l wr) 2 ) dS R- 


for k = (n +3 ) , (n +4 n. ; r=n +2 , n n =2n +2 
' w ' w ' ' ' 1’ w ' 1 w 


The components 9" 2 and are 


(t) = f h (r) exp (-a . (t-r) dr 

n J 1 


l mi 


* -C \ LUJ- 

- ( h 0 * EA J 2 


(3.97) 


(t) = J h m (r) exp (-as (t-r) dr 


l k - 1 


= £ hi ^ J i h. E c. 

i=l ' r=l r=l 


(r+1) J- 


(3.98) 


r 3<t) 


= f v(r) 2 exp {- (t-t) ) dr = 


2 (m-1) 

C 0 J 1 + £ C r J 2 


(3.99) 
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with 


II 

H 

h) 


r 

1 - exp ( 

-&•_ t) 


r 

(-l) k r! 

. r-k 

t 

+ V 

t 

oc . 

l 

+ A 

(r-k) ! 

k+1 

k=l 

a i 


and 


( - 1 ) r — ^ exp(-oc.t) 

a . 

1 


(3.100) 


The infinite limits 

considered finite as ii„ (r) 

F k 


of the integrals 3 r 1 , 3" 2 and 9^ are 

is zero for x< 0 and available only up 


to x s t 

Heave-Pitch-Model 

In this model, there are two wheels. Each of them is located on 
the centre of front and rear axle. Both wheels follow the same 
profile along the centre line of the track although input to one 
wheel is time delayed version of the other. The expressions for 
the coefficients of the mean force and integrals are given 
separately for the models with rigid and flexible fuselage. 


(a) Rigid Fuselage 

The coefficients A, are 

kp 


A kp = < 


0 for k =1,2,5, ... ,n 1 and p =1,2 ; k =3 and p=2; 
k=4 and p=l 


K 


'ul 


u2 


for k=3, p=l 
for k=4,p=2 


' ( 'k-n 1 (x l ) 

mP 

k~n n 


^'k-n 1 (x 2 ) 


K 


Ul 


for k = n. +1, . . . , n and p = 1 


k-n„ 


K 


u2 


for k=n 1 +l,...,n and p = 2 


( n i= 4 + 2n w ) 


(3.101) 
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The coefficients B kp are exactly similar to except K U1 and 

K u 2 are to be replaced by and respectively. The 

coefficient is given by 


C k= 


0 for k = 1, . . ,4, n 1+ l, . . .,n 

L T 
wL 

< L 0 / "SV [ I { "£-4 < S L> -wL< S L>^t4< S L>} 


{ 


2 / „ 2 


a 1/2 


wL 


1 " s £ 7 L wl| ds L 


* (M 0 7 M k-4> J C h (S L> # k-4 (S L> ds l 


0 


for k =5,6, . . , m . ( m, = 4 + n ) 

1 i, - w 


wR 


(L 0 / M 


k-m> [ | { W k E mi 


{ 


1 s R / ^ , w p j r ds n | + (M 


A 1/2 


(S R ) " e wR (s R^k-m 1 (S R^ 


wR 


3 r] + (M 0 / Cm^ J C h (S R> ^k-m^V ds 


for k = m^+1, . . . , n^ 


R 

(3.102) 


The integral components 3^ , and 7^ are 


l k 


W t] - [ h 0 < c i,0>p ) J 1 < C i.r>p J 2 < 3 ' 103 


l k - 1 


r 2p (t) ( c iil)p J, (c i,r + l>p (r+1) J 2 


(3.104) 


r 3 , and J 2 are given by eqs . (3.99) and (3.100) 


(b) Flexible fuselage 


The coefficients A^. , and are as 
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0 



for k=l,2 / ...n b , n b +3,...,n b +2+2n w and p=l,2; 

k= n b +l and p=2; k=n b +2 and p=l 

K for k=n b +l and p=l; k=n b +2 and p=2 
np n 

^k-n (x p } ^p/ M k-n for k = n ! +1 '*--' n and P=l/2 

(n = n,+ 2+ 2n ) 
in w 

(3.105) 


The coefficients are exactly similar to A Vr ^ except the 


kp 


stiffness term K is to be replaced by damping term C. The 
coefficients are 


C k= < 


0 for k =1,2,..., n b +2 ; n +1, . . .n 


wL 


VVJJ 

<V M k-r>[ J {"tr (s L> - e „L (s L>fk-r (s L>} ^ 


wL 


+ <V M g-r > K (S L> *k-r< s L> ds l 


for k =n b + 3 , . . . (n b + 2 +n w ) and r= n b +2 
L wR 

‘V M k R r>[ J {’’E-r'V - e wR (s R ) *k-r <s R>} f 1 -^ ^j^R 

L _ 

wR 

+ <V Mf r , J c £ (Sr) # R_ r(SR , dSR 


for k - ( n b +n w +3) , . . . ^and r= n fa + n w +2 (3.106) 
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The integral 5 r 1 , 3" 2 and V are same as eqs (3.103), (3.104) and 

(3.99) 

Heave-Pitch-Roll Model 

In this model, each wheel follows different path along the 
track. Due to variation of longitudinal and transverse roughness, 
input to each wheel is different. Coefficients of mean forcing 
function and the integral for the mean response for two types of 
models are given below. 

(a) Bending of fuselage with rigid body roll 

0 for k=l, 2 , . . .n b +l,n b +5, . . . ,n b +4+2n w for all p, s 

K for k=n, +2 and p=l, s=l; k=n. +3 and p=2, s=2; 

ups b £> 

k=n b +4 and p=2, s=3 

ih . (x ) K / M, for k = n, +5+2n , ...,n and either 
Hc-j v p ; ups ' k-3 b w 

p=s=l or p=2 , s=2 or p=2, s=3 

(j=n b +4+2n w ) 

(3.107) 

The coefficients ®kps are exactly similar to A^p S except the 
stiffness term K is to be replaced by damping term C. The 

coefficients are 



87 





0 for k =1,2, . . .,n, +1, n h + 5 + 2 n , . . . , n 


wL ) ( 

<V M k-j>[ J {^-j (s L> - e wL (s L^k-j (S L>} I 1 " 


2 o ^ 1 / 2 
S L 7 L wL| ds L 


wL 


+ (M 0 / M^j) | c£(s l ) ^-j (S L } dS L 


for k =n b +5, . . . ,n b +4+n w ; j=n fa +4 


wR 


( V M k R j ) [ J {"k-j vo R' * c wR v °R ( r k-j VkJ R 
0 


1/2 


W? 4 (s„) -e, lB (s D ) jl-Sg/ 1- ds 


R 


J wR 


+ <V M k R r> J C h (S R> *k-j (S R> dS R 


for k — n^+n^+5, . . .,n-^+4+2n^; 2=nj^+4+n^ 


(3,108) 


The components of integral 1^ for the mean response are 
given as 

l Ik 

^lps" [ {h 0 } ps + i ^ 1 (h i ) ps (c i,0 ) ps ] J 1 + i J 1 (h i ) ps r^ps J 2 


( 3 . 109 ) 


k -1 


r 2ps i } ps (c i,l } ps J 1 + i J 1 (h i ) ps r ^ 1 (c i,r+l ) ps (r+1) J 2 


where 7^ and and J 2 are same as other models . 


88 



(b) Bending and torsionally flexible fuselage 

( 0 for k=l, 2, . . • n j :i +n t , n^+n t +4, . . . , n^+n t +3+2n w for all p,s 


A kps ’ 


K ups ^ or k =n b +n t +1 an< ^ P =1 » s=1 ; k=n b +n fc +2 and p=2, s=2; 


k=n +n, +3 and p=2, s=3 


^k-j (x p } K ups 7 M k-j for k = n b +n t +4+2n w n and 


either p=s=l or p=2, s=2 or p=2, s=3 

j =n, +n, +3+2n 
J b t '■ 


w 


(3.110) 


The coefficients B kps can be obtained from eqs. (3.110) simply by 
replacing the stiffness term K by damping term C. The coefficients 
C k are 


C k= 


0 for k =1, 2, . . . ,n b +n t +3, n b +n t +4+2n ur , . . . , n 


w 


<V H k-j> 


wL 1/2 

'K-j (S L>- 5 W S L ) *tj (s L>} {^L 7 L S} ds L 


wL 


(Mq/ C h (S L ) ^k-j (S L } dS L 


for k =n b +n t +4, . . . ,n b +n t +3+n w ; j-n b +n t +3 


<V M k R j } 


wR . \l/2 

"£-j< s R> - e wR (S R ) # k-j (s R ) }{ 1 ' S R / K\ ds R 


wR 


w, 

{ 


(M 0 7 M k E j>l c h (s R )# k-j (s R )ds R 


for k =n b +n t +4+n w , . . . ,n b +n t +3+2n w ; j-n b +n t +3+n w 

(3.111) 
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The components of integral for the mean response are 


given by eqs . (3.109) 


3.7.2 Response Covariance 

Covariance of the response for the system generalised 
coordinates can be written as 


K „ (t-, t 9 ) = E 

‘Ji'Jk 1 2 


’ { W - "q. (t l> } { 3k (t 2> - U q k (t 2> } _ 


(3.112) 

Upon substitution of eqs. (3.88) and (3.89) in eq. (3.112) and 
after simplification yields 

2nn2nnnn 

K q.q ^ t l ,t 2^ = .^ ^-i U i ^ U £r mrs U k ^pg m gw 

^x^k 1=1 r=l p=l s=l g=l w=l n n 3 


X { I s ; g (t l ,t 2 ) I £s (t l ) I pg (t 2 ) } 


(3.113) 


where i= i+n, k = k+n, I. and I are defined in eq.(3.94) and 
n n Is pg ^ 

integral I g ^ is given by 

00 CM 

I s,g^ t l ,t 2^ = J J H p^ W 2 ,t: 2^ E [ ds ( F s ( w ]_)) 

~ CO “00 


J jH^(w 1 ,t 1 ) H p( w 2 ,t: 2^ $ 


F F„ (< V W 2 } dW ! dw 2 
s g 


(3 .114) 
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The evaluation of the above integral requires the 
description of the generalised PSD of the input process 
F. It is assumed that aerodynamic forces are uncorrelated with the 
track unevenness. Thus input PSD can be expressed in terms of 
track profile PSD only. 

Various expressions for track PSD have been used in 
literatures and most of them are close to the measured value in 
the range of interest of the spatial frequency. To illustrate the 
present approach, the track roughness is assumed as a homogeneous 
process in space domain with the PSD [73] 

(£2) = A exp (-k r Q 2 ) (3.115) 
n R 

where A r and k r are roughness and correlation constants for a 
particular class of track. The corresponding generalised PSD can 
be put as 


Vr 1 " 1 ' 


n 2> 


$ 


h R (n i } 


5 (Q. 


fi 2 ) 


(3.116) 


With the description of input PSD in terms of track 
characteristics, the integral I s can be evaluated in complex 
domain using Cauchy' s residue theorem. The input PSD and integral 
I for different models are given below 

s,g 


Heave Model 

The input PSD are 
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The integrals for base values other than 2 can be obtained from 
equations (3.117) and (3.119) 

Heave -Pitch Model 
(a) Rigid fuselage 
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The input PSD for the two point input rigid model are given as 
I 0 ; for s or g =1,2 and s * g 


C ul U 1 U 2 + ^ C ul K ul ((J l" W 2 ) + K ul } f h 1 h 1 (a) l' W 2 ) 


for s= g = 3 . 


C ul C u2 (J l W 2 + ^ ( C ul K u2 W l " K ul C u2 U) 2' 


+ K ul K u2 } \h 2 (w l' W 2 J ; f ° r 5=3 ' 9=4 


F F 
s g 


u2 u2 1 2 


C u2 W 1 W 2 + ^ C u2 K u2 (c V W 2 } + K u 2 } $ h 2 h 2 ((J 1 ' W 2 } ; 

for s= g = 4 . 


( 1/M p, { V n x ( X l } f F 3 F 3 (£J l' W 2 } + 


Vn (X 2 } $ F F (w l' U 2 ) j f ° r S=3 ' 9= n l +1 '- 
X 3 4: 

n- =4+2n„ 

1 w 

< ^Vg-V { VV x i> s f 3 f 4 <“i' “ 2 > + 

fg-n < x 2> *P,P <“!' “2> } for s - 4 ' S= n l +1 '- 


p,s-n. p,g-P. 


[ ^s-n-^l 5 ' / 'g-n 1 (x l ) $ F 3 F 3 (w l' w 2 } + 


^s-n_ (x l } ^g-n {x 2 ) + Vn (x 2 } (X 1 } H 3 F 4 (W 1' W 2 } 


Vn^V Vn^l 1 $F 4 F 4 (Wl ' " 2> ] ^ S = ^ n 

and g=n 1 +l, . . . ,n (3.120) 
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Ln 

gj_j}£,g the heave pitch model, two wheels follow the same patii 
along the centre line of the track surface and h 1 and h 2 abre 
representatives of the same process, it can be tacitly assumed 
that 


W" 1 ' 


w 2 } 




“ 2 > = 


V 2 (Ul ' 


U 2 ) = 


V> (<d, / 

^R h R 1 


CO, 


• g 

With the PSD of the input process known, the integral I g 
can be evaluated in complex domain using Cauchy's Residue theorem. 
The integral for the base value s = 3, g =4 is given as 


I3 , 4 ( t^ , t 2 ) — 


2 * A r V l 

v 


r 

f a £ 1 *. _ , 1 

I exp 

{- v r ivr ^1 } 


. 2 
k„a„ n 


{ 


, x 1 f r£ ) 

exp ( -a^t 1 - a p t 2 ) J exp ^ j 

} 


v r K ul K u2 + a £ ( K ul C u2 " v r C ul K u2 } “ C ul C u2 a £ 


a £ + v r a p 


(3.121) 


iy 

The other two basic terms I 3 3 and I 4 4 have forms closely 3 
related to I 3 4 and can be obtained by simply changing subscript^ 
of the damping and stiffness coefficients. I 3 3 is obtained frontf 
eq. (3 .l^i) by replacing the subscripts 2 of K u and C u with 1 every u 
where, while I 4 4 is given by replacing subscripts 1 of and C u 
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with 2 . 


(b) Flexible Fuselage 

The input PSD for this model are expressed as 


0 for s,g = n^+3,..., n^+2+2n w 

and s s g 

{ C ur C uk“l“2 + JSeVV “2> + K ur K uk}*h r h^“l' u 2> 

for s=g; s=n^+l, r=k=l; s=g, s=n^+2, r=k=2; 
s=n^+l, g=n^+2 f r=l, k=2 




(x ) _ (u 1 f w 0 )+ 

A 1 X *s *j+l 


g-nr-l / *F ey F. 4j .. 1 ' w l' w 2' ■ *'g-n:“2'-F_F 


(x„) #_ „ (o> ,u ) 

1 2 F s F j + l 1 2 


"F F 
s g 


for s = n b +l, n b +2 and g=n i +l,...,n 


'^"p.s-n^p.g-n^ [' i s-n< x l ) '"g-n, (x l> *F .P . > 

+ {' # s-n i (x l ) '"g-n i <X 2 )+ Vn/V Vn, (X l’ }*F jFj +; <“l 


+ ip (x„) \p_ _ (x,)$„ F (w-1'^2^ 
^s-n^ 2' v g-n a 2 F j+1 J +i 1 ^ J 


for s = n^+1/ n i +2 , •••, ft 

and g=n i +l ; n^+2, n ' 3 = n £> +1 


(3.122) 


95 



The integral I s,g 
Cauchy' s residue theorem, 
evaluation are 


is evaluated in complex domain using 
The basic terms needed for covariance 


I s,g (t l' t 2 > * '“r’l'Vh { ' “r 1 V*' ^ 1 } 


exp ( -a g t 1 - tt p t 2 ) 


exp (k r <Xj / I | v r K ur K uk 


* “I < K ur C uk - Vur K uk> ‘ C ur C U k “? } ' <V V r“p> (3 ' 123) 


for s=g, s=n b +l, r=k=l; and s=g, s=n b +2, r=k=2; 
s=n b +l, g=n^+2, r=l, k=2 

The integral for other base values can be found with the help of 
eqs. (3.122) and (3.123). 


Heave- Pitch-Roll Model 


(a) 


Bending of fuselage with rigid roll 
The input PSD for this model are 
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(p-p p (w-|^/ u >2 ) 0 foie s , g — 1,2,... n b + 1 i n- b + 5 , . . . , n^+4+2n w 

s g 

= [ C urk C ur'k' u l u 2 + j { w l C urk K ur'k' ' u 2 I Sirk C ur / k / } +K urk K ur'k'] 
h rk h r'k' 1 2 


for s =g =n^+2 , . . . , n^+4 

s=n b +2, r=r / =k=k , =l; s=n b +3 , r=r' =k=k' =2 and 
s=n b +4, r=r / =2, k=k'=3 

for s*g, s=n b +2, g=s+l, r=k=l, r'=k'=2; s=n fa +2, 

g=s+2, r=k=l,r'=2, k'=3; s=n, +3 , g=s+2 , r=k=2 , 


r'=2, k'=3 


' (1/M p,g-k >[ * 9 -k (x l> f F P<“l'“2> + Vk (X 2’ {* 


F F. (cJ l' W 2 ) 
s 3+1 


+ $ F F. (U 
j+2 


1 , 0 ) 2 ^} 


for s=n b +2, . . .,n b +4; g=n b +5+2n w , . . . , n; k=n b +4+2n w ; j=n b +2 


(1/M. 


p,s-k M p,g-k } [^s-k (x l } ^g-k (x i } $ f ,F . (tJ l' u 2 ) 


+ K-k (X 1 } Vk (X 2> + *s-k (X 2 > Vk (X 1 } } { $ F. F.^l'^ 


+ *s-k (X 2> * S -k (X 2* { 


+ 2 V F. + *F. .„F,. '"l'"2> 


F j+l F j+2 1 Z 


j+2 3+2 


for s,g = n, +5+2n , . . . ,n 


(3.124) 


97 



The integral I has the following basic terms needed for 
covariance calculations 


I (t- , t ' 
s,g 1 2 


(27iA r V 1 / v r ) 



a , 


v t, - 

r i 



exp ( -a f t. 


- at) 

P 2 


exp 


,, 2 
(k r“« 


/ v‘) 


{ 


v r*Sj.rk*Sir' k' 


+ a. 


^ K urk C ur'k' v r C urk 


K ur'k' } - 


C urk C ur'k ,a £ 


I 


/ 


( v 


v a ) 
r p 


for s=g, s=n b +2 , . . . , n fa +4 

s=n b +2, r=k=r'=k'=l/ s=n b +3, r=k=r'=k'=2 and 
s=n fa +4 , r=r'=2, k=k'=3 


(3.125) 


The integral for other base values can be found with the help of 
eqs. (3.124) and (3.125) 


(b) Fuselage with bending and torsional flexibility 
Input PSD in terms of track characteristics are 
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<P F F (u 1 , oj 2 ) =0 for s,g = 1, 2 , . . .n b +n t ; n^+n t +4 , . . . , n j-, +n t +3+ 2n w 
s g 


C , C ,, , on + j iu>_ 
urk ur'k' 12 J | 1 

h (a V (J 2 ) 
rk r'k' x 2 


b t ' b t 
C urk K ur'k / _w 2 K urk C ur 


'k' } 


+K . K ,, , 
urk ur k' 


for s =g =n b +n t +l, . . . , n b +n t +3 

s=n K +n t .+l, r=r'=k=k'=l; s=n b +n t +2 , r=r' =k=k' =2 and 


To t 

s=n^+n t +3, r=r'=2, k=k'=3 


for s*g, s=nj_ ) +n t +l, g=s+l, r=k=l, r'=k'=2; s=n^+n^+l, 

g=s+2 / r=k=l,r'=2, k'=3; s=n^+n t +2 , g=s+2 , r=k=2 , 
r' =2 , k' =3 


= (1/M, 


k )| ift 1 ( X— ) $ TTi (tl)f f (Ja y ~r yj . \ ^ x p — \ w-. , u/ « 

L 9 -ic 1 1 z 9 -JC z l 1 z 


) + V (X,) jfp w (w n ,Wo) 


s j + 1 


+ F . ^1 

s 3 + 2 


'"2>} 


for s=n^+n t +l, . . . ,n^+n t +3; g=nj 3 +n t +4+2n w , . . . ,n; 


k=n b +n t +3+2n t-i ; jr^n^+n^+l 


w 


b t 


(1/M 


p,s-k M p,g-k ) [/s-k V 1 Vk (x l ) *FjFj (w l ’V 


* {*s-k <x l> Vk <x 2> * Vk (x 2> Vk V 1 } {* F j p j + i“ 1 ’"2 


+$ F.F. ^ W 1 
3 3+2 


'"2>} + Vk (X 2 ) Vk ‘V { + 

W } 


2 „ (to-,,0)^) + $p _ p 


*p F. v 1' 2 J 
3+1 j+2 


j +2 j +2 


for s,g = n t +n t +4+2n w , . . . ,-n 


(3.126) 
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The basic terms of the integral ^s,g 
calculations are 


nedded for covariance 


I s ,g<tl-t 2 > * <277A r V 1 /V r ) [exp { - - 


- V t. - t_ 

v 1 r 1 2 

r 


- exp ( -a £ t 1 - a p 


t 2 ) ] exp ( / vj ) { v r K urk K ur , k . 


+ oc i ( K ur ] c C ur'k / v r C urk K ur'k'^ C urk C ur , k /CX £ } ^ v r a p^ 


(3.127) 


for s=g, s=n^+n t +l, . . .,n b +n t +3 

s=n b +n t +l, r=k=r'=k' =1; s=n b +n t +2, r=k=r'=k'=2 and 
s=n b +n t +3, r=r'=2, k=k' =3 


The integral for other base values can be found with the help of 
eqs. (3.126) and (3.127) 

Using expressions for the coefficients for mean forcing 

function and integral 1^ and input PSD function with integral 

I for covariance function, displacement mean and covariance of 
s , g 

the system generalised coordinates can be evaluated in different 
vehicle models. In these two expressions only lower half of the 
matrix U are used. Utilisation of the upper half of the matrix U 
would yield the velocity response mean and covariance . The 
acceleration mean and covariance are to be determined by time 
differentiation of the corresponding velocity characteristics . 
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CHAPTER 4 


NONLINEAR VEHICLE MODEL 

All mechanical and structural systems have some 
nonlinearities that arise in various forms. In vehicle systems 
these may be due to suspension, structural stiffness and damping 
nonlinearities as well as due to friction in the suspension. 

4.1 BRIEF DISCUSSION ON THE AVAILABLE METHODS 

Exact solution of randomly excited nonlinear systems, is not 
always possible except for some cases modelled by a Fokker-Plank 
equation with specific conditions. To handle practical problems, 
various approximate methods have been developed. These may be 
broadly grouped as equivalent linearisation, perturbation and 
Monte Carlo simulation 

In the equivalent approach a linear model is developed for 
the nonlinear system based on some performance criterion. Energy 
dissipation per harmonic cycle between nonlinear system and an 
equivalent linear system may be balanced to obtain linearised 
parameters . The linearised parameters can be used in the linear 
model for response evaluation. This method is suitable for 
systems subjected to stationary input. 

In statistical linearisation, the nonlinear system is 
replaced by an auxiliary equivalent linear system. The replacement 
is made so as to minimize the response error between the original 
and the auxiliary systems. Although statistical linearisation has 
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high degree of flexibility in its application to stationary and 
nonstationary excitations, it has some basic limitations, in this 
approach, the probability density function of the response has to 
be assumed a priori. In most applications, it is taken to be 
Gaussian to simplify the mathematics . When the response is not 
having the assumed probability structure, it becomes ineffective 
in providing the correct estimates of the response statistics. 

The perturbation technique needs a power series expansion of 
structural forces with nonlinear terms scaled by small parameters. 
It works well in first order application to small degrees of 
freedom. The technique becomes extremely tedious for system with 
larger degrees of freedom and for orders higher than one. 

An alternative approach for estimating the response 
statistics of randomly excited nonlinear systems is based on 
random digital experiments, popularly known as Monte Carlo 
simulation. The theoretical basis of the approach is that the 
stochastic environment of the system can be interpreted as a 
sufficiently large set of samples confirming to the known 
statistical description. Random environment may consist of 
randomness in the input or system parameters. During experiment, 
the input sample environment is used in deterministic form and 
corresponding response is obtained by structural analysis to give 
a response function. Repeating the experiment generates an 
ensemble which is analysed to provide the response statistics. 

Monte Carlo simulation, in general can tackle the practical 
problems in many real situations. However, it becomes less 
attractive in terms of computer storage and processing time. 
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4.2 VEHICLE-TRACK MODEL 

A single point input model of the aircraft is considered (see 
Fig. 4.1) . Vehicle mass is assumed to be a sprung mass lumped at 
its centre of gravity. A part of the landing gear mass, wheel and 
tyre is modelled as lumped unsprung mass. Shock absorber 
characteristics are treated as nonlinear with combination of air 
spring, orifice and Columb friction damper. Tyre stiffness and 
damping are assumed linear. Two wings are modelled as cantilever 
beams clamped to the fuselage. The mass, stiffness and damping of 
the wings are function of spatial coordinate along the span. 
Runway pavement is modelled by a uniform beam with surface 
unevenness and supported by elastic and dissipative subgrade along 
its length as in the linear model. 

4.3 SYSTEM EQUATIONS 

All displacements are constrained to be in the vertical 
direction and assumed positive upward. Sprung and unsprung mass 
vertical displacements z 1 and z 2 are measured from their 
respective static equilibrium positions. Deflection of wing 
relative to sprung mass displacement along the span is w(s,t). 
Track centre line deflection along the longitudinal axis x at any 
time instant t is denoted by y(x,t) . 


103 




Fig. 4.1. Nonlinear Aircraft Model 
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The equations of motion of the sprung and unsprung masses are 


M S Z 1 + F S (Z 1' Z 2 ' z i' z 2 } “ si" 


E wL I wL (s L ) 


8 2 w t 


as: 


s T =0 


s 

ss 


S 2 w. 


R 


E „ I _ ( s_ ) 
wR wR R 


R 


as' 


(4.1) 


R S R =0 


M us Z 2 - F s (z l- Z 2 ' z l- z 2> - C us< z 2 - y (x o' t> - h(x c> > 


K ug { Z 2 " y {X n' t] ~ h(X r- } ) =0 


(4.2) 


where F is the shock strut force. This force is due to the 
s 

compression of enclosed air in the pneumatic chamber of shock 
strut providing a spring effect (F ) , discharge of compressed oil 
from hydraulic chamber through an orifice (F^) and friction 
between sliding parts (F f ) producing damping. Thus 



F + 
a 


where F 
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V. 


P A 
ao ao 


ao 


V ao~ 


A s 
ao 


n 


(4.3) 


F d - C nL I z I * + -7- “s F a <4 - 4) 

s 

where s= z--z 0 is the shock strut stroke, fe is the stroke 
x z 

velocity, A & is the pneumatic area, P & is pneumatic pressure, 
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is the air volume (o refers to initial condition for the fully- 
extended strut) , n is polytropic exponent for air compression, 4 
is the coefficient of seal friction and C nL is the hydraulic 
damping coefficient which is given by 

C nL = p A h ^ 2 C d A n (4 ' 5 > 


where p is the specific mass of the fluid, is the hydraulic 
area, A n is the orifice area and C d is the coefficient of 
discharge . 

The oscillation of the flexible tapered wing can be 
represented by partial differential equation as 



v„ (s) 


a 2 w \ 
as 2 J 


+ m (s) 
w 



c , . ( S ) 
w 


aw 

at 


- f w ( s , t ) 


( 4 . 6 ) 


where the impressed force per unit length of the span is given by 

f w (s,t)= -m w (s) Zl - c w (s)z 1 -| M k ( w ( s , t ) + Z][ ) S(s-s k ) + V w (s,t) 

k— - 1 

( 4 . 7 ) 

where M k is the kth concentrated mass at span location s k and V ac 
is the aerodynamic lift force. The aerodynamic lift force can be 
expressed as a function of vehicle forward velocity, exposed 
surface characteristics and aerodynamic coefficient as discussed 
in chapter 2 . 

The track is subjected to wheel load that varies with space 
and time. The equation of motion of the track beam elastically 
supported by disspitatve subgrade is given by 
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(4.8) 


a 4 y a 2 y dy 

Vpi + + c f-^- + k f y - V x ' tJ 


where the impressed load per unit length is given by 


f p (x,t) =- C us j z 2 -y(x,t) -h(x)| +K us |z 2 - y (x, t) -h(x)| 


5 (x-x c ) 


(4.9) 


4.4 DISCRETISED EQUATIONS FOR CONTINUOUS SYSTEM 

The wings and tracks are continuous members with distributed 
parameters . The partial differential equations of motion can be 
discretized in terms of normal coordinates by adopting mode 
superpostion technique described in chapter 3 . 

Assuming wing response as a superimposition of the normal 

modes 

00 

w ( s , t ) = £ W i (s)T] i (t) (4.10) 

i=l 

where W^(s) is the ith bending mode shape and 77 ^ (t) is the 
corresponding generalised coordinates . 

Using eq. (4.10) in eq. (4.6) and invoking the orthogonality 
property of the normal mode, yields the modal equations for the 
left and right wing as given in eq.(3.69) and (3.70). 

Following a similar procedure, the discretised equation of 
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the track is found and is given in eq.(3.71) for the linear heave 
model with single point input. 


4.5 SOLUTION APPROACH 

The set of coupled ordinary differential equations (4.1), 
(4.2), (3.69), (3.70) and (3.71) have nonlinearities and can be 

solved by some numerical integration scheme. Of the continuous 
members only first n bending modes of each wing and n for the 
track have been retained in order to limit the size of modal 
equations. The integration of the system equations requires the 
track input. 

The random track profile is generated following Shinozuka's 
algorithm [83] . The track unevenness is homogeneous process in 
space. The zero mean random process is simulated in space using 
its statistical description as 


h R (x) = 


a (2/N) 


1/2 


N 

E cos 

k=l 


(Q k x + 


f k ) 


(4.11) 


where <x is the standard deviation of the track roughness h R (x), Q k 
are independent random variables identically distributed with the 
density function g(Q)sg(Q ) obtained by normalising $^(Q). The <p k 
are identically distributed independent random variables with 
uniform density l/2?r between 0 and 2n. N, the number of terms 
selected for the series has to be large. 

The deterministic part is then superimposed on this to obtain 
complete profile. The track height experienced by the wheel as it 
moves along the track is known against time parameter with the 


108 



help of eq. (2 .1) 

Integration of the nonlinear equations of motions in time 

domain yields system response and the shock absorber stroke and 

velocity at each discrete time step. These are used in eqs.(4.3) 

and (4.4) to obtain the shock strut force F and F , . At each time 

a d 

step t^, differentiation of eq.(4.3) yields instantaneous spring 
constant as 


K g (t k ) = 9F a / as (4.12) 

The expression for the damping force in equation (4.4) is not 
analytically differentiable in the form used. However, being a 
real process its derivative should exist at all times. To obtain a 
usable derivative with velocity a limited order polynomial is 
fitted locally in the force- velocity relation and instantaneous 
damping constant is determined as 

c s (t k ) = (sf^/ at) (at/ as) (4.13) 

/ 

where F^ is the damping force obtained by a polynomial fit. 

The system is then assumed to have a instantaneously linear 
behaviour at instant t fc that sustains in the interval t=t fc to 
t=t k+1 - This assumption is accurate if the time interval used is 
small. The system then can be represented by a set of coupled 
linear differential equations in the small time interval as 


M q ( t) + C(t) q (t ) + K (t ) q(t) = F(t) 


(t k * t* t 


k+1 


) 


(4.14) 
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where M, C and K is the mass, damping and stiffness matrices of 
the system at the instant t fc . The dimension of the matrix is n , 
where 

n = 2+ 2 n w + n p 

Instantaneous response of the linear system can be determined 
after decoupling of the system equation (4.14) by modal approach. 
Response of the system generalised coordinates may be expressed as 


2n 


q^t) - £ X oi (t k > u m exp<-a.(t fc )(t -t fc ) > 


i=l 


2n n n f 

u m i (t k> S “ ir (t k ) J i rg j H i (»,t)<J8(F g <» 
i_i i r-x y- x 


) ) 


(4.15) 


m=l, 2 , . . . , n; m = m+n 


where X . (t, ) are the constants of integration to be determined at 
Ol X. 

each time step to match the response at the end of the previous 
time step, a^(t^) and u^ (t^) are the eigenvalues and eigenvectors 
of the instantaneous state matrix A 


A(t k ) 


m 1 c(t k ) : M _1 K(t K ) 

-I 1 0 


(4.16) 


Further, u^ r are the elements of the inverted modal matrix U 
formed by eigenvectors u^ arranged as columns. EL (w,t) are the 
transient frequency response functions for the linear system in 
the time interval, given as 

H i (u ' t) = J(A« ■ ( t" ~ y { ex P ( ^ t} -exp(- ai (t k ) (t-t k ))} (4.17) 
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4.5.1 Response Statistics 

The mean and covariance of the response can be found by 
considering the expectations of equation (4.15) . Mean response may 
be expressed as 


<^ t) “ta i'V exp(- a .(t k )<t-t k >) 


+ iSl U % i<tk ’ £ «rg 


(4.18) 


where I ±g (t) = A g I 1 (t) + B g I 2 (t) + C I 3 (t) 


(4.19) 


Expressions for the coefficients A , B and C have been qiven in 

9 9 9 

eqs.(3.95) and (3.96). Components of the mean response integral 
1^, I 2 and I 3 are given as 


W = (h o + i J 1 h i c i,0 )J l t i J 1 h i r ? 1 c i,r J 2 

1 2 <t k ) h i c i,l J l + .£ h i E c i r +l <r+1) J 2 

1=1 1=1 r=l 


2 (m-1) 

W J 1 + A c i J 2 


” here J i - ex p(-“i (t k ) (t - t k )) } 


(4.20) 


(4.21) 
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r r 


2 a. 


2 <-» p r! 


p=l 


(r-p) ! p+1 


- r ~P 


oc 


{- 


a ■ 


I ( - 1 )P -G^rrr —5a tiT 5 } ex P<-“i<t-t k )) 

p=l 


(4.22) 


a . 
1 


and 


r 

c i,0 = a 0 ; c i,r= (1/ra 0 ) v ^ ki " r+k)a k c i,r-k for r^l, i=0,l,. 

iC — _L 


• ,1 


c 0 = (a 0 )2; c r =(1 / ra n } ^ (3k-r) a^c;_ w - a£«(k+l)a 


k=l 


k r-k' “k~ “k+1 

r^l and k=0, 1, . . . ,m 


The covariance response of the system generalised coordinates 
may be expressed as 


2n n 2n n n 


n 


q q f 2 ^ ~ X E E E E E u r u.» m i u u ,ra 

r s 3=1 1-1 p = i g=i b=l w=l ^ 3 l 9 s n p pb h 


w 


{ I g < b (t l' t 2 ) ~ I jg (t l ) I pb (t 2 ) } 


(4.23) 


where r R =r+n, s n = s+n; r, s=l, 2 , . . . , n 


t k~ t l' t 2" fc k+l 


The integral 1 i s 


I 3,b (t l' t 2 ) = J I VW 


7^ 

H p (cJ 2 ,1: 2 )$ R F (w 1 ,w 2 )dw 1 dw 

§ b 


(4.24) 


Description of the input PSD are required to evaluate the 
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above integrals. These can be expressed in terms of the track PSD 
and are given by eqs. (3.117) and (3.118). Using the track 
generalised PSD in eg. (3.116), the integral I can be evaluated 
in complex domain with Cauchys' residue theorem. For base value 
g=b=2, the integral becomes 


^ 2 , 2 ^ 1 ' ^2 ^ 


2rrA r V 1 exp(k r a^ / vjj 
v r ( “j + Vp> 


[{ 


exp (-a. I v r t 1 -t 2 |/v r ) 


- expC-*.^- « p (t 2 -t k )) 
v r ex P (k r (v r a p-“j )/V l> 


} { K us v r + c us K us a j (1 - v r> - C u“j } + 
{exp(- a .( tl - t k )-« p t 2 ) -exp(-«.(t 1 -t k ) 


vw* 


„2 2 

C us v r a p + 


C K a (1-v ) 
us us p r' 



(4.25) 


in which v r = vyv 2 where and V 2 are vehicle forward velocities 
at the two time instants and t 2 defined in small interval 

tk+i-tk* The integral for other base values can be obtained with 
the help of I 2 2 in eq.(4.25) and expression for input PSD in 
eqs. (3.117) and (3.118) . 

It may be noted that only the lower half of the matrix U has 
been utilised in equations (4.18) and (4.23) to obtain the mean 
and the covariance of the displacement response. Utilisation of 
the upper half yields the mean and the variance of the velocity 
response. Acceleration characteristics have to be obtained by 
differentiation of the velocity characteristics . 
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The above sequence of operations from instantaneous 
linearisation to evaluation of response characteristics are 
repeated from step to step to cover the required time range. 
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CHAPTER 5 


RESULTS AND DISCUSSIONS 

The present study aims at analytically obtaining the response 
statistics of aircrafts to track induced excitations. In different 
aircraft models considered, slender flexible members like wings 
and fuselage have been idealised as continuous elements with 
nonuni formly distributed parameters. Their natural frequencies and 
mode shapes have been obtained with the application of the 
approach developed in Chapter 3. The natural frequencies and mode 
shapes are used in response analysis of different vehicle models 
with linear shock absorber by adopting modal superposition 
technique. Nonlinearities in the shock absorber behaviour have 
been handled by "instantaneous linearisation" of the shock strut 
behaviour . 

Results have been presented in three sections in this 
chapter. Section 1 contains the theoretical and experimental 
natural frequencies for nonuniform beam and comparison with 
results available in literature. This also serves as a validation 
study for the approach developed. Section 2 presents the 
vehicle-track response statistics for different linear models. In 
section 3, the response statistics of the nonlinear vehicle model 
obtained by instantaneous linear approach of Chapter 4 have been 
presented. Monte Carlo simulation studies have also been carried 
out for the same nonlinear models and results are compared with 
that obtained the by instantaneous linear approach. 
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5.1 FREE VIBRATION OF NONUNIFORM SECTION BEAM 

This section presents results for natural frequencies of 
nonuniform slender members in bending and coupled bending- torsion 
mode of oscillations . 

Results are presented for bending vibration using the 
analytical approach proposed in the present study. These are also 
compared with results available in literature. Experiments have 
been conducted with variable section beams having classical end 
conditions. The experimental results have been compared with those 
predicted by the analytical approach as further validation of the 
proposed theory. To give insight into the behaviour of variable 
section beams, a study is conducted on the sensitivity of its 
natural frequency to changes in beam taper and externally added 
mass . 

Results are also presented for natural frequencies in coupled 
bending- torsion mode for uniform and variable section beams and 
compared with numerical and experimental results available in 
literature . 

5.1.1 Natural Frequencies in Bending Vibration 

The analytical results obtained by using the present approach 
have been compared with published results and are presented in 
Tables 5.1 and 5.2. The natural frequency is nondimens ionalised to 
frequency number as (u 2 m Q L 4 /EI o ) 1/2 , where m Q and I Q refer to s=0 
end. The larger end is mentioned first in describing the beam end 
conditions. The ratios of the end widths b(L)/b(0) and end heights 
h(L)/h(0) are termed as taper ratios and denoted by t^ and t^ 
respectively. The ratio bf the concentrated mass to the beam mass 
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is denoted by m^. and the concentrated mass location as a fraction 
of the beam length is denoted by s r . 

Table 5.1 presents the comparison of nondimensional frequencies 
(frequency numbers) of beams with no concentrated mass loading and 
various end conditions. Table 5.2 is for beams with one or more 
concentrated masses . Published results are shown from various 
references indicated in the tables. 

Comparison indicates good agreement between the analytical 
values obtained and published results. The agreement is better in 
lower modes . This may be so because most of the references have 
used approximate methods and the system becomes more sensitive to 
errors in modelling with the increase in modal sequence. 

5.1.2 Experimental Evaluation of Natural Frequencies in Bending 
Vibration 

Experiments have been conducted on variable section beams with 
three classical support conditions clamped-free (Cl-Fr) , 
pinned-pinned (Pn-Pn) and clamped- clamped (Cl -Cl) . A schematic 
diagram for the experimental set up is shown in Fig. 5.1 

The clamped end condition has been simulated by bolting the 
wider end of the beam between two M.S. plates by six bolts of 6mm 
size. The base plate was 1 0 Ommxl 0 Ommxl Omm in size and was welded 
to a heavy frame. The cover plate was 10 Ommxl 0 0mmx6mm in size. The 
pinned end condition has been simulated with the help of a fixture 
carrying two shafts one vertically above the another. Each shaft 
has been mounted with three sets of ball bearings. The lower shaft 
is held in position while the upper shaft can be raised or lowered 
with the help of a screw arrangement . The beam end was 
accommodated between the upper and the lower ball bearings . The 
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Table. 5.1 


Nondimens ional Natural Frequency (w 2 m o L 4 /EI o ) 1/2 of 
linearly tapered beams vith no concentrated mass loadings 


End 

Cond. 

Taper 

Result 

Source 1 

Mode Sequence 

2 3 4 

5 

Cl-Fr 

0 .2, 1. 0 

* 

[43] 

[44] 
[50] 

5.3973 

5 .3969 
5.3976 
5.3976 

25.6555 

25.656 

25.6558 

25.6560 

65.7473 

65.7470 

65.7470 

125.2584 

125.2600 

204.3782 

204.5500 


1.0, 0.2 

★ 

[43] 

[44] 

4.2920 

4.2926 

4.2925 

15.7425 

15.7420 

15.743 

36.8848 

36.886 

68.1160 

68.144 

109 . 9079 

110.06 


0 .2, 0 .2 

* 

[43] 

[50] 

6.1965 

6.1972 

6.1964 

18.3852 

18.384 

18.386 

39.8335 

39.8370 

71.2413 

71.2880 

112.8328 

113.3300 

Fr-Cl 

0.2, 1.0 

* 

[44] 

2.1709 

2.1709 

18.9589 

18.9589 

58.5574 

58.5574 

117.5956 

196.5318 

Cl-Cl 

0.2, 1.0 

* 

[44] 

21.4707 

21.4707 

60.3316 

60.3321 

119.3342 

119.3396 

198.1076 

295.2434 


1.0,1. 2 

* 

[46] 

24.5634 

24.5634 

67.7044 

67.7048 

132.7371 

132.7240 

219.2361 

324.8328 

Cl-Pn 

0.2, 1.0 

* 

[44] 

16.5075 

16.5074 

51.0248 

51.0250 

105.4035 

105.4009 

179 . 0118 

272.0306 


1.0, 1.2 

* 

[46] 

16.5029 

16.5029 

54.4611 

54.4614 

114.0600 

114.0516 

195.1918 

296.6478 

Pn-Cl 

0.2, 1.0 

* 

[44] 

13.4846 

13.4845 

47.7671 

47.7672 

101.8402 

101.8389 

175.7147 

268.8754 

Pn-Pn 

0 .2, 1 . 0 

* 

[44] 

9.6722 

9.6721 

39.5851 

39.5851 

89.0714 

89.0704 

158.2421 

246.7976 


1.0, 1.2 

★ 

[50] 

10.8267 

10.8270 

43.3565 

43.357 

97.5357 

97.535 

173.3681 

173.38 

271.4061 

271.41 


Note: * represents the results from the present work. 
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Table. 5.2 


Nondimensional Natural Frequency (w 2 m L 4 /EI ) 1/2 of cantilever 

" o o 

beams with concentrated mass loadings 


Taper Result Mode Sequence 

t b , t h Source 123 45 


(a) Single Concentrated Mass 



m : 

=0.5, s r =0.3 




r 




* 

3.4513 

17.4415 47.3072 

113.8837 

196.6725 

[50] 

3.4512 

17.430 47.302 

113 .83 

196.35 


m = 

=0.5, s =1.0 




r 

r 



* 

2 . 0156 

16.9036 51.7071 

106.1008 

180.2362 

[50] 

2.0163 

16.901 51.701 

106.06 

180.18 

[47] 

2.0163 

16.901 51.701 

106.06 

180.12 

(b) Three equal Concentrated Masses 



m r = 0 . 01 ; 

S r =l/3, 2/3, 1.0 



★ 

3.4256 

21.3949 59.4744 

118.5629 

194.3881 

[49] 

3.4262 

21.3943 59.4734 

118.5838 

194.4162 


m = 0 . 1 ; 
r ' 

S r =l/3, 2/3, 1.0 



* 

2.8406 

17.5515 47.7382 

109.3044 

171.7031 

[49] 

2.8406 

17.5502 47.7528 

109.2080 

171.6661 


3 

H 

II 

o 

to 
*«■ • 

S r =l/3, 2/3, 1.0 



★ 

2.4460 

15.1708 41.1835 

105.3713 

162.1602 

[49] 

2.4460 

15.1688 41.1555 

105.3552 

162.1501 


Note: * represents the results from the present work. 
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upper shaft was lowered to just make bearing contact with the top 
of the beam surface. This allowed freedom of rotation to the beam 
end with no vertical displacement. The arrangement was housed in a 
rigid frame constructed with 150mmx75mmxl0mm M.S. channel section. 

The beam was excited harmonically by an electromagnetic shaker 
(capacity 62.5 Kg) through a power amplifier. The power amplifier 
had a built in sinusoidal waveform generator which provided a 
variable amplitude, variable frequency, sine wave out put. The 
frequency of the wave could be adjusted from 1 Hz to 2 0 KHz in 
seven over lapping ranges. In the experiment, driving frequency 
was monitored through the power amplifier and was varied to match 
the resonant frequency. Accelerometers weighing 0.0296 kg were 
mounted at planned locations on the beam and were treated as 
concentrated masses. The output signals from the accelerometers 
were passed through a signal conditioner and observed on a 
dual -trace oscilloscope to identify the event of resonance and the 
corresponding frequency was recorded. 

Two different beam specimens made of Aluminum ( p = 2744 
3 

kg/m , E = 70 GPa) have been tested. The beams have linearly 
varying width with constant height. Details of the specimens and 
their support conditions are listed below. 

Specimen B1 ; 

L =0.3 m, b (0) =0.05 m, h(0) =0.00325 m, b (L) =0.02 m. End 
conditions: clamped- free and pinned- pinned. Concentrated 
mass: m 1 = 0.31, s^^ = 0.5 L. 


12 0 



Specimen 

I Accelerometer 



Power ampli - 

tier with Electromagnetic 

Signal gene- shaker 

rator 



Signal Oscilloscope 

conditioner 


r 

Fig. 5 . 1 . Block diagram of experimental set-up 
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0.00385 m, b (L) 


0.04 m. End 


Specimen B2 : 

L =0.60 m, b(0) = 0.1 m, h(0) = 

conditions: clamped- clamped; Concentrated masses: m^nu^m^ =m^ 
= 0.067, s 1 =0.225 L, S 2 = 0.492 L, S 3 = 0.710 L, S 4 = 0.980 L. 

The experimental values of the first three natural frequencies 
are given along with the theoretical values in Table 5.3. These 
show that first three natural frequencies match within acceptable 
limits for clamped-free and clamped- clamped conditions. However, 
for the simply supported condition the matching is not so good as 
in the other cases. The discrepancies may be due to the 
limitations of the experimental set up not being able to provide a 
ideal clamped and pinned end supports used in analysis. The 
disagreement is greater in the pinned ends case as the restraint 
against rotation at the ends was not completely eliminated by the 
fixtures . 


Table 5 . 3 

Experimental validation of Natural frequency 
( t fa = 0.4, t h = 1.0) 


Specimen 

No. 

Support 

condition 

Source 

of 

Result 

Natural 

f l 

Frequencies (Hz) 

f f 

Z 2 3 

B1 

Cl-Fr 

Theory 

27.62 

187.79 

444.69 



Experiment 

27.00 

190.00 

451.00 


Pn-Pn 

Theory 

103.33 

325.00 

770.00 



Experiment 

107.00 

337.00 

790.00 

B2 

Cl -Cl 

Theory 

46.74 

130.05 

250.18 



Experiment 

43.00 

125.00 

241.00 
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5.1-3 Effect of parameter variation 


on natural frequencies in 


bending vibration. 

The sensitivity of the bending vibration to beam taper and mass 
location have been studied by varying these parameters . The 
proposed analytical approach is used to obtain the natural 
frequencies and mode shapes for variable section cantilevered 
beams . Numerical results have been generated with the following 
data : 

p = 2744 kg/m 3 , E = 70.0 GPa, L =0.60 m, b(0) = 0.10 m, h(0) = 

0.00385 m, c(0)/pA(0) = 0.56 s - '!' 

The boundary conditions of the beam are assumed as clamped at s 
= 0, and free at s = L. Two cases of variation of beam width and 
depth have been considered. These variations are chosen to 
represent a linear and nonlinear convex (parabolic) taper. The 
width b(s) at a distance s from the origin, can be written as 

b(s) = b(0) 1 - w 1 (s/L) J (linear) 

(5.1) 

b(s) = h(o) £ 1 - w 1 (s/L) 2 J (parabolic) 

in which = 1 - b(L)/b(0) . 

Similarly the thickness of the beam can be expressed by 
appropriately replacing b(s) by h(s), b(0) by h(0), b(L) by h(L) 
and w^ by w 2 where w 2 = 1 - h(L)/h(0) . 

The polynomial functions f- L (s) and f 2 (s) for the stiffness and 
mass and damping variation along the span are expressed as 
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(i) Linear 

f^s) = 1 - (w 1 + 3 w 2 ) (s/L) + 3 (w 2 + w^) (s/L) 2 

- (w 2 + Sw^ 2 ) (s/L) 2 + w -j . w 2 (s/L) 4 

(5.2) 

2 

f 2 (s) = 1 - (w 1 + w 2 ) (s/L) + w 1 w 2 (s/L) 

(ii) Parabolic 

f (s) = 1 - (w x + 3 w 2 ) (s/L) 2 + 3 ( w 2 + w^) (s/L) 4 

- (w 2 + 3w 1 w 2 ) (s/L) 6 + w 1 w 2 (s/L) 8 

(5.3) 

f 2 ( s ) = 1 - (w 1 + w 2 ) (s/L) 2 + w x w 2 (s/L) 4 

The natural frequencies are obtained for the following forms 
of tapered cross sections - width varying with thickness constant, 
thickness varying with width constant and width and thickness 
varying simultaneously. The effect of these three forms of taper 
on natural frequencies have been studied without any external mass 
loading as well as with one concentrated mass at different 
locations. The eigen values os are complex in presence of damping 
and appear in complex conjugate pairs. The real part of as gives a 
decaying factor associated with each mode and the imaginary part 
gives the corresponding damped natural frequency. 

5. 1.3.1 Effect of taper ratio on natural frequencies 

Nondimensional damped natural frequencies of cantilevered beam 
in first five modes for different combinations of taper are 
plotted against taper ratios and shown in Fig. 5.2 


124 



(a) Width taper, height constant 

Fig. 5. 2 (a) presents the variation of the first five frequency 
numbers with the two types of taper in width only. The fundamental 
frequency of the beam is seen to decrease with the increase in 
taper ratio for linear as well as parabolic shapes . At lower taper 
ratio, parabolically tapered beam has slightly lower natural 
frequency compared to a linearly tapered beam. In the second mode, 
the frequency number gradually falls with the increase in taper 
ratio. The two types of tapered beams show similar behaviour. 
Third and fourth frequency numbers of the beams show decreasing 
values with taper ratios for both linear and parabolically tapered 
beams. Out of the two, the parabolically tapered beam has slightly 
higher values with the difference decreasing with increase in 
taper ratio. The fifth bending mode frequency for the linearly 
tapered beam has a pattern similar to the third and the fourth 
bending modes. Parabolically, tapered beam, however, has a smaller 
frequency number compared to the linearly tapered beam at low 
values of taper ratios. It has a more moderate decrease in 
frequency with change in taper ratio. 

(b) Height taper, width constant 

The variation of first five frequency numbers for the 
cantilever beam with height tapering and constant width is 
shown in Fig.5.2(b). In case of the first mode, the fundamental 
frequency for tapered beams decreases with increase in value of 
the taper ratio. Second, third, fourth and fifth frequency numbers 
increase with increase in taper ratio. The natural frequency for 
the parabolic tapered beam is seen to be higher than the linearly 
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tapered beam for all the modes considered. 


(c) Height and width taper 

The variation of frequency numbers when both width and height 
taper simultaneously in the same proportion, are presented in 
Fig. 5. 2(c). The first mode frequency is seen to decrease with 
increase in taper ratio for both types of beams showing nonlinear 
trends. The second, third, fourth and fifth frequency numbers for 
linearly tapered beam are seen to increase with taper ratios. This 
is also observed for parabolically tapered beams in the third, 
fourth and fifth modes. However, in the second mode, the 
parabolically tapered beam has reduction in frequency with 
increase in taper ratio for small values of taper ratios . With the 
increase in taper ratio beyond 0.4, the change in natural 
frequency is very small for the parabolically tapered beams . 

5. 1.3. 2 Effect of mass location on natural frequencies 

The effect of location of single concentrated mass ( m^= 0.27) 
on the nondimens ional natural frequency number of the tapered beam 
has been shown in Fig. 5.3 for different combinations of taper. 
Only one taper ratio is considered to illustrate the results. 

In Fig. 5.3 (a), the first five frequency numbers of the beam 
tapered in width only are shown for linearly and parabolically 
tapered beams with taper ratio t b = 0.4. In the first mode, for 
both types of taper, frequency increases as the point mass moves 
from tip towards the support. In the second, third, fourth and 
fifth modes, frequency shows a strong dependence on the position 
of mass, but no regular pattern can be observed. Both types of 
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Nondimensional Natural Frequencies 




u 


t h 


tbi th 


Fig. 5. 2. Effect of taper on natural frequencies, (a) Width 
taper (b) height taper (c) width and height taper. 

Key: Linear taper , parabolic taper 
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taper, in general show a similarity in frequency behaviour. In the 
first mode, the natural frequency of parabolically tapered beam is 
higher than the linearly tapered beam when the mass is near the 
tip and is lower with the mass closer to the root . In all other 
cases, for most of the location of concentrated mass, the natural 
frequency of parabolically tapered beam seems to be higher than 
the linearly tapered beam. 

Fig. 5. 3(b) presents the variation of frequency numbers for the 
beam tapered in height only with width constant. The first mode 
frequency decreases as the location of the mass shifts towards the 
tip. Second, third, fourth and fifth mode frequency numbers in 
linearly tapered beam do not reveal any regular pattern of 
variation. The parabolically tapered beam follows similar pattern 
in all the higher modes except in fifth mode. In the fifth mode, 
frequency number of the parabolically tapered beam decreases only 
marginally as the mass is moved towards the tip. For every 
location of concentrated mass, the parabolically tapered beam 
shows higher value of natural frequency in the first five modes 
considered. 

The variation of frequency number of the beam tapered in width 
and height simultaneously in the same proportion is shown in Fig. 
5.3(c) for various locations of single concentrated mass. Linear 
as well as parabolically tapered beams show decreasing values of 
the fundamental frequency with the shifting of mass towards free 
end. Second, third, fourth and fifth frequency numbers do not show 
any predictable pattern of variation with location of mass for the 
linearly tapered beam. Parabolically tapered beam shows more or 
less similar behaviour in second, third and fourth modes. However, 
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Nondimensional Natural Frequencies 
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6ig.5.3. Effect of mass location on natural frequencies. 

(a) Width taper (b) height taper (c) width and height taper. 
Key same as Fig. 5. 2. 

















in the fifth mode, no significant change in frequency with the 
change in location of mass can be observed. 

5.1.4. Natural Frequencies in Coupled Bending-Torsion Vibration 

5. 1.4.1 Uniform Beam 

Results of a thin walled open section beam [69] and uniform 
section cantilever wing beam [70] have been selected for 
comparison purpose. Present theory has been used to obtain the 
results of uniform section beams as a special case of variable 
section beam for which the functions f 1 (s) , . . . , f g (s) in eq.(3.24) 
representing the variation of flexural stiffness, mass, torsional 
rigidity, mass polar moment of inertia and shear centre offset are 
taken unity. 

Table 5.4 presents the natural frequency of coupled bending 
torsion vibration of thin walled open section beam with the 
following characteristics [69] 

L=5 . 0 m, m=2 . 45 kg/m, EI=5.8xl0 4 N.m 2 , GJ=78.3 N.m 2 , I =0. 02 

kg.m, e=0.08 m. 

First five natural frequencies for clamped- free, 
clamped- clamped, pinned-pineed and free- free boundary conditions 
obtained by the present approach are compared with the results 
quoted in Ref. [69] . In the reference both approximate and exact 
results for the uniform beam are given. These results are only 
compared with the results obtained by present theory in Table 
5.4. 
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Natural frequencies of a uniform cantilever wing is shown in 
Table 5.5 along with the results given in Ref. [70], wing 
characteristics are 

L=6.0 m, m=35 . 75 kg/m, EI=9.75xl0 6 N.m 2 , GJ=9.88xl0 5 N m 2 I = 

' a 

8.65 kg.m, e=0.18 m. 


Table 5.4 

Natural Frequencies of Thin Walled Open Section Beams 
in Coupled Bending-Torsion (Hz) 


End 

Cond. 

Source 

of 

result 

1 

Mode 

2 

Sequence 

3 

Number 

4 

5 

Cl-Fr 

★ 

2.41 

7.16 

11.83 

14.95 

21.55 

Exact [69] 

2.41 

7.16 

11.83 

14.95 

21.55 

Approx. [69] 

2.40 

7.13 

11.78 

14.85 

21.27 

Cl-Cl 

★ 

6.06 

12.31 

18.45 

24.81 

31.07 

Exact [69] 

6 . 07 

12.32 

18.45 

24.81 

31.07 

Approx. [69] 

6.06 

12.26 

18.28 

24.40 

30.27 

Pn-Pn 

★ 

6.21 

10.82 

18.27 

22.31 

26.39 

Exact [69] 

6.21 

10.82 

18.27 

22.31 

26.39 

Approx. [69] 

6.20 

10.77 

18.09 

22.07 

26.12 

Fr-Fr 

★ 

11.13 

12.89 

18.84 

24.29 

30.99 

Exact [69] 

11.12 

12.89 

18.84 

24.29 

30.99 

Approx. [69] 

11.08 

12.81 

18.68 

23.87 

30.19 
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Table 5.5 


Natural Frequencies of a Uniform Section Wing Beam 
in Coupled Bending-Torsion (u rad/s) 


Source of 
result 

1 

Mode 

2 

sequence 

3 

Number 

4 

5 

* 

49.6 

96.9 

248.8 

355.4 

451.7 

[70] 

49.6 

97.0 

248.9 

355.6 

451.5 


First five natural frequencies show excellent agreement with 
the published results. 


5. 1.4. 2 Nonuniform beam 

Actual wing of aircraft can be taken as an example of variable 
section bending- torsion coupled beam. Cross sectional properties 
of a wing of a small military aircraft have been measured at 
discrete span location [88] . Power series have been fitted in to 
the data to enable theoretical evaluation of natural frequencies 
using the proposed analytical approach. The following functions 
have been found to fit the actual data [Fig. 5. 4] 

m(s) = m(0) [1-0.39 (s/L) ] 

El ( s) =EI ( 0 ) [l-l.ll(s/L)+0.3 (s/L) 2 ] 

GJ ( S ) =G J ( 0 ) [l-1.21(s/L)+0.33 (S/L) 2 ] (5.4) 

e (s) =e (0) [1- . 75 ( s/L) + 0.27(s/L) 2 ] 

I a (s)=I a (0) [1-0 . 87 (s/L) ] 
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Fig. 5. 4. Polynomial fit of wing cross sectional properties. 
Key: Measured data [88]++ + + + , polynomial fit 








where the nondimens ional parameters at the root section are 

m 0 e 0 2/ x a o = 1 - 23 and EI 0 /GJ 0 = 2 ' 18 

Table 5.6 gives the comparison between theoretical and 
experimental values. The closeness of the agreement indicates that 
the proposed approach may be employed for a general section for 
accurate evaluation of dynamic characteristics. 


Table 5.6 

Comparison of theoretical and experimental values of 
natural frequency (Hz) 


Mode No 

★ 

Theory 

Experiment [88] 

1 

6.87 

6.50 

2 

21.73 

20 . 00 

3 

36.33 

39.75 

4 

55.17 

52.00 

5 

69.07 

67.00 


5.2 RESPONSE STATISTICS OF LINEAR VEHICLE MODEL 

A passenger plane has been selected for the study. The 
following system data have been adopted to generate numerical 
results : 

Aircraft .* 

Total mass (M^) : 1.469 xl0 5 Kg; Pitch Moment of inertia (I ) : 
7.178 x 10 8 Kg m 2 ; Roll moment of inertia: 6.325 x 10 5 Kg.m 2 ; 
Fuselage length: 44.35 m,* Main gear distance from aircraft 

c.g(l 1 ) :1.32 m; Nose gear distance from c.g (1 2 ) : 16.68 m; Main 
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gear tread 6.3 m ; Flexural rigidity of fuselage (E b I b (0)) 9.135 
10 2 

X 10 N. in ; Torsional rigidity of fuselage at origin 

2 

^b^b^^ : 11 * 6 x 1( - 1 N * m - Airframe viscous damping 

(CbtOJ/^CO) ) : 0.298 s~* 

Landing Gears : 

Nose gear mass (M^.^) : 0.248 x 10 2 Kg; Main gears masses 

^u22 ' ^u23^ : ^ * -^-■'■05 x 10 Kg; Nose gear suspension spring 

constant (K g ^) : 2.36 x 10^N/m; Two main gears spring constants 

6 

(K S22' K s23 ) : 14 • 5 x 10 N/m; Nose gear suspension damping 

4 

( c s il) : 8.24 x 10 ; Main gears suspension dampings (C s22 , C g23 ) 

4 

5 0.1 x 10 N.s/m; Nose and main wheel tyre stiffnesses 

(K ull ,K u22' K u23^ :1 - 34 x 1C)6 ' 8 - 5 x 106 and 8 . 5 x 10 6 N /m ; Nose 
and main wheels tyre damping ( C ull /C u22 , C u23 ):0.98 x 10 4 , 7.86 
x 10^ and 7.86 x lO^N.s /m 
Wing: 

Span (L w ) : 20.0 m; Aspect ratio (A^J : 7.08; Mean chord (c^CO)): 

3.0 m; Flexural rigidity (El (0)): 7.174 x 10 6 N.m 2 ; Torsional 

w w 

6 2 

rigidity (G J (0)) : 4.662 x 10 N. m ; Mass moment of inertia 
w w 

2 

about elastic axis (I (0)): 8.643 kg. m /m; Mass distribution 

at the root section (m w (0)): 42.73 kg/m; Shear center offset 

_ 1 

(e (0) ) : 0.25 m; Linear damping (c w (0)/ m w (0)): 0.16 s ; 

Torsional damping (il (0) /I (0) ) : 0.0 08 s 1 ; Aerodynamic 

coefficients for lift and moment (C L , C mac ) : 0.8 and 0.1; Air 

3 

density (p ) : 1.12 kg/m . 

3 . 

Track and Pavement : 

Runway length (L ) : 1400.0 m; Track mass (m Q ) : 3620 kg/m ; Track 
flexural rigidity (El) s 1.38 x 10 7 N.m 2 ; Foundation damping 
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(c / m ) : 0.04 s” 1 ; Foundation spring constant (k f ) : 1.705 x 

v f 0 5 

10 8 N/ m 2 ; Roughness constant (A r ) : 0.5025 x 10 ; Correlation 

index (k r ) : 0.1012; Track mean profile as in eq.(2.15) with h Q 

= h = h = = 0 and h ± = 0.001 + 2TtA Q / W i; h 3 = - 2 ttA 0 /(3! 

W 2 ); h 5 = 2 tiA 0 /(5! W^) etc. where W x = lOOn m and A Q = 0.075 m 

for nose wheel path, 0 . 04 m for left main wheel and 0.06 m for 

right main wheel path. This represents a combination of uniform 

slope and sinusoidal profiles followed by different wheels . 

The response statistics of the aircraft model has been obtained 

for constant velocity, accelerating and decelerating runs during 

taxi, takeoff and landing on the runway. The response selected for 

the study are displacement, velocity and acceleration of the 

system generalised coordinates. 

For the taxi run, three uniform aircraft velocities 4 0 Km/h, 

60 Km/h and 80 Km/h have been considered. Response plots are shown 

for initial 20 seconds period of motion. In the takeoff run, 

responses are plotted for three forward acceleration of the 

aircraft - 1.6 m/s 2 , 1.8 m/s 2 and 2.0 m/s 2 . Time history plots are 

given from the start of forward motion till the aircraft reaches 

lift off velocity at 216 Km/h. In the landing run, time histories 

are shown from touch down instant to the stoppage of the forward 

motion of the aircraft. Responses are presented for sink 

velocities 0.6, 0.9 and 1.2 m/s. The glide velocity and uniform 

2 

slowing rate are assumed to be 215 Km/h and 1.5 m/s . 

The response in the three types of runs are being presented 
for the different models used in the study one by one. 
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5.2.1 Heave Model 


5 . 2 . 1 . 1 Taxi Run 

The aircraft is assumed to move forward at constant velocity 
in the taxi run. This corresponds to a situation of the vehicle 
moving at constant speed over a smooth surface and then suddenly 
encountering a rough terrain at t = 0, without any change in the 
forward velocity. 

The mean displacement, velocity and acceleration of the 
generalised coordinates are presented in Fig. 5. 5. The mean 
displacements of the sprung and unsprung masses and the track 
modes in the later phase of taxi run follow the pattern of the 
mean track profile. The small initial transience gets damped out 
quickly. The magnitude of the mean displacements, during early 
phase of travel, do not show much variation with the vehicle 
velocity. However, with elapse of time, higher velocity seems to 
produce increased displacement. Due to presence of a rising 
gradient in the track mean used for the study, vehicle wheel 
experiences more height as it travels forward and more dynamic 
load is imposed on the pavement . 

The mean velocities of the lumped masses and the track modes 
reveal stronger and more persistent transience. Except for the 
wing normal modes, which show very low values, the other steady 
state velocity responses have the same general pattern as the 
displacements with a difference in phase. The magnitude of the 
mean velocity is higher for higher vehicle forward velocity. 

The mean acceleration of the lumped masses follow a pattern 
somewhat similar to their mean velocity, except for the difference 
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in phase. The magnitude of transient accelerations for all the 
response are higher in comparison to their steady state values. 
This may be due to the large force transmitted to the masses 
suddenly at transition from smooth to rough surface. The wing 
normal mode mean accelerations have been found to be sensitive to 
the vehicle forward velocity after elapse of some time. The mean 
acceleration of the track first three normal coordinates have 
initial transience on entry of vehicle over the track which is 
quickly damped out. 

The displacement, velocity and acceleration mean responses of 
the wing first three normal coordinates have initial high 
transient values, diminishing to lower values at later stage of 
the taxi run. The higher vehicle velocity is found to produce 
larger amplitude for the wing response in its first three bending 
modes . The track input frequency seems to be dominant in the wing 
response in later stage. 

The lumped mass mean response show that the track input 
frequency is dominant in the later stage of taxi run. The input 
frequency is also modified by the vehicle forward velocity as 
seen by the shift of peaks towards left with higher vehicle 
forward velocity. In all the modes, it is found that higher 
vehicle velocity is responsible for higher response. In most of 
the cases, the vehicle forward velocity is seen to modify the 
response frequency in the later part. 

The response variances are presented in Fig. 5. 6. The 
displacement variances show that in early stage, the response is 
oscillatory with a high peak at the onset of vehicle motion. 
Subsequently the response subsides to reach a steady asymptotic 
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value. The steady state covariance magnitudes for displacements 
are higher for higher vehicle speeds . The track mode 
characteristics have low frequency variations in the steady state. 
The variances of the velocity response follow a trend similar to 
the displacements . The acceleration variances indicate dominant 
values at the initiation of motion. The steady state condition is 
reached within shorter time than that for displacement and 
velocity without any oscillations, showing progressively smaller 
peaks. The steady state values, however, are very low compared to 
the initial values. 


Increase 

in 

the vehicle speed is seen 

to 

increase the 

magnitude 

of 

the 

response variances for all 

the 

coordinates . 

Comparison 

of 

response mean and variance in 

the 

first three 


flexible wing and track bending modes indicate that response 
magnitude decreases as the order of mode increases . This conforms 
to the usual assumption that the lower modes dominate the higher 
ones and inclusion of the first few lower modes in response 
computations is adequate . 

5. 2. 1.2 TaReoff Run 

During this run aircraft is assumed to start from rest, 
increasing its speed with uniform acceleration to lift off. The 
time history plots are from start of aircraft forward motion till 
it lifts off. 

Fig. 5. 7 presents the mean response. Initial portion of the 
vehicle response show low magnitude transience. The accelerating 
vehicle takes some time to attain a certain level of velocity for 
the track input to become effective in setting up vibration in the 
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vehicle structure. Subsequently, mean response magnitude grows 
with increase in vehicle speed. Influence of the mean track 
profile is apparent in the vehicle mean response pattern - 
increasing level with sinusoidal variation. The constant frequency 
transient oscillation, as in taxi, is not observed in the takeoff 
run. Mean acceleration response of unsprung mass indicates a 
departure and has high frequency transience in the initial phase. 

The mean displacement responses of the track beam for its 
first three normal coordinates show low initial response that 
grows as the vehicle accelerates to takeoff . The magnitude of the 
displacement is quite low compared to that in constant velocity 
run. This may be as the moving load is passing fast over the track 
beam, the track structure has little time to react to this load. 
The mean velocity has stronger initial transience and its 
magnitude is comparable to the response in the later phase of 
takeoff run. The mean acceleration of the first three normal 
coordinates of the track have high value prior to takeoff . 

The response variances in takeoff run are presented in 
Fig. 5. 8. The variances for the vehicle response show very low 
values in the initial stage. As the vehicle accelerates from rest, 
its response variance show a steady increase till takeoff. The 
magnitude of the variances are seen to increase with the forward 
acceleration of the vehicle. The unsprung mass response variances 
are smaller compared to the sprung mass in displacement but larger 
in velocity and acceleration. 

The variances of the first three track normal coordinates 
have an oscillatory pattern about an increasing mean level . There 
is a gradual increase in amplitude of oscillation till take off. 
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Higher vehicle accelerations induce increased response variances 
of track s first three normal coordinates. The oscillatory trend 
increases from displacement to velocity to acceleration response . 

5 . 2 . 1 . 3 Landing Run 

The response quantities studied during landing run are the 
same as in taxi and takeoff. The shock strut at the time of landing 
being under no load, is fully extended. This is utilised to 
obtain the initial displacements of the sprung and unsprung 
masses from their static values . The initial velocities for 
sprung and unsprung masses are the sink velocity of the aircraft . 
The initial conditions for the wing relative to the sprung mass 
and for the track are assumed to be zero. Response characteristics 
plots have been obtained from the instant of touch down till 
stoppage of the aircraft forward motion. 

Fig. 5. 9 presents the mean displacement response of the 
vehicle-track system. The response pattern shows initial high 
value due to touch down impact, increasing with increase in sink 
velocity. After the initial impact energy is dissipated, the 
steady state response becomes dependent on the track roughness . 
The mean track input frequency is seen to play a dominant role in 
response during the later phase of landing run. As the vehicle 
looses its forward velocity, the effectiveness of input from the 
track decreases, reducing the system response. The magnitudes of 
response during impact phase are larger than that of the taxi and 
takeoff runs. The displacement and velocity means of the wing's 
first normal coordinates show the development of largest peak 
subsequent to impact. The velocity mean for the wing second normal 
coordinate indicates the presence of more than one frequency in 


141 



the oscillation during the impact phase. 

Response variances are presented in Fig. 5. 10. Strong 
oscillatory patterns are present for all the displacements and for 
the vehicle velocity responses in the touch down phase which 
subsequently decrease gradually with the slowing of vehicle 
forward motion. The magnitude of the first peak seems to be higher 
with higher sink velocity. The wing third mode velocity variance 
rises to a peak after impact and then decreases with a slow rate 
showing little transience oscillations . The initial oscillations 
are absent for velocity variances of track modes also. These, 
however, indicate the presence of oscillation in the later phase 
of landing run. The initial impact oscillations are subdued in the 
acceleration variance showing only one or two peaks . 

5.2.2 Heave Pitch Model with Rigid fuselage 

For this model, response statistics are presented for ten 
coordinates -fuselage rigid body heave, pitch, nose and main gear 
heave degrees of freedom and first three generalised coordinates 
for the one wing and the track. 

5 . 2 . 2 . 1 Taxi Run 

The mean system response are presented in Fig. 5. 11. At 
transition from smooth to rough surface, the vehicle experiences a 
shock at t=0. This is indicated by strong initial mean for all the 
response coordinates considered. Influence of the mean track 
pattern is apparent in the vehicle lumped mass response. Increased 
forward velocity induces larger mean response magnitudes . 
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run. Key same as Fig. 5. 7 
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Key: Aircraft sink velocity, 0,6 m/s , 0.9 m/s 
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The input track frequency sensed by che vehicle is modified 
by ins — o r wa r d velocity and are higher for higher velocicies . The 
mean response show low frequency components superimposed wich 
higher frequency components. The lower frequency elements indicate 
sensitivity to the vehicle forward velocity. However, the higher 
frequency components undergo little modification with the forward 
velocity. 

Mean responses of the wing normal coordinates show distinct 
single dominant frequency at transient stage which persist for a 
long duration. The response behaviour of track normal coordinates 
show low amplitude vibration at track input frequency which change 
with the change in the vehicle speed. 

The response variances, presented in Fig. 5. 12, show large 
number of peaks, gradually decreasing with time. The distinctness 
of the pattern is reduced from displacement to velocity, co 
acceleration. Higher forward velocity induces stronger response. 
The higher frequency component does not show significant variation 
with change in forward velocity. The acceleration response attains 
steady state in a very short duration. 

5 . 2 . 2 . 2 Taheof f Run 

Fig. 5. 13 shows the mean response during takeoff. These have 
initial low values, subsequently rise with the increasing speed of 
the vehicle. Influence of mean track profile is apparent on the 
aircraft lumped masses showing low frequency components 
superimposed with small amplitude high frequency components . The 
low frequency component indicates sensitivity to the magnitude of 
the forward motion. The mean displacement and acceleration 
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decrease prior to takeoff. Tlie aircraft pitch, rotations are found 
to be oscillatory in pattern. The track normal coordinates also 
have oscillatory pattern that gradually increase in magnitude with 
passage of time. The longest takeoff time induces highest track 
mean response . 

Fig. 5. 14 represents the response variances during takeoff. As 
the vehicle starts from rest, increasing its speed gradually, it 
takes some time for the track input to become effective. When the 
effective speed level is attained, variances are found to increase 
with the increase in speed showing an oscillatory pattern over an 
increasing mean value. The high frequency components reduce with 
the time derivative of the response for the continuous member, 
being lower in velocity and acceleration . 

5 . 2 . 2 . 3 Landing Run 

Mean responses are presented in Fig. 5. 15. The effect of impact 
at touch down, which increases with increase in the sink velocity, 
is reveled in the mean response. With the progress of landing run, 
initial energy gets dissipated and response -shows dependence only 
on the track roughness inputs. Pitching motion of the sprung mass 
is found less dependent on the ‘sink velocity. The nose gear mean 
displacement and velocity show much slower dissipation of impact 
energy. Main gear response magnitudes are less compared to the 
nose gear mean characteristics . Track normal coordinate 
displacement means are found to be more sensitive to the change of 
sink velocity compared to the velocity and acceleration response. 
The track acceleration shows a high value at touch down impact 
which steeply falls with time to a very low value. 
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Response variances during landing run are presented in 
Fig. 5. 16. The variance pattern is oscillatory in general. The 
effect of touch down impact is also seen here with initial high 
peak response which subsequently decreases with reduction in the 
vehicle speed. The magnitude of the peak is higher for higher sink 
velocity. Effect of sink velocities is seen to be present in the 
wing higher normal modes. Track normal coordinate displacement 
variances show pronounced fluctuations compared to velocity and 
acceleration variance . Displacement variance magnitude has greater 
dependence on the sink velocity. 

The wing has coupled bending - torsion behaviour and a space 
distribution of the bending - torsion displacement ratios in the 
first three modes are shown in Fig. 5. 17. The figure reveals strong 
coupling between the two motions in the second and third modes 
compared to the first one. The first mode response is mainly 
dominated by the bending deformation towards the tip for the 
system parameters selected in the study. 

5.2.3 Heave-Pitch Model with Fuselage flexible in Bending 

For this model, first five flexural modes of the fuselage are 
considered out of which first two represent rigid body heave and 
pitch degrees of freedom. The other three represent elastic 
bending modes. Besides these, two lumped masses heave degrees of 
freedom, wing and track first three generalised coordinates are 
considered as in the previous rigid model . The response 
characteristics of the thirteen generalised coordinates are 
presented for the taxi, takeoff and landing run. Results have been 
presented for one wing only. 
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Time (s) Time (s) Time (s) 

Fig. 5. 11. Heave-Pitch Model (Rigid fuselage)-Mean response 
in taxi run. Key same as Fig. 5. 5 
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Fig. 5.1 1 (continued). Heave-Pitch 
response in taxi run. Key same 
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Fig. 5. 13. Heave-Pitch Model (Rigid fuseiage)-Mean response 
in takeoff run. Key same as Fig. 5. 7 
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response in takeoff run. Key same as Fig. 5. 7 
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Fig. 5. 14. Heave-Pitch Model (Rigid fuselage)-Response 
variance in takeoff run. Key same as Fig. 5. 7 
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Fig. 5.1 4. (continued). Heave-Pitch Model (Rigid fuselage)- 
Response variance in takeoff run. Key same as Fig. 5. 7 
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Fig. 5. 1 5(continued). Heave-Pitch Model (Rigid fuselage) -Mean 
response in landing run. Key same as Fig. 5. 9 
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■ Time ( s ) Time (s) Time (s) 

Fig.5.16. Heave-Pitch Model (Rigid fuselage)-Response 
variance in landing run. Key same as Fig. 5. 9 
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Fig. 5. 1 6(continued). 
Response variance 








5 . 2 . 3 . 1 Taxi Run 

The mean system response are presented in Fig. 5. 18. Rigid 
body vertical motion and bending deformations of the fuselage show 
the influence of mean track profile. Track input frequency is 
modified by the vehicle forward velocity. Flexural response have 
low values at the end of taxi run compared to rigid body vertical 
heave motion. Pitch motion exhibits high frequency components 
initially/ which diminished with passage of time. Wheel response 
also have more or less a similar pattern. Acceleration of the nose 
gear shows high initial value which gradually decreases to reach a 
low steady state value. Initial oscillations do not show much 
variations with the change in vehicle velocity. This points to 
their association with structural modes rather than the external 
excitation. 

Mean response of wing's normal coordinates shows less 
sensitivity to track influence in the lower modes. Higher mode 
response, however, indicates a increase in frequency with increase 
in the vehicle speed. Track normal coordinate displacement and 
velocity means show an increase in magnitude with increasing level 
of track mean height experienced by the wheels . Track 
accelerations have higher values at the onset of motion. This may 
be due to the shock on entry of the vehicle over the track. 

Response variance during taxi run are presented in Fig. 5. 19. 
Variances are oscillatory in the initial stage with a high peak at 
the initiation of motion. These subside to reach steady 
asymptotic levels with passage of time. The acceleration variances 
achieve steady value within a shorter time compared to 
displacement and velocity. Track mode characteristics show low 



frequency variation in the steady state part. Higher vehicle 
velocity is fo un d to increase the response magnitudes . The 
response of the elastic members indicate substantial contributions 
from the lower modes also. 

5. 2. 3. 2 Takeoff Run 

Fig. 5. 20 presents mean response during takeoff. Fuselage 
generalised coordinates have dominant low frequency components 
with small amplitude high frequency components superimposed over 
it. The response amplitudes show an initial low value which 
gradually builds up with increase in speed and track mean level. 
Higher forward acceleration generates higher response. Differences 
are observed at take off instant depending on its location in the 
low frequency cycle. Influence of mean track profile is apparent 
in the air frame as well as main wheel response. Nose wheel 
response during takeoff is found to be smaller compared to the 
main wheel. The airframe and wing lower modes show marginal 
reduction in amplitude with progress in time. Track normal 
coordinates mean response also show a dominant low frequency 
component with gradually increasing amplitude . The velocity 
response have superimposed high frequency small amplitude 
components while the acceleration response have some low frequency 
components also. 

Fig. 5. 21 represents the response variances during takeoff. 
When the effective speed level for the aircraft is achieved, 
variances are found to increase steadily with the increase in 
vehicle speed. Oscillatory pattern over an increasing mean value 
is noticeable in most of the cases. A gradual increase in response 
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variance is also observed in the track normal coordinates. The 
oscillatory pattern over increasing mean level and amplitude is 
more pronounced for the track modes. 

5 . 2 . 3 . 3 . Landing Run 

Mean responses are presented in Fig. 5. 22. The effect of 
impact at touchdown, which increases with increase in sink 
velocity, is revealed in the mean response. As in the previous 
model, with the progress of the landing run, impact energy gets 
dissipated and response shows dependence only on the track 
roughness input. Higher bending modes of the fuselage show 
increasing mean of the response pattern with progress of the 
landing run. Main wheel displacement response has faster 
dissipation of the impact energy compared to the nose wheel. 

Wing generalised coordinate mean responses have much higher 
peaks at impact phase compared to the response amplitudes in the 
later part of the landing run. Response in the wing's higher modes 
are found to damp out much faster than the lower modes for 
velocity and acceleration. Mean responses of the track generalised 
coordinate also reveal the effect of landing impact. The magnitude 
of the track response at touchdown instant is higher than the 
constant velocity and accelerated runs . 

Response variances for the landing run of the aircraft are 
presented in Fig. 5. 23. Oscillatory patterns are present for all 
the displacements and velocity response in the touchdown phase. 
Magnitudes of the^first peaks seem to be higher with higher sink 
velocity. Response variances of the wing's second and third normal 
coordinates show single dominant peak at touch down phase which 
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steeply falls to lower values during the landing run. Track normal 
coordinates displacement variance show oscillations with higher 
frequency variation while velocity and acceleration responses have 
dominant low frequency components . 

5.2.4 Heave-Pitch-Roll Model (Flexible fuselage in bending with 
rigid roll) 

In this model, response of five generalised coordinates of 
fuselage flexure , rigid body roll of fuselage, three unsprung 
masses vertical bounce motion, first three normal coordinates of 
left and right wing coupled bending-torsion modes and first three 
flexural modes of the track have been considered. The first two 
fuselage generalised coordinates give the rigid body heave and 
pitch motion. Plots have been presented for eighteen generalised 
system modes in taxi, takeoff and landing runs of the aircraft. 

5 . 2 . 4 . 1 Taxi Run 

Mean system response are presented in Fig. 5. 24. Fuselage 
rigid body bounce motion shows increase in displacement with 
increasing height of the track mean profile in presence of rising 
gradient selected for this study. Mean velocity and acceleration 
response show high transience set up at the onset of vehicle 
forward motion which tends to diminish with passage of time. The 
response magnitude is higher at the transient stage compared to 
value at the later stage. Pitch response reveals the presence of 
high frequency components which persists in transient and the 
steady state part. Change of vehicle forward velocity shows minor 
effect in case of pitch velocity and acceleration. 
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Fig. 5. 19. Heave-Pitch Model (Flexible fuselage in bending)- 
Resonse variance in taxi run. 

Key same as Fig. 5. 18 
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Mean response of third, fourth and fifth normal coordinates 
of fuselage elastic modes show the influence of mean track profile 
as in the case of rigid body heave motion. However, response 
magnitudes are smaller compared to the rigid body heave motion. 
Fuselage rigid body rolling pattern shows steady state condition 
at constant velocity run. Effect of change of vehicle forward 
velocity has more pronounced effect in the roll compared to pitch 
response. Wheel mean displacements are seen to be influenced by 
track mean variation and vehicle forward speed. Nose wheel 
velocity and acceleration mean has low steady state response. The 
two main wheel response characteristics show a difference, 
indicating the effect of fuselage rolling. Similar is the case for 
the left and right wing's generalised coordinate's mean response 
which differ in magnitude as well as phase. 

Track normal coordinates mean response grow in magnitude with 
increasing track mean level and vehicle speed. Higher vehicle 
forward velocity is seen to induce higher response in the track. 
Acceleration mean of the track show large value on entry of 
vehicle over the track. 

Variance of displacement, velocity and acceleration of the 
system generalised coordinates are presented in Fig. 5. 25. Variance 
of displacement and velocity of the fuselage rigid body heave and 
pitch motion has oscillatory pattern in the early stages which 
tend to an asymptotic value with time. Pitch acceleration variance 
has initial high value which decreases to a low asymptotic value 
with low amplitude fluctuation. Variance of fuselage third, fourth 
and fifth normal coordinates, rigid body rolling and wheel heave 
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degrees of freedom reveal similar characteristics . Response 
variance of left and right wing normal coordinates show 
fluctuation in early stage which subside to reach asymptotic 
values. Steady value is reached much faster in case of velocity 
and acceleration response. 

Track normal coordinates response variance indicate low 
frequency variation in later part of taxi runs which become more 
pronounced in velocity response at higher forward velocity. 

5. 2. 4. 2 Takeoff Run 

Fig. 5. 26 presents mean response during takeoff of the 
aircraft. Rigid body heave motion show gradual building up of the 
response with increase in vehicle forward speed. Longer takeoff 
run is seen to induce higher displacement magnitude for the heave 
motion. Pitch acceleration during takeoff is seen to decrease 
after initial transience as aircraft proceeds to takeoff. Pitch 
acceleration is found to have less dependence on the change of 
vehicle forward acceleration. Mean response of fuselage elastic 
bending reveals high frequency components in the initial stage of 
motion. High frequency oscillation show low sensitivity to vehicle 
forward acceleration. Rigid body roll of the airplane indicates 
low value at the initiation of motion. Higher vehicle forward 
acceleration induce higher roll acceleration in the airplane. Nose 
wheel displacement and velocity response show high frequency 
components which tends to diminish as the takeoff run progresses. 
Response of two main wheels are found different. At takeoff 
instant mean displacement of right wheel is seen to have larger 


194 



magnitude compared to that of left wheel. Mean response of left 
and right wing's first and second normal coordinate show more 
persistent transience compared to that in the third normal 
coordinate . 

Track normal coordinate mean response show growing up of 
amplitude with increasing speed. High frequency components present 
in the response indicates less sensitivity to vehicle forward 
acceleration . 

Response variances are presented in Fig. 5. 27. Variance of 
fuselage rigid body heave motion shows gradual building up of the 
amplitude with increase m vehicle speed. Acceleration variance 
show very low value which steeply increases during takeoff time. 
Pitch displacement and velocity have a fluctuating pattern which 
is not present in acceleration response. Fuselage normal 
coordinates response variance in elastic bending modes follow a 
pattern similar to the rigid body heave mode. Roll variance also 
indicates a rising tendency with the increase in vehicle forward 
velocity. Oscillations are more pronounced at higher vehicle 
forward acceleration. Nose and main wheel response variance follow 
a pattern similar to the fuselage response. High frequency 
components are more prominent in wheels acceleration response 
variance . 

Left and right wing response variance magnitudes increase 
with vehicle forward speed. Acceleration variance of right wing's 
second normal coordinate, however, indicate the presence of 
asymptotic level at takeoff instant. Track normal coordinate 
response variance reveal low frequency oscillation which is more 
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apparent in velocity and acceleration. Acceleration response of 
third normal coordinate is found to have a steep rise at takeoff 
instant . 

5. 2. 4. 3 Landing Run 

Mean system response in landing run have been presented in 
Fig. 5. 28. Fuselage rigid body modes have high frequency 
oscillation during the impact phase. Response is gradually damped 
out with the slowing of vehicle forward motion. Higher sink 
velocity induces higher heave response. Pitch response also show 
the effect of landing impact. Fuselage third and fifth 
coordinate's mean displacements reveal low frequency oscillation 
in the later stage of landing run. Mean displacement and velocity 
of rigid body rolling has second peak higher than the first one 
after touch down whereas the peaks progressively reduce in the 
acceleration response. Wheels response also reveal the effect of 
impact at landing touch down which gradually decays and the 
response become dependent on track mean input . 

Wing mean response show dominant values at impact phase with 
a cluster of closely spaced peaks . Effect of landing impact dies 
out soon showing very low response in the later phase of the 
landing run. Mean response of the track generalized coordinates 
capture the high frequency oscillation after the vehicle causes 
impact on the track. Response magnitude decreases with the slowing 
of the vehicle forward motion. 

Response variance during landing run are presented in 
Fig. 5. 29. Out of the fuselage first five normal coordinates, rigid 
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body heave motion show less number of peaks compared to other- 
normal coordinate response variance. Acceleration variances of 
second, third, fourth and fifth generalised coordinates are found 
to decay faster compared to displacement and velocity. 
Displacement and velocity variance of rigid body rolling has first 
peak of much higher magnitude compared to subsequent peaks . Two 
main wheels response variances follow the pattern similar to the 
heave response variance. Nose gear displacement and velocity show 
slower dissipation of impact energy compared to the two main 
wheels. Nose gear acceleration variance at touch down instant has 
higher value compared to main wheels. Wing's response are 
oscillatory in the early stage. The oscillatory nature reduces 
from displacement to velocity to acceleration. Wing acceleration 
variance in the first mode indicate two peaks, first peak being 
much higher compared to second one wheras acceleration variance in 
second and third mode indicates single dominant peak. 

Track normal coordinates displacement variance show 
oscillatory pattern with high frequency components in the impact 
phase. Velocity and acceleration show low frequency variation in 
the later stage of the landing run. 

5.2.5 Heave-Pitch-Roll Model (Flexible fuselage in bending and 
torsion ) 

First five flexural modes and first four torsional modes of 
the airframe have been considered in this three point input model 
out of which first two in flexure represent rigid body heave and 
pitch while first one in torsion is for the rigid rolling of the 
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Ng.b.24(continued;. Heave -Pitch-Roll Model (Flexible 
fuselage in bending with rigid roll)-Mean response in taxi 
run. Key same as Fig. 5. 18 
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Fig. 5.25. Heave-Pitch-Roll Model (Flexible fuselage in 
bending with rigid roll) -Response variance in taxi run 
Key same as Fig. 5. 18 











Fig. 5. 26. Heave-Pitch-Roll Model (Flexible fuselage in 
bending with rigid roll)-Mean response in takeoff run 
Key same as Fig. 5. 20 
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Fig. 5. 27 (continued). Heave-F 
fuselage in bending with ri 
takeoff run. Key same as F 
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Fig. 5. 28. H eave - P i tch - Roil Model (Flexible fuselage in 
bending- with rigid roll)-Mean response in landing run. 
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airframe. First three coupled bending - torsional modes of both 
the wings and first three flexural modes of the, track have been 
taken into account. Along with the heave response of the three 
unsprung masses, results have been obtained for twentyone 
generalised coordinates for this model. 

5 . 2 . 5 . 1 Taxi Run 

Mean response are presented in Fig. 5. 30. Rigid body heave, 
pitch and roll displacement of the airframe show influence of 
input process mean shape. Initial transience are more prominent in 
velocity and acceleration mean. Out of the two angular rigid body 
modes, rolling motion is more affected by the change in vehicle 
forward velocity. Fuselage first three elastic modes in bending 
and torsion reveal the presence of high frequency components in 
mean response for the initial phase which gradually diminish to 
reach steady state patterns . Single dominant frequency in the 
steady state part which gets modified by the change in vehicle 
forward velocity can be observed. High frequency oscillations in 
the transient stage, however, are less sensitive to the change in 
vehicle velocity. 

Unsprung masses response exhibit the similarity in pattern 
with the rigid body heave motion. Nose wheel acceleration mean 
shows high transience initially and takes considerable time to 
reach steady state. Two main wheels mean response characteristic 
are not found exactly identical. This signifies the effect of 
rolling of the vehicle on account of transverse variation of the 
track level. 
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Left wing mean displacements exhibit clear steady state 
pattern in constant velocity run which is subdued in mean velocity 
and acceleration. The velocity and acceleration mean are highly 
effected by the transience set up at the onset of motion. Right 
wing mean response behaviour is found to follow the trend of the 
left wing response except that the difference in phase in later 
phase of taxi run. 

Track normal coordinates displacement and velocity mean show 
an increase in magnitude with increasing level of track mean 
height experienced by the wheels . Track accleration mean has 
strong transience at the onset of motion and low steady state 
value. This may be due to shock on entry of the vehicle over the 
track . 

Response variances of the system generalised coordinates in 
taxi run are presented in Fig. 5. 31. The displacement variances 
have an oscillatory pattern with high peak at the onset of the 
vehicle motion. Subsequently the response subsides to reach a 
steady asymptotic value. The variances of the velocity response 
follow a trend similar to the displacements. The acceleration 
variances indicate dominant values at the initiation of motion. 
The asymptotic value is reached within shorter time than that for 
displacement and velocity, showing progressively smaller peaks. 
The asymptotic values, however, are very low compared to the 
initial value. 



5. 2. 5. 2 Take off Run 

Fig. 5. 3 2 presents the mean system response during takeoff. 
The vehicle is assumed to start from rest, runs with gradually 
increasing velocity till liftoff. Initial portion of the mean 
response shows low transience. Mean response magnitudes are seen 
to build up gradually with the increase in vehicle forward speed 
and track mean level which has a rising gradient. However, as the 
aircraft proceeds to takeoff, mean response shows decrease in 
amplitude. Rigid body motion in roll seems to be more pronounced 
with incensing velocity. Response behaviour of the right main 
wheel is different from the left main wheel showing strong 
interaction of rolling in the accelerating runs . 

Mean response of the generalised coordinates of the left wing 
are found to be different from the right wing. Less transience is 
observed in the right wing mean response. However, phase 
difference is noticed between the response of the two wings. The 
track normal coordinate mean response shows initial low values 
which subsequently increase as the vehicle vibration picks up with 
increase in speed imposing more dynamic loads . 

Response variances of the system generalised coordinates are 
presented in Fig. 5. 33. The variances for the vehicle response show 
very low values in the initial stage. Variances are seen to 
increase with forward acceleration of the vehicle, but at takeoff 
instances, differences in magnitudes are not appreciable. 

Variances of the first three track normal coordinates have 
low frequency components in their response characteristics, which 
are more prominent in velocity and acceleration. For the slower- 
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takeoff acceleration, the track mode response variance seems to 
have higher magnitudes . 

5. 2 . 5. 3 Landing Run 

Fig. 5. 34 presents the mean system response during landing 
run. The effect of impact at touchdown, which increase with the 
increase in sink velocity , is revealed in the mean response of all 
the degrees of freedom considered m the model . With the progress 
of landing run, impact energy gets dissipated and response shows 
dependence only on the track roughness input. The effect of the 
mean input excitation is seen in fuselage elastic bending modes 
in displacement. However, in the airframe elastic torsional mode 
mean response the effect of input mean excitation is subdued. 
Torsional modes of the fuselage have faster dissipation of energy 
compared to rigid body mode in rolling. 

Nose and main wheel response mean are found to have similar 
pattern. However, impact energy dampens quickly in case of right 
gear mean response. Response mean of wings generalised coordinates 
have very closely spaced peaks caused by the landing impact . These 
dampen out as the vehicle runs at decreasing speed. Wing's higher 
mode response mean shows slight dependence on the track roughness 
in the latter phase of landing run. Mean response of the track 
generalised coordinates also reveal the effect of landing impact. 
The magnitude of the track response at touch down instant is 
higher than for constant velocity and takeoff run. 

Response variances for the landing run of the aircraft are 
presented in Fig. 5.35. Oscillatory patterns are observed for all 
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the rigid body and flexible modes of the airplane during touchdown 
impact. Rigid body pitch and roll response variances are quite 
similar in pattern. Number of peaks in roll response variance are 
smaller compared to pitch variance. The effect of impact persists 
for slightly longer duration for the roll acceleration variance. 
Nose wheel response variance shows slower dissipation of impact 
energy than the main wheels response . 

Wing response variance pattern resemble the vehicle response 
variance. Two peaks are seen in acceleration variance of both left 
and right wing's first normal coordinates. Amplitude of the first 
peak is much higher than the second peak. Track normal 
coordinate's displacement variance show oscillations with high 
frequency components while velocity and acceleration response have 
dominant low frequency components . 

5.2.6 Comparison of Some Linear Models Behaviour 

In response analysis, selection of model is very important. 
The purpose of the present section is to present a comparative 
study of response behaviour of different vehicle models so that a 
approximate guidelines fqr the selection of model can be framed. 
Four different vehicle models have been considered. These are (1) 
Model 1- Heave -Pitch rigid aircraft model (2) Model 2 - Heave - 
Pitch flexible aircraft model with bending flexibility of the 
fuselage. (3) Model 3 - Heave - Pitch - Roll model with fuselage 
bending and rigid roll. (4) Model 4 - Heave - Pitch - Roll model 
including fuselage bending and torsional flexibility. 
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bending and torsion)-Mear 
Key same same as Fig. 5," 























C\J oo -tf 

CO CM T- 

CO ''t OJ 

LO o LO 

CO OJ T- 

o o 

o o o 

odd 

T- 1 T - o 

o o d 

^/ T xn 

^/ C J 9-0I 

Wz^e-OT 

> s /z m p- 0 l 



O) CO CO 

OJ 00 Tj- 

CO OJ 

h- 00 CD 

XT CO 03 

odd 

t- d d 

odd 

oi o 

oi-r-'o 

01 

z s /z J i - 01 

gff/z^-OI 

Z S/ ux 01 

z s ! 5^9-01 


in o to 

T“ T““ C) 


CD <M 

odd 


CO Tt CM 
CO oJ t-* 


CD CD CO 

odd 


xr CO CO 

• mm 

C\J 'T- O 


Fig. 5.31. Heave-Pitch-Roll Model (Flexible fuselage i 
bending and torsion)-Response variance in taxi run. 
Key same as Fig. 5. 18 
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Fig. 5, 33. Heave-Pitch-Roll Model (Flexible fuselage in 
bending and torsion)-Response variance in takeoff run. 
Key same as Fig. 5. 20 
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Fig.5.3 3 (continued). Heave-Pitch-Roll Model (Flexible 
fuselage in bending and torsion)-Response variance i 
takeoff run. Key same as Fig. 5. 20 
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eoff run. Key same as Fig. 5. 20 











bending and torsion) -Mean response in landing run. 
Key same as Fig. 5. 22 
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Fig. 5. 35. Heave-Pitch-Roll Model (Flexible fuselage in 
bending and torsion) -Response variance in landing run. 
Key same as Fig. 5. 22 
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fiuselage in bending and torsion) -Response variance in 
landing run. Key same as Fig. 5. 22 









Displacement means and variances of rigid body heave, pitch, 
roll with first two elastic bending modes of the fuselage, nose 
and one main wheel response, first two coupled bending - torsion 
modes of one wing and track first normal coordinate have been 
presented for comparison. Taxi, takeoff and landing run have been 
studied with vehicle speed 8 0 Km/h in taxi, forward acceleration 
2.0 m/s 2 in takeoff and sink velocity 1.2 m/s in landing runs. 

Fig. 5. 3 6 presents the displacement means and variances in 
constant velocity runs . Response mean of the rigid body modes of 
the vehicle shows steady state oscillation representing track mean 
profile for all the models considered. Rigid vehicle model (Model 
1) predicts higher response mean and variances for the heave and 
pitch degree of freedom. In flexible structures, a part of the 
dynamic energy is being utilised to cause elastic deformation. 
As the vehicle body becomes more flexible, displacements in rigid 
body modes are further reduced. Roll response mean and variance of 
the Model 4 shows slightly lower value compared to that of Model 3 
which is assumed torsionally rigid. High frequency ocsillation in 
the nose gear response is more prominent in Model 1 with higher 
amplitudes than other models considered. In case of main wheel 
response, slightly higher value is obtained in three point input 
models (Model 3 and 4) . Fuselage elastic bending modes response 
mean and variance has less value in Model 4 . Wing mean and 
variance in two point input models (Model 1 and 2) show 
predominant high frequency oscillation which are found to be less 
in three point input models (Model 3 and 4) . Influence of mean 
track profile is more apparent in the higher mode mean response of 
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wing in Model 1 and 2. Track first normal coordinate's mean 

displacement reveals high value of the transient response f or 
Model 1. In Model 3 and 4, the effect of roll motion is seen to 
reduce the response mean and variance of track deflections. 

Fig. 5. 37 presents the comparison of response mean and 
variance in takeoff run. Rigid model shows high value of response 
mean and variance for the heave and pitch degree of freedom. 

Variances are found oscillatory in the rigid model. Torsional 
flexibility is found more effective to reduce the response 

magnitude in case of rigid body roll motions induced in three 

point input model (Model 4) during takeoff run. Nose gear 
responses have more high frequency components in Model 1 compared 
to other models. Main wheel response for all the models do not 
reveal much difference in various models except the variance is 
found to fluctuate about an increasing mean in Model 1. Fuselage 
elastic bending response in flexible three point input models 
(Model 3 and 4) retains similar characteristics. Flexible wing's 
first and second normal coordinate mean response in Model 3 and 4 
show appreciable deviation from that of Model 1 and 2 . Track 
response builds up gradually with increase in speed in flexible 
models (Model 2, 3 and 4) 

Mean displacement and variances in landing run of the 
aircraft are compared in Fig. 5. 38. Heave and pitch motion of the 
aircraft has slightly higher peaks in the impact phase for Model 1 
and 2 . Mean and variance of the rigid body rolling of the vehicle 
in Model 3 and 4 have similar behaviour at touchdown impact with 
little difference in the amplitude. In case of fuselage second 
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elastic bending mode mean displacement, peak developed at touch 
down instant has higher magnitude in Model 3 and 4. However, 
variance shows peaks of low magnitude and faster dissipation of 
impact energy compared to two point input models (Model 1 and 2) . 
Flexible wing's normal coordinate mean and response variance in 
rigid model show much significant effect of the landing impact. In 
the latter phase of landing run, response behaviour in all the 
models are guite similar. Track behaviour are similar in all the 
models with higher peaks developed due to impact of rigid vehicle 
on the track. However, some departure in response behaviour is 
noticed during the progress of landing run. 

Comparative study shows that, in general, rigid vehicle 
response magnitudes are higher than the flexible models although 
evolution of mean and variances with time predicts similar nature. 
Inclusion of elastic torsional mode of the vehicle shows minor 
influence on the response behaviour as compared to elastic bending 
mode . 

5.3 NONLINEAR MODEL RESPONSE STATISTICS 

Results of the heave model of aircraft with tapered flexible 
wing are obtained with the system parameters mentioned earlier. 
Shock absorber characteristics of the nonlinear suspension are 
assumed as 

Initial air pressure (P a Q) : 3 ‘ ^ x ^ N/m , pneumatic area 
and volume (A aQ , V aQ ) : 0.65 x 10 -2 m 2 and 0.2518 x 10~ 2 m 3 , 
hydraulic and orifice area (A h , A n ) * 80.5 x 10 4 m 2 , 3.044 x 
10" 4 m 2 , coefficient of orifice discharge 0 . 6 , fluid 
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Fig. 5. 37. Comparison of mean 
takeoff run. 

Key same as Fig. 5. 36 
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density (p) : 850 kg/m , polytropic index (n) : 1 . 4 . 

Response mean and variance are obtained for the system generalised 
coordinates in taxi, takeoff and landing run by two methods - 
instantaneous linear approach and Monte Carlo simulation. 

5.3.1 Taxi Run 

A constant vehicle forward velocity of 22.2 m/s has been 
assumed. Displacement, velocity and acceleration mean and 
variances are presented in Figs. 5.39, 5.40 and 5.41. 

The displacement means of the two lumped masses show the 
characteristics of the track mean profile. The figures reveal 
similarity in response behaviour obtained by the two methods . The 
responses show the presence of high frequency components in the 
simulated results which are not always mimicked by the response of 
linearised model. However, the underlying low frequency pattern 
and the response in later stage of taxi run in the two cases are 
similar. The agreement in the variance is better than the mean 
between the two methods . The agreement improves with the progress 
of the motion. This is seen, in general, in all the three 
responses- displacement, velocity and acceleration. 

5.3.2 Takeoff Run 

2 

Takeoff corresponds to an accelerated run at 2.0 m/s 
starting from rest till liftoff velocity of 59.72 m/s is attained. 
Mean and variance of displacement are shown in Fig. 5. 42. The high- 
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frequency component in the wing normal m odes for the simuiated 
response means are not fully reflected by the instantaneously 
linear model. Otherwise the two methods have a close match for the 
mean an< ^ the variance for the inodes considered 

Displacement variances of sprung and unsprung mass show 
gradual building up of amplitude associated with the low frequency 
component. This may be due to effectiveness of ground input with 
increasing velocity. However, prior to takeoff, displacement 
variance indicate a decrease in magnitude. Wing response in both 

the models show growth of the amplitude at much faster rate in the 
initial period. 


Velocity response mean and variance for takeoff run are shown 
in Fig. 5. 43. The lumped mass mean velocities show close agreement 
between the two methods except for a shift in the positions of the 
peaks . The continuous member mean velocities reveal high frequency 
components which are not reflected in the linearised model. 
However, the two solutions follow similar general trend. This is 
true for the variance response also. 

Fig. 5. 44 presents the mean and variances of acceleration 
response. As for displacement and velocity response, the 
linearised model does not reflect the high frequency pattern of 
the simulated results completely. The low frequency behaviour in 
the two cases are similar. 


5.3.3 Landing Run 


Results have been obtained for the landing run of the 



aircraft with sink velocity 1.2 m/s, glide velocity 6 0 m/s and 
constant deceleration of 1.5 m/s 2 . The results have been presented 
for the first 20 seconds of the run. 

Mean and variances of displacement are shown in Fig. 5. 45. The 
linearised sprung and unsprung masses response again miss the high 
frequency oscillation during transience caused by the landing 
impact. The other response components also exhibit a similar 
general behaviour. During the later phase of the landing run, the 
agreement improves between the two approaches. 

Displacement response variance have a better agreement 
between the two methods compared to the mean characteristics. 
Variance of the unsprung mass shows a slower dissipation of impact 
energy compared to variance of the sprung mass. Response variance 
of other degrees of freedom presented show the formation of peak 
caused by landing impact. In most of the cases, a single peak is 
seen. Peaks are found to have some differences in magnitude and 
position for the two models . 

Mean and variance of velocity and acceleration response are 
presented in Figs. 5.46 and 5.47 respectively. The velocity and 
acceleration, in general, follow the pattern of the displacement 
response. However, dependence of mean response on the track 
roughness input reduces from displacement to velocity to 
acceleration. Results of the instantaneously linear model are 
found to be close to the simulated results when the vehicle loses 
a significant amount of its forward velocity. 
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Fig. 5. 39. Nonlinear Aircraft Model-Mean and variance of 
displacement in taxi run. 
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Fig. 5. 41. Nonlinear Aircraft Model 
acceleration in taxi run. 

Key same as Fig. 5. 39 








































CHAPTER 6 


CONCLUSIONS 

The conclusion based on the results obtained are being 
resented in this chapter. Following the preceeding chapter, the 
onclusions are also placed in three sections. 

.. 1 NONUNIFORM BEAM VIBRATION 

An analytical approach has been outlined to study the free 
'■ibrational characteristics of a nonuniform beam. The beam mass 
ind stiffness can be of general nature with distributed damping 
)roportional to mass distributions. The approach is applicable for 
ill types of boundary conditions. Some conclusions based on the 
study are : - 

(1) Natural frequency of a nonuniform beam is sensitive to the 
taper ratio and magnitude as well as location of the 
concentrated mass loading present in the beam. 

(2) The polynomial shape function derived may be useful to 
analyse nonuniform flexible structurs modelled by finite beam 
elements. The use of the analytical nonuniform beam shape 
function instead of the traditional beam shape function is 
expected to bring out rapid convergence with less number of 
elements . 
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6.2 RESPONSE STATISTICS OF LINEAR VEHICLE MODELS 

An approach has been outlined to analyse the flexible 
vehicle-track coupled dynamic behaviour for track irregularity in 
variable velocity run. The method uses linear suspension 
behaviour. A general pattern of vehicle forward motion and spatial 
nonhomogenity of track roughness can be included in the analysis. 
Some major conclusions drawn from the study are:- 

(1) Mean response of the vehicle and track is primarily 
influenced by the mean track profile inducing mean input 
excitation . 

(2) Periodicity present in the mean track profile is reflected in 
the vehicle response. Response frequency is, however, modified 
by the vehicle forward velocity. 

(3) Variance characteristics of the vehicle response is 
influenced by track second order statistics. ' 

(4) Vehicle response is greatly dependent on the track roughness 
and vehicle speed. It is important to maintain track in proper 
condition to keep roughness low in high speed operation. This 
is expected to increase undercarriage and vehicle fatigue 
life . 

(5) System response is found to be higher during touch down 
impact compared to taxi and takeoff runs. 

(6) Flexible wings have strong transience set up at landing touch 
down and the system response suggests that inclusion of wing 
flexibility in the aircraft model is important for the 
analysis of possible metal fatigue damage as well as for 
catastropic failure. Flexible wing has appreciable response in 
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taxi and takeoff runs also. Thus inclusion of slender members 
in the vehicle model is important from their design and 
performance considerations. 

■ 7 ) increased speed of the vehicle is found to yield large 
response in the track bed. However, response magnitudes are 
smaller for the system parameters selected in the study. This 
low level of track deformation may not affect vehicle response 
significantly but its inclusion in system model may be 
important for track fatigue life and rehabilitation. 

(8) Unsymmetrical distribution of roughness across the track 
causes a significant system response not only in the roll mode 
but also in the heave and pitch modes due to strong coupling 

in the system. 

(9) The modelling predicts higher response magnitudes for rigid 
vehicles as compared to flexible vehicles. 

(10) Bending flexibility of aircraft type of vehicle is important 
for accurate description of the response behaviour. The low 
tread makes the torsional flexibility less significant for 
the prediction of vehicle response. 

(11) Consideration of rigid vehicle model may be adequate for 
preliminary design of undercarriages and vehicle structural 

components . 

6.3 NONLINEAR AIRCRAFT MODEL 

An instantaneously linearised approach for evaluation of the 

response statistics for nonlinear vehicle model has been presented 

for nonstationary track inputs. The method differs from other 
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linearisation techniques in the way that instead Qf f . nding 

equivalent set of linear parameters at one time, it introduces 
' instantaneous ' spring and damping coefficients evaluated from 
the response of the original nonlinear system at discrete time 
instants very close to each other. The method combines numerical 
technique with the linear random vibration theory. Accuracy of the 
method is tested with results obtained by Monte Carlo simulation 
of the nonlinear model. Some conclusions from the study are 

(1) The simulated response characteristics indicate high 
frequency components which are not always reflected by the 
linearised model . The simulation characteristics may contain 
digital noise introduced due to numerical integration and 
small sample size of 100 . The noise is not present in the 
linearised solution leading to the above difference. This 
argument gains support from the presence of high frequencies 
in some linearised characteristics and absent in others. 

(2) Response characteristics show a good match between the two 
methods in the later stage of the vehicle run. 

(3) The variance matching is better compared to the mean between 
the two methods . 


(4) The present method can handle general combinations of 
nonlinearities found in vehicle suspension. It does not 
require a priori the probability density function of the 
response which is an essential feature of some other existing 
linearisation techniques . 

(5) The method proves cost effective as compared to digital 
simulation in terms of computer processing time and storage 



space. Ratio of CPU time for the present linearised model and 
Monte Carlo simulation results per hundred samples is 
approximately 1 : 60. 

(6) The instantaneous linearisation approach has a potential for 
application to general nonlinear dynamics problems . 

6.4 SUGGESTIONS FOR FUTURE WORKS 

The present study has application to vehicle-track dynamics 

problems. The following scope exists for its extension: - 

(1) Vehicle forward motion and initial conditions may be treated 
as random variable. 

(2) Random variation of track subgrade properties with 

discontinuous jointed pavement can be included. 

(3) Effect of shear deformation, rotary inertia and warping 
deformation can be included in vibration of nonuniform beam 
element . 

(4) Vibration control studies can be carried out with the 
proposed models. 

(5) A proceedure for fatigue life estimation can be developed for 
undercarriage and flexible vehicle components . 

(6) Instantaneous linearised approach for nonlinear models can be 
extended to multi wheeled input models including tyre 
nonlinearities . 
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